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SECTION 1 


INTRODUCTION 

Over the past several years, Acurex developed the HET catalytic 
combustor computer model as an aid in the solution of problems related to 
catalytic combustion applications for various systems such as boilers and 
gas turbines. The HET model differs from previous models in that it can 
handle the high temperature effects of catalytic combustion where bed radia- 
tive heat transfer and "flame type" phenomena are important by including gas 
phase reactions and wall radiation terms. For a balance between accuracy 
and economy of operation, the code utilizes a quasi-one-dimensional model. 

In the present project, the, HET code is used as a basis for examining dif- 
ferent concepts utilizing catalytic combustion for Stirling engine heater 
heads. 

Current heater head designs for Stirling engines depend on convective 
heat transfer from the combustion mixture to heat the working fluid. This 
approach results in a nonuniform circumferential and axial heat transfer 
distribution on the heater tubes. Since the tube temperature distribution 
is somewhat dependent on the local heat transfer rate, a temperature dif- 
ference will then also exist on the heater tubes. The maximum allowable 
material temperature then dictates the maximum allowable heat flux. Tech- 
niques which result in uniform tube temperatures or heat fluxes will maxi- 
mize the net amount of energy that can be transferred to the working fluid. 

Two advanced concepts that have the potential for more uniform heat 
fluxes are catalyzed heater head and radiative energy transfer from catalytic 
reactors to the working fluid heat exchanger. The catalytically coated 
heater head can react the fuel and transfer the heat of reaction directly 
into the engine working fluid in a single component. With this concept. 



heating would be more spatially uniform allowing a higher average head tem- 
perature and thereby increased engine efficiency. Also, combustion tempera- 
ture and extent of reaction could be tailored to achieve very low exhaust 
emission levels. Finally, construction would be simplified by incorporating 
the combustion surface into the head design. In the radiatively cooled con- 
cept, heat would be transferred radiatively as well as convectively from the 
catalytic reactor to the working fluid heat exchanger. The advantage of this 
concept is that catalyst surface temperatures could be maintained considerably 
above the working fluid temperature resulting in efficient combustion. With 
this concept, catalyst and heat exchanger temperatures could be more easily 
tailored for combustion and heat exchanger efficiency as well as emissions. 

Section 2 will present the basic HET model and the modifications and 
additions introduced by various heater head concepts. Section 3 will describe 
the code input instructions. 

This work was performed in support of the Department of Energy, HEAT 
Engine Highway Vehicle Systems Project under DOE/NASA Interagency Agreement 
DE-AI01-77CS51040. 
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SECTION 2 


THEORY 

This section briefly summarizes the theoretical foundations upon which 
the calculation methods of the Stirling Engine HET Code program are based. An 
attempt has been made to present these fundamentals in a fashion most appro- 
priate to a user unfamiliar with the code. To this end, and in the interest 
of brevity, explanations tend to be nonrigorous and perhaps somewhat heuristic. 
In discussing each type of problem treated by the program, the pertinent govern- 
ing equations are developed and it is shown that these relations are sufficient 
to determine the unknown quantities for that particular problem. Since all 
problems related to the Stirling engine are essentially solved by a modifica- 
tion of the basic HET code developed earlier by Acurex (Reference 2-1), this 

is then 
Stirling 
and trans- 


section will begin with a description of the basic HET code. This 
followed by discussions on how the basic code is modified for each 
concept. Calculations and correlations for radiation view factors 
port properties are then discussed at the end of this section. 
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2.1 BASIC HET CODE 

The basic model describes the reactions of a prevaporized and premixed 
fuel and air stream inside a catalytically coated channel. The following 
phenomena are treated: 

• Homogeneous gas phase chemical reactions 

t Heterogeneous surface chemical reactions 

• Radial heat and mass transport 

• Axial bed conduction 

• Radiative heat exchange between walls 
The following assumptions are made: 

• Negligible gas phase axial heat and mass transport diffusion 

• Wall fluxes can be treated by transfer coefficients which 
directly relate fluxes to the driving forces 

• Radial conductive heat transfer is neglected 

The resultant quasi one dimensional governing equations are: 


Species 


dY. 

m = AW. - C 
ds 1 w w. 


( 1 ) 


where 


m = mass flowrate/channel 

= mass fraction of species i in the bulk gas 
s = axial distance 

A = cross-sectional area of flow channel 

= perimeter of flow channel 

W^ = chemical production rate (gas phase) of species i 

J, , = flux of species i at the wall 

W1 


4 



Energy 


where 


m dh _ p _ 

m 3“ = - C Q 
ds w^w 


w 


enthalpy in the bulk gas phase 
heat flux at the wall 


( 2 ) 


Species Flux at the Wall 

The flux of species i at the wall is determined by the rate of trans- 
port of i to the wall and the rate of reaction of i at the wall. These are 
given by: 




( 3 ) 


and 







) 


where W = rate of destruction of species i at the wall 
'"i 

Y = mass fraction of species i at the wall 

C = dimensionless mass transfer coefficient 

(C = — , k = mass transfer coefficient 
m. V c 

V = free stream velocity) 


( 4 ) 


Energy Flow at the Wall 

The wall heat flux is governed by convective heat transfer, chemical 
reaction, radiation exchange between walls and conduction in the solid phase. 
These are given by: 
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( 5 ) 


A 

„ _ „ s K w 

^ ds'^ 


m — Y ^ 

%-i<^n f’-V*lT;f'’w'''i-V> 

^ n 1 1 


where 


q = radiative heat loss 

^r 

A = cross-sectional area of the solid bed 
s 

k = conductivity of the solid bed material 
s 

T = wall temperature 

w 

= dimensionless heat transfer coefficient 
h. 

(Cr = » h = heat transfer coefficient 

P p = density 

Cp = heat capacity) 

F = enthalpy of edge gas at wall temperature (defined as 
w 


h =yy,-h ). 

w Zrf 1 w . 


h^ = enthalpy of species i at wall temperature 

Chemical reaction rate constants are calculated by Arrhenius type rate 
expressions of the form: 


kp = aT exp(-E/RT) 


where kp = reaction rate constant 
E = activation energy 


a,b = constants 
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Radiation heat exchange is calculated using a view factor approach 
where qp., the wall radiative heat flux at station j is given by 



(B) 


where K is the channel segment view factor 
k denotes all other stations except j 
r^r2 denote upstream and downstream reservoirs 


Derivation for the view factors for different geometries are given in 
Section 2.6. Transfer coefficient calculation/correlations used in the code 
are given in Section 2.7. 

The boundary conditions required to solve the system of equations are: 


= = “ • ''i “ ’'o, • T ■ ds 


w 


= 0 


s=0 


dT , 
_w| 

ds 


= 0 


s=L 


The above equations are solved simultaneously using a finite difference 
implicit numerical scheme. The results will provide both temperature and 
species distributions in the gas phase and at the catalyst surface. In the 
case when the catalyst wall temperature distributions are known, one can also 
over specify the problem by inputting this information into the program. The 
program will then automatically bypass the catalyst wall energy equation. 


These correspond to insulated edges. If heat transfer occurs from the ends 
of the channel to the surroundings, it can easily be taken into account by 
equating the conductive heat flux to heat loss from the ends, say by con- 
vection and/or radiation. 
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2.2 EXTERNALLY COOLED CATALYTIC CYLINDRICAL REACTOR/FLAT PLATE 



This problem is essentially the same as the basic HET problem except 
that the heat is also removed externally from the tube. Heat conduction must 
therefore occur transversely as well as axially in the tube walls. To in- 
corporate this effect, a two-dimensional heat conduction equation was imple- 
mented into the formulation. This is given by: 

^ + ^=0 (9) 

9s ar'^ 


where r = radial distance 

r-| = inside tube radius 
V 2 = outside tube radius 

The above equation is not exact since it does not account for area variations 
due to changes in radius. It is, however, a good approximation where r^ » 
(r 2 ■ 'f'^) which is typical for cylindrical catalytic reactors. 

Four boundary conditions are needed for this equation. For the radial 
direction, these are given by heat exchange between tube and coolant and tube 
and fuel . 
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( 10 ) 


h.(T„ 


r=r. 


dT 


r=r. 


dT , 

q = q - C k 

s dr ' 


r=r. 


where ha, = external heat transfer coefficient 


w 


inside tube perimeter 


T = average external coolant temperature 


( 11 ) 


The heat transfer occurring from the outside of the tube to the coolant 
is represented by an external heat transfer coefficient h^. Heat transfer may 
also occur by radiation. If heat transfer by radiation from the surface of 
the tube is important, then h^ can include this effect by linearizing the tem- 
perature term of radiation and assuming h^ is an effective heat transfer coef- 
ficient representing both convective and radiative heat transfer. 

q^, the wall heat flux, is the same as that defined in the basic model 
(Equation (5)), except the axial conduction term is replaced by a conduction 
term at the surface of the wall. The boundary conditions for the axial direc- 
tion remain the same as in the basic model with either insulated ends (no heat 
transfer at the ends of the tube) or conductive heat transfer at the ends of 
the tube balanced by radiation and convection. 

The input conditions for this problem are similar to the typical 
monolithic reactor problem except that the external heat transfer coefficient 
and coolant temperature must be inputted. The values of these input parameters, 
however, need not be constant. 
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2.3 FLOW ALONG FLAT PLATE 



For the flat plate problem, the bulk gas temperature and composi- 
tion remain constant and the boundary layer is growing as the gas streams 
travel down the tube. The basic HET model equations are easily adapted 
to this problem and much simplified. Equations (1) and (2) are not 
necessary since the bulk composition and temperatures remain constant 
Wall compositions and temperatures can be solved by equating Equations 
(3) and (4), and (5) and (6). Thus, the flat plate problem, in a way, 
is a subset problem of the HET code. Transfer coefficients, however, 
must be modified accordingly and also there is no radiation exchange 
between walls (only radiation heat loss to the ambient). 
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2.4 COANNULAR TUBE RADIATIVELY COOLED CONCEPT 


T, air ♦ fuel 

In coolant 


Cataj^5t_coatea wan 




Coolant tube 1 

VT 

’’2 '"l 

T 





In this problem, coolant and fuel flow in concentric tubes with 
coolant in the inner tube and fuel in the annular region. The inner 
walls of the outer tube are catalytically coated. Heat is transferred 
to the inner tube by convection and radiation. The following assump- 
tions are made: 

t The inner coolant tube has a thickness Ar which offers some 
resistance to heat transfer 

• Axial conduction in the inner tube is neglected 

0 Radiation from the inner tube wall to the catalyst coated 

wall is neglected since the temperature of the inner tube 
wall is much lower than the catalyst wall 

For this problem, the basic equations of the HET model can be 
adapted with a few minor modifications. Equation (2) has to be modi- 
fied to include heat exchange with the inner tube in addition to the 
catalyst wall. This is given by: 


m 


dh 

ds 


= - C q 
w^w 


- 


r=r. 


( 12 ) 


where 


C = outer perimeter of inner coolant tube 
”2 

h. = heat transfer coefficient from bulk gas to wall surface 
of coolant tube 

T' = wall temperature of coolant tube (function of both 
^ axial and radial directions) 
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To solve for the wall and coolant temperature distributions, the 
following heat balances on the coolant tube are needed: 


- ■r; 


dT‘ 

) + k—^ 
^ dr 


r=r. 


r=r. 


+ q' 

^r 


= 0 


( 13 ) 




r=r. 


dT 


) - k — ^ 
> ^ dr 


= 0 


r=r. 


(14) 




r=r 


- T,) = 0 


1 


(15) 


where 


C 

^r 


■h 


H 


conductivity of the coolant tube 

inner perimeter of coolant tube 

heat transfer by radiation from catalyst wall to 
coolant tube wall 

heat transfer coefficient of coolant 

mass flowrate and heat capacity of coolant, respectively 


If the axial condition is neglected, the radial conduction heat 
flow in the coolant tube is given by: 


-kA 


dT' 

w 


dr 




r=r. 




r=r^ 




) (16) 


r=r. 


where A = tube surface area for heat transfer 
L = tube length 
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Combining Equations (13), (14), and (16), and substituting into 
Equation (15) gives: 


dT 


ds 


H 


- U,C 

1 Wi 




(17) 


boundary condition 


= T|^ at S = 0 


where U-j = overall heat transfer coefficient 


1 


Ai Ai^n(r2/r^) _ ^ 

2nkL h 


^2^t 


H 


Ao = surface area at r = r^ 


A^ = surface area at r = r 


1 


If we assume that the tube wall is thin, can be simplified to 


Ui = 


1 


i_+ i_ 

^ k h, 


(19) 


The second term on the right hand side of Equation (12) can also 
be expressed as functions of T and only to give: 


. dh 
d? 


C q - C 
w^w w. 


Ai A, 


( 20 ) 


Equations (17) and (20) can be incorporated into the basic HET 
equations and solved simultaneously. Once and T are solved, the 
coolant wall temperature profile can readily be derived from Equa- 
tions (13) and (14). 
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2.5 MONOLITHIC REACTOR RADIATING TO UPSTREAM AND 
DOWNSTREAM HEAT EXCHANGERS 


Coolant 



The problem of a monolith placed between two banks of coolant 
tubes is rather straightforward in that no modification is necessary 
to the basic HET equations. The problem of heat transfer to the cool- 
ant is decoupled from the solution to the HET equations since heat 
transfer to the tube banks does not affect the HET equations. Heat 
transfer to the tube banks occurs by two mechanisms. The first is by 
convective flow of the fuel over the tube banks. The second is by radi- 
ation from the monolith walls to the tube banks. Back radiation from 
the tube banks to the monolith walls can be neglected since the tube 
walls are at significantly lower temperatures compared to the catalyst 
walls. The equation describing heat transfer to the tube banks is simi- 
lar to Equations (13), (14), (15), and (17) except that convective heat 
transfer occurs by crossflow over the tubes and radiation heat transfer 
occurs from the whole monolith to the tube bank. The key equation to 
be solved is: 



UnC 
1 w. 



* (T - T„) 


boundary condition 



at X = 0 


( 21 ) 
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where x = axial distance along coolant tube 
T. + T 

T = 2 ^-j ■ temperature of fuel air mixture 

approaching coolant tube 

T = temperature of fuel air mixture 
after passing over coolant tube 

q = radiation heat flux from the monolith to the tube bank, 
assumed uniformly distributed along the tube length 

If the inlet fuel/air temperature (temperature after passing over 
the first tube bank) is known, the equations of the basic HET model can 
be solved without modification to give the wall and gas temperatures in 
the monolith. The monolith wall temperatures can be combined with suit- 
able view factors to give the radiation heat transfer flow to the tube 
banks. The exit temperature of the gas stream from the monolith is the 
temperature of approach to the second tube bank. Equation (21) can be 
solved analytically by assuming the radiation heat flux term (q^) is 
uniform over the length of the tube to give the hydrogen temperature as 
a function of distance x: 


where 


T = — 

'h 


Bid - e 




" V 


-A^x 


(22) 


U,C’ 
1 w 


1 


H 


Uiqr 


C. 


w 


+ U,T 


1 


I'/VVo 

” Ph, 


The wall temperature distribution of the coolant tubes can easily be 
solved once T^^ is known. 


15 



2.6 CALCULATION OF RADIATION VIEW FACTORS 


Calculation of radiation view factors for different geometries is 
relatively straightforward. Radiation heat transfer between walls of a 
cylindrical tube is given by (Reference 2-2): 

09 


f>» 

id 



Surface 2 


Y 




22 


= 0 


.5 |x 



(23) 


Radiation heat transfer between concentric tubes is given by (Refer- 
ence 2-2): 


R2 



A = R^ + - 1 

0 

B = - R^ + 1 
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F 


23 


0.5(1 - F2^ - F 22 ) 


(26) 


In the case of the monolith radiating to a bank of tubes, we can assume 
that the face of the monolith is very close to the tube bank for maximum heat 
transfer. The two faces can therefore be treated as semi-infinite plates with 
the view factor given as (Reference 2-2): 


II II ! / 

r 

h 


! II II ! 

• — w ► 


H = h/w 


Fi 2 = /l + - H 


(27) 
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To correct for the fact that one of the plates is really a bank of 

tubes of diameter D^. with center to center distance D^, the actual view factor 

t _ c 

from monolith to tube bank (Reference 2-3): 

% = y ^}2 

y = -0.059Z^ + 0.049Z +1.01 1 1 Z < 2 

(29) 

= 0.0145Z^ - 0.224Z +1.26 2 < Z 


where 


Z 
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2.7 TRANSFER COEFFICIENTS 


Values for gas properties and transfer coefficients are based on 
correlations and theoretical equations. The fuel/air side consists of ana- 
lytical equations and correlations already built into the HET code. They 
are more refined than the accuracy of the present problem requires but since 
they are already there, little effort is made to convert them to simpler cor- 
relations. Of the equations and correlations used for the fuel /air side, the 
following will be described in more detail. 

The fuel/air side consists of a complex mixture of reactants and pro- 
ducts that changes in concentration with distance down the channel . Trans- 
port properties such as diffusion, viscosity, and thermal conductivity must 
therefore make use of correlations based on mixtures. The development of 
expressions for these is discussed in detail elsewhere (Reference 2-8). A 
brief summary of this development, and the resulting expressions, is presented 
here. 

Correlations for binary diffusion coefficients are given in terms of 
a bifurcation approximation discussed in detail in Reference 2-8. This is of 
the form: 



(30) 


where D is a reference diffusion coefficient and F^ are diffusion factors. 

In the present code, oxygen is chosen as the reference species with Fq^ = 1.0. 
From the molecular theory of gases and liquids (Reference 2-9): 


0 • 2.628 X 10"^ 


T(T/H 

;::2 ;rrTT7 


2 

(cm /sec) 


(31) 



with T in °K, P in atmospheres, and collision cross section, o in A. For 

o 

O 2 as the reference species, a is equal to 3.467 A. Using the data from 
Reference 2-9, the integral expression for transport properties is approxi- 
mated by: 

1.07[T/(e/k)]‘°*^^^ (32) 

* J 

where the maximum energy of attraction function, e/k, for O 2 is 106.7 and 
thus: 

D = 0.172 X 10'^ j1.659 /p (cm^/gec) (33) 

For system viscosity, use is made of the correlation suggested by 
Buddenberg and Wilke (Reference 2-10) and endorsed by Hirschfelder et al . 
(Reference 2-9), namely 


mix 






1 X, 


+ 1.385 


RTy, 




PM, - f 

J 

J7i 


where is the mole fraction for component i 

M^. is the molecular weight of component i 
is the viscosity of pure component i 


(34) 


By introducing the bifurcation relations, taking from Reference 2-9 the rela- 
tion for pure component viscosity 


= 


6At. 


p.D. . 

n 


(35) 
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where is a ratio of collision integrals. Assuming A+^ = 1.12 (actually 
varies from 1.10 to 1.14 in the temperature range of interest), and adjusting 
1.385 to 1.344 for simplification, there is obtained 


- ^2 

^mix " 1.3448^M 


(36) 


where 8-| = 


The thermal conductivity of a mixture of polyatomic molecules can be 
written (Reference 2-10) as: 


+ k. . (37) 

mix mono-mix int 

where k is the mixture thermal conductivity computed neglecting all 

internal energy modes and is the contribution of the internal energy 
modes of the molecules to the mixture conductivity. An approximate relation- 
ship for k . developed by Mason and Saxena (Reference 2-11) is: 

^ mono-mix r j 


mono-mix 


I 


X,k. 


mono 


RT 

Xi + 1.385 


X.. 




j?*i 


1.1 


(38) 
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where kj is the thermal conductivity of the pure species i, neglecting 


'mono 

all internal degrees of freedom. 
(Reference 2-9) as: 


Expressions for ki can be written 

'mono 


k. 

^mono 


15 _R. 
4 M, 


The expression for k then becomes: 

mono-mix 


(39) 


mono-mix 



r X 


J. 15 R 
F. 4 M 


1.344 


(40) 


To complete the definition of k, the expression for the internal energy mode 

contribution to thermal conductivity must be added to k 

mono-mix 

From Equations (59) and (77) of Reference 2-12, the following relation 
can be obtained: 


M. 


*^int 


= 1- 


M {%. ~ 2 


I 


(41) 


^int 




‘■“Xi T (s, ■ I 

e,F, 


(42) 


The value of k used is then: 


_ pD 


Cp.0.29R63 


(43) 


where 


1 ^i 1 ^ ^ Pi 

^3 " M 2 F7 ^p " M 5! ~F^ 


(44) 
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The mass transfer coefficient is developed from the heat transfer 
coefficient by 


= C,(Le.) 


2/3 


(45) 


where Le. = Lewis number. 

On the coolant side, the gases used are either hydrogen, helium or air. 
Since there are no reactions in the coolant, gas properties and transfer 
coefficients can be treated as a one component gas stream and simple expres- 
sions can be derived. 

Gas viscosity can be represented by the Hirschfelder, Curtiss and Bird 
equation (Reference 2-9): 


]i = 26.69 



(46) 


where 


y = viscosity (yp) 

T = temperature (K) 

M = molecular weight 

O 

a = hard sphere diameter (A) 
= collision integral 


, - ( 1-161 \ 0.525 , 2.162 

V " \ 0.1501,^^ exp(0.773T|^)^ exp(2.438T|^) 


(47) 



(48) 


£ _ characteristic energy of interaction between molecules 
k Boltzmann constant 
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The values for a and e/k for the coolant used are given below: 


Coolant 

ojAl 

e/k (K) 

Hel ium 

2.551 

10.22 

Hydrogen 

2.827 

59.7 

Air 

3.711 

78.6 


Thermal conductivity can be calculated by the relation (Reference 2-4) 



If we know the thermal conductivities of the gas at two temperatures, 
the value of n can be obtained from Equation (49). 

Heat capacity of the coolant is correlated by the general expression: 

Cp = a + bT + CT'^ (50) 

T = temperature (K) 

Cp = heat capacity (cal/g mole °C) 

Values for a, b, c for each gas can be readily obtained from any 
handbook. 
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SECTION 3 


CODE INPUT INSTRUCTIONS 

This section describes the punched card input needed to activate the 
problem options described in Section 2, A number of comments are included to 
help guide the user in setting up input card decks. 

3.1 CARD INPUT DECK 

Card 1 (Format 20A4) 

1-80 Title of run 

Card 2 (Format 713) 


1 

- 3 

IS 

Number of species 

4 

- 6 

NL 

Number of grid points 

7 

- 9 

NIT 

Number of iterations allowed 

10 

- 12 

ILOSS 

Heat loss option, 


0 - adiabatic reactor option 

1 - traverse heat loss option 

2 - assigned catalyst temperature option 

13-15 IGEOM Reactor geometry option, 

0 - monolith or cylindrical reactor option 

1 - flat plat option 

2 - coannular tube radiatively cooled 

reactor option 

3 - monolith reactor radiating to up- and 

downstream heat exchangers problem 
option 


26 



16 - 18 

ICOEF 

Input external heat transfer coefficient 
option 

>0 - invoke this option and bypass the 
built-in correlation 

19 - 21 

KR7 

Diagnostic output option. 



>3 - invoke diagnostic output option 

Card 3 


(Format 5F10.0) 

1 - 10 

DTUB 

Effective cell diameter (cm) 

11-20 

VF 

Void fraction of monolith 

21 - 30 

DIAMS 

Diameter/height of reactor (cm) 

31 - 40 

AK 

Thermal conductivity of reactor (cal/s/cm-K) 

41 - 50 

EMIV 

Emissivity of reactor 

Card 4 


(Format 3F10.0) 

1 - 10 

AM 

Mass flowrate (gm/s) for IGEOM f 1. Mass 
flux (gm/cm^-s) for IGEOM = 1. 

11-20 

TI 

Inlet temperature (K) 

21 - 30 

P 

Pressure (atm) 

Card 5 


(Format 2F10.0) 

1 - 10 

TRES(l) 

Upstream reservoir temperature (K) 

11-20 

TRES(2) 

Downstream reservoir temperature (K) 

Card Set 


Spatial grid parameters 

Card(s) 1 

(Format 8E10.5) 

1 - 10 

S(N),N=1,NL 

Selected grid locations (cm) 

11-20 
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Card Set 7 


First guessed/assigned wall temperature 
profile 


Card(s) 1 


(Format 8E10.5) 

1 - 10 

11 - 20 

TL(N),N=1,NL-1 

Wall temperature at node point N (K) 

Card Set 8 
Card(s) 1 


(Format 4012) 

1 - 2 

IKIN(N),N=1,NL-1 

If the node is a noncatalytic surface, 

3-4 


enter zero. If it is a catalytic sur- 

• mm 


face, enter 1 

Card Set 9 


Parameters for ILOSS >0 and IGEOM >0 
options 

Card 1 


(Format 5F10.5) 

1 - 10 

VL 

Coolant initial velocity (cm/s) 

11 - 20 

AMH 

Coolant flowrate (gm/s) 

21 - 30 

TTI 

Coolant inlet temperature (K) 

31 - 40 

PT 

Coolant pressure (atm) 

41 - 50 
Card 2 

WT 

Coolant molecular weight (gm/gmole) 

1 - 10 

AKT 

Thermal conductivity (for coolant tube) 
(cal/s-cm-K) 

11-20 

DR 

Coolant tube thickness (cm) 

21 - 30 

DT 

Coolant tube diameter (cm). Leave 
blank if IGEOM >2. 

41 - 50 

CC 

Distance between monolith and tube bank 


(cm). Leave blank unless IGEOM = 3. 
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Card 3 


(Format 4F10.5). (See Equation 49.) 


1 - 10 

TKl 

Lower reference temperature for thermal 
conductivity of the coolant (K) 

11 - 20 

TK2 

Upper reference temperature for thermal 
conductivity of the coolant (K) 

21 - 30 

GKTl 

Coolant thermal conductivity corresponding 
to TKl 

31 - 40 

GKT2 

Coolant thermal conductivity corresponding 
to TK2 

Card 4 


(Format 5F10.5) 

1 - 10 

SG 

1 Parameters for calculation of coolant 

11 - 20 

EK 

1 viscosity (see Equations 46 to 48) 

21 - 30 

Cl 

Coefficients for coolant specific heat 

31 - 40 

C2 

^ correlation (see Equation 50) 

41 - 50 

C3 

) 

[ C-| = a, C 2 = b, C^ = c 

Card Set 

JO 

Initial species concentrations, wall species 
concentration guesses, and species diffusion 
factors 

Card(s) 1 

1 to IS 

(Format 2A4, 2X, 2E10.3, F10.4) 

1 - 8 

mm 

Species name utilized in kinetic reaction 
input. (These names must also be compatible 
with the names on the thermochemical input 
data file.) See Section 3.2 for a descrip- 
tion of thermochemical data. 

11 - 20 

ALPHF 

Initial species concentration in relative 


mole concentrations 
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21 - 30 ALPHE Estimate of final wall species concentra- 

tions in relative mole concentrations 

31 - 40 BPA Bifurcation diffusion factors, 1/F^, for 

each species. (See Appendix A for the 
definition of these factors) 

Card Set 11 Kinetic data 

Reactions and their associated rates are given in the form: 



u. B. 

m 



and 


Jin P. Jin p.-Jln K 

g m . e 


where y are the stoichiometric coefficients on species B^ , P^ is 

the partial pressure of species i, and Kn is the corresponding 

^m 

equilibrium constant. The forward reaction rate coefficient, , 
has the Arrhenius form 

m 


The backward reaction rate coefficient is taken to be equal to the 
forward reaction rate coefficient divided by the equilibrium constant. 



Card 1 


(Format 213) 


1 - 3 

MT 

Total number of reactions input 
ing both gas phase and surface 

(includ- 

reactions) 

4 - 6 

MGAS 

Number of gas phase reactions 
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Card(s) 2, 

4,...,2MT 

(Format 10(A4,1X), 3E10.4) 

1 Name of reactants (must be compatible with 
Card Set 10 

1 - 

4 

NA(1) j 

j 

21 - 

24 

NA{5) J 

1 

26 - 

29 

NB(1) ^ 

1 Name of products (must be compatible with 
j, Card Set 10 

46 - 

49 

1 

NB(2) 

1 

51 - 

60 

FKF(M) 

Pre-exponential factor of the mth reaction 
(in mole, cm, sec units) 

61 - 

70 

EXK(M) 

Temperature exponent 

71 - 

80 

EAK(M) 

Activation energy (kcal/gmole-K) 

Card(s) 3, 

5 2MT+1 

(Format 5F5.0, T41 , 5F5.0) 

1 

1 - 

5 

AMU(l) > 

• • * 


1 

1 

1 R 

y. associated with NA(i) 
1 * 

21 - 

25 

AMU(5) J 

1 

41 - 

45 

BMU(l) 'j 

1 

. . . 


1 

i 

P 

)■ associated with NB(i) 

61 - 

65 

BMU(5) J 

1 


3.2 THERI-10CHEMICAL INPUT DATA FORMAT 

For the program operation, thermochemical data (i.e., specific heat, 
enthalpy and entropy) for all species included in Card Set 10 must be present 
on mass storage devices designated by unit 11. The order of this data on the 
storage device is not important. However, it is important that the species 
names input on Card Set 10 match identically with the species names on the 
thermochemical data files and the kinetic reaction input data. Card Set 11. 
The thermochemical data files can be full libraries of data, from which only 
a limited number of species will be selected for each problem of interest. 

2 

Curve fits of data are in the form Cp = F^ + F^T + F^/T and cover a 
lower and upper temperature range. Cp units are cal/mole/K. 
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Card 1 


(Format lOX, 2A4, 2X, 2F10.0) 

n - 18 

NAMA 

Species name utilized in initial concen- 
tration and reaction input 

21 - 20 

MW 

Molecular weight of species 

Card 2 


(Format 6E9.6, 6X, F6.0, 6X, A4) 
Low temperature curve fit. 

1 - 9 

»^1 

Heat of formation at 298K (cal/mole) 

10 - 18 

•^2 

Change in enthalpy from 298K to 3000K 
(cal/mole) 

19 - 27 


Coefficient in above Cp expression 

28 - 36 

^4 

Coefficient in above Cp expression 

36 - 45 


Coefficient in above Cp expression 

46 - 54 

*^6 

Entropy constant at 3000K (cal/mole) 

61 - 66 

TU 

Upper temperature limit of low tempera- 
ture range curve fit (K) 

73 - 76 

NAMA 

Name of species used in reaction set 

Card 3 


(Format 6E9.6, 6X, F6.0, 6X, A4) 
High temperature curve fit 

1 - 9 

^1 

Heat of formation at 298K (cal/mole) 

10 - 18 

^2 

Change in enthalpy from 298K to 3000K 
(cal/mole) 

19 - 27 


Coefficient in above Cp expression 

28 - 36 

^4, 

Coefficient in above Cp expression 

36 - 45 

>^5 

Coefficient in above Cp expression 

46 - 54 

^6 

Coefficient in above Cp expression 

61 - 66 

TU 

Upper temperature limit of low tem- 
perature range curve fit (K) 

73 - 76 

NAMA 

Name of species used in reaction set 


3 ? 



SECTION 4 


SAMPLE PROBLEMS 

Several sample problems were run with the Stirling Engine Combustor 
Code to show the variety of problems that can be solved. In this section, 
each problem is briefly described, input card data is shown and select out- 
put results are presented. 
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SAMPLE PROBLEM NO. 1 
ADIABATIC CYLINDRICAL REACTOR (CERAMIC) 


Description: This problem represents a fuel /air mixture flow in a catalytic ceramic tube with no 

heat loss to the environment. This problem is also the same as catalytic combustion 
in an adiabatic honeycomb reactor. 


Fuel : 


Propane/air 

Equivalence ratio 
Mixture flowrate 
Preheat temperature 
Pressure 



0.2 

0.0394 g/s 
BOOK ^ 

3 X 105 Pa 


Geometry : 

Cell diameter 
Reactor length 
Void fraction 


Dfub “ 0.32 cm 
l = 20 cm 
Vp = 0.87 



iiAMPt-t 

PROBLEM A 

- uNContEn 

CAT alt ST MOfgOLITH 

P/A EOUIV 

«ATIO=:0.2 


5 **2 

10 0 

0 0 






n. 1^875 

0.67 

0.17 

u.ouus 

U.ft 




0,0394 

000.0 

3.0 






«no.o 

-1.0 







u.o 

0.0b 

0.10 

0.1b 

0.2 

0.25 

0.3 

0.35 

n.4 

0.4b 

0.5 

0.7b 

1.0 

1.5 

1.7b 

2.0 

ii.b 

3 . 

3.b 

4 . 

4.5 

5. 

5.b 

6. 

b.b 

7. 

7.b 

a. 

tt.b 

9. 

9.5 

10. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200, 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200, 

1200. 

1200. 

1200. 

1200. 

1200. 

1 1 1 

11111 

111111 

1 1 1 1 1 

1111 

1111 

1111111 

1111 

N2 

.70334 

.77666 

1.0347b 





02 

.20033 

.16529 

1.10045 





C3H8 

.00033 

2.0 E- 

031,2756 





H20 

1.0 E 

-10 .03306 

.84406’ 





C02 

1.0 E 

-10 .02479 

1.27446 





2 1 








C3H6 02 


C3H8 02 

C02 

9.1 

E+9 

40.61 

.25 1. 

04 



-.75 - 

3.964.0 

3,0 


C5H8 


02 

H20 C02 


1.1 

E+09 0. 

17.6 

1.0 




-5.0 4 

. 0 3.0 





SAMPLE PRORLtM A - I'rjforM KCi CATALYST MON^>HTH r/A kHUIV HATI0=0,2 


* integral PARAMETERS ♦ 


IS ( NUMBER OF SPECIES I = 5 
NL ( GRID POINTS ) . = 

Nil ( number of overall ITERATIONS ) =: 1 
ILOSS ( HEAT LOSS OPTION ) = n 
IGEOM ( PROBLEM GEOMETRY OPTION ) = n 
ICCEFF ( HEAT TRANSFER INPUT OPTION ) = 0 
KR7 ( PRINT OPTION ) = 0 


* NONINTEGRAL PARAMETERS ♦ 


♦ ** REACTOR [JAlA ♦♦♦ 


CO 

CT» 


monolith honeycomb catalytic 

NUMBER OF TUBES = 
VOID fraction =. 
niAMETER OF combustor r 
EFFECTIVE CELL niAMETER = 


COMBUSTOR 
1.00 
• B70 

.17 (CM) 
.lb (CM) 


INLET TEMPERATURE 

PRESSURE 

MASS FLOW RATE 

CONDUCTIVITY 

EMISSIVI lY 


BOO. 00 (K) 

3,00 (aTM) 

.394-001 (G/S) 
.500-003 ( CAL/S-CM-K ) 

.ao 


♦** reservoir temperatures 

UPSTREAM = 800,00 

OOUNSTHEAm = 1.0000 


GRID POINTS CHOSEN (S) 


,0000 

.5000 

4,500 

«^,500 

,5000-001 , inuo + OOIJ 

.7500 l.UOU 

5.000 5,500 

10,00 

. 1500 
1,500 
fa, 000 

.2000 
1 .750 
6,500 

.2500 

2,000 

/.ooo 

,3000 

2.500 

7.500 

.3500 

5.000 

6.000 

• •♦oon 
3.500 
H.500 

.4500 

4.000 

9,000 

INITIAL GUFSSEO WALL TFMpe;h ATURFS ( TL ) 
1200, 1200, 1200. 1200. 
1200. 1200. 1200. 1200. 

1 ? n 0 , 1 2 0 (1 . 12 0 0, 1200, 

1200, 
1200, 
1200 . 

1200, 

1 2 n . 

1 2 Ij 0 . 

1200, 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

1200. 

12U0. 

1200. 

1200. 





♦ SPKCltS IMiri/vL MOLf H P AC T I ONS , P 1 KST GUCi>SF_StAHn niFFUSlUTNl FftCTOHS ♦ 


sflcie: 

ALPF 

^2 

.763+000 

C2 

,2064000 

C3M6 

.833-002 

H20 

.loo-ooy 

C02 

.100-009 


ALPt DIFFUSION F ACTOH 

.777+000 1.0347 
.lGt3 + 00O 1.1004 
.200-002 1,2756 
.331-001 ,6441 
.246-001 1,2745 



• THCm»inCHE"1ISTHY PATA » 


CURVE fit of data IN FORM CP=HU+nC*T+RU/ ( T*T I < C AL/ « HOLE*K ) 


HF(CAL/r*OL> 

H(CAL/MnU 

Rrt 

PC 

KO 

S(CAL/MOL/K) 


TU(K) 

name: 

N2 

79UUO+000 

28,000 

,23995+005 

,61391+001 

.16415-002 

,2<T/27+005 

.64548+002 

300. 

1000, 

N2 

79UUO+000 

.22168+005 

,63049+001 

.23357-003 

7^:006+006 

.63774+002 

lOOO. 

2500. 

N2 

02 

00000 

32.000 

,25541+005 

.67044+001 

.16795-002 

-.19H19+005 

,68908+002 

300. 

1000. 

02 

OUOOO 

.234504005 

,80725+001 

.50513-003 

-.21419+006 

,67976+002 

1000. 

2500. 

02 

H20 

57798+005 

18.000 

.31716+005 

.66290+001 

.29575-002 

.276H3+00b 

.69022+002 

300. 

1000. 

H20 

57798+005 

,30303+005 

,99696+001 

.12774-002 

-.19426+007 

,68463+002 

1000. 

25U0, 

H20 

C02 

94054+005 

44.000 

,41438+005 

,95776+001 

,36020-002 

-.16469+006 

,81997+002 

300. 

inoo. 

C02 

94U54+005 

.36573+005 

,13972+002 

,36449-003 

-.l'»‘*12 + 007 

,79867+002 

1000. 

2500. 

C02 

C3H8 

,248204.005 

44,000 

•10866+006 

.32040+002 

,11928-001 

-.20472+007 

,15940+003 

300. 

2500, 

C3H8 

,24820+005 

,10888+006 

,32090+002 

•11928-001 

-.20472+007 

,15940+003 

2500. 

6000, 

C3H8 



* KiMfTiC HLACTlON DATA ♦ 


TOTAL UUlillFH OF KtACTIOfJS ^ 

AuMbLK OF OAS PHASF KLACTiUrjS 1 


h 

TYPF 


RTACTION 




PRE FXP FAClOP 

TEMP ExP 

ACTIVATION 








(rtOLL-CM-SI 


(KCAL/MOLE) 

1 

GAS 

5.0 0? 

♦ 1.0 C3HB — st.O 

H20 

+ 3,0 

C02 

.mno'foiu 

• 000 

H0.6100 

2 

SURF 

5.0 02 

♦ 1.0 C3H6--=‘t.O 

H2o 

+ 3.U 

C02 

.llOO-fOlu 

• 000 

17.6000 


Ca) 

iO 




AXIAL 

LilSTAMCC 

At.Ofir, JTH , CYLI'^DrK OR 

PLATE 

S CCM) 


‘ 



. ij'sn 

.10(1 

.15ii 

.2U0 

.250 

.300 

.350 

.400 

.450 

• 5U0 

BULK tLHP |K) 


0?i,h 

832.3 

640,2 

ft47.2 

853.8 

659. P 

885.6 

871.1 

876.3 

WftUL TLMP (K) 

1197. fl 

121b. rt 

1226.0 

1231.7 

1235.3 

1257.5 

1259.1 

1240.2 

1241.5 

1243.8 

MOLL FRACTTOM 











BULK SPtCitS 











N2 

.7A314 

.7fl29’l 

.79266 

,76276 

.78266 

.78257 

,78248 

.78240 

.78233 

.78226 

02 

.20697 

.20^99 

.20517 

.20445 

.20360 

.20320 

.20264 

.20211 

.20161 

.20114 

C3H« 

,00607 

.00789 

.00772 

,00758 

.00745 

.U0734 

.00723 

.00713 

.00703 

•00694 

H20 

•00104 

.OOlflO 

.00243 

.00298 

.00348 

.00394 

,00437 

.00478 

.00516 

.00553 

C02 

,00076 

.00135 

.00102 

,00224 

.00261 

.00296 

,00528 

.00358 

,00367 

,00414 

WALL SPECIES 











N2 

.77667 

.77687 

.77687 

,77687 

.776^7 

.77687 

.77687 

.77687 

.77687 

,77687 

02 

.16^31 

.16531 

.16S30 

,16530 

,16530 

.16530 

.16530 

.16530 

,16530 

,16530 

C3H»i 

,00000 

.00000 

.00000 

,00000 

,00000 

.00000 

.00000 

.00000 

.00000 

,00000 

H20 

.03304 

.03304 

.03304 

,03304 

.03304 

.03304 

.03304 

.03304 

.03304 

•03304 

C02 

,0?47R 

.02978 

.02478 

,0247ft 

,02478 

.tiR478 

,02478 

.02478 

,02418 

,02478 




AXIAL 

.750 

r>isTAi'jn: 

1.000 

Al 0|\*r 40W0I. IlH t 1 

1.500 1.750 

CYLINIJCH oh 

2.ono 

PL>^Tf 

2.500 

S (Ch) 
3.000 

3,500 

4,000 

4.500 

bULK Tf.np (K) 

6*98.4 

916.0 

9S0 .6 

965.6 

979.7 

1004.9 

1027.5 

1048.0 

1066,8 

1084.0 

WALL TErP (K) 

1247.3 

1249.4 

1231.1 

lPbO.7 

l?3l.5 

1253,2 

1254,0 

1254,7 

1255.3 

1255.8 

MOLL FRACTION 

BULK SPECIES 

N2 

.76196 

,78169 

.78125 

.7810H 

.76085 

.78050 

.78019 

.77990 

.77964 

.77940 

02 

.19913 

.19737 

•19441 

,1930‘» 

.19175 

.18945 

.18737 

.18547 

.16374 

.18215 

C3H8 

,00^55 

,00621 

.00564 

,00537 

.00512 

,00467 

,00427 

,00391 

.00357 

.00326 

H2o 

.00706 

,008**2 

.01069 

.01175 

.01274 

.01450 

,01610 

,01755 

.oieea 

.02011 

C02 

•00530 

.00631 

.00802 

.00861 

.00^55 

.01088 

.01207 

,01317 

.01416 

.01508 

WALL SPECIES 

N2 

.77647 

,77667 

.77687 

,77687 

.77687 

,77687 

.77687 

.77687 

.77687 

.77687 

02 

.16530 

.16530 

.16530 

,16530 

•16530 

.16530 

,16530 

.16530 

.16550 

.16530 

CiHfl 

•00000 

.00000 

.00000 

,00000 

•00000 

,00000 

.00000 

.00000 

.00000 

.00000 

H20 

.03304 

.03304 

.03304 

,05304 

•03304 

,03304 

,03304 

,03304 

.03304 

•03304 

C02 

.02474 

.02474 

.02478 

.02478 

.02478 

,02478 

,02478 

.02478 

.02478 

.02478 




AXIAL 

B.OOll 

DISTANCr 

5.500 

ALOnr, MdMULITH • cyi inofjh oh 
fa. non 6,500 7,000 

PLATE 

7,500 

S (CM) 
6.000 

8.500 

y.ooo 

9.5U0 

Bulk temp «K) 

1099, n 

1114,4 

1127.9 

1140.4 

1152.1 

116«f.9 

1173.0 

1182,5 

1191,3 

1199,7 

WALL TEMP (K) 

12i)&.a 

1250. fa 

I2b7,0 

1257.3 

1257,6 

1257.0 

1256.0 

1258, 2 

1258.4 

1259,2 

MOLL FRACTION 

BULK SPECIES 

N2 


.7709P 

,77679 

.77H61 

• 77‘>44 

,77029 

.77015 

.77601 

.77700 

.77777 

02 

• Ifl 0 fa7 

.17931 

.17005 

,17687 

.17578 

.17475 

.17379 

.17290 

.17206 

.17127 


,00290 

.00271 

.00247 

,00224 

•00203 

.00103 

.00164 

,00147 

.00131 

.00116 

H20 


.02229 

.02326 

.02416 

• 025U0 

.02579 

.02652 

.02721 

.02706 

.02646 

C02 

,01b93 

.01671 

.01744 

,01812 

.01875 

.01934 

.01969 

.02041 

.02009 

.02135 

WALL SPECIES 

N2 

.776»57 

,77607 

.77667 

.77667 

.77607 

,77687 

,77687 

,77687 

,77607 

,77687 

02 

.16530 

,16530 

.16530 

, 16530 

.16530 

,16530 

,16530 

,16530 

.16530 

.16530 

C3HB 

•unonu 

,00000 

.oonnu 

,00000 

.unouo 

,00000 

,00000 

.00000 

.00000 

,00000 

H20 

,U33n^ 

,03304 

.03304 

.03304 

.03304 

,03304 

,U3304 

,03304 

.03304 

,03304 

C02 

.U2H 

,0247P 

.02476 

.02476 

,02476 

,02478 

,02478 

.02478 

.024/6 

,02478 



AXIAL LISTAUCF. AL^fiC, MOKOLIIH , CYLI'^iULK OP HLaTL 


S (Cn) 


1 0 • ouu 

JULK TLMP (K) 1207. fc 

WALL TErtP (K| 12f»0,B 

MOLL FRACTION 

bulk species 


N2 

.77766 

02 

.170b3 

C3HB 

•ooial 

H20 

. 02*^03 

C02 

.02177 

WALL SPECIES 

N2 

.77607 

02 

.16b50 

C3H« 

.00000 

H20 

.03304 

C02 

.02**76 


••****«*•«•*««•****«***** 

ITERATIONS = 1 TLMP ERROR = .11^31-001 





SAMPLE PROBLEM NO. 2 

NONADIABATIC CYLINDRICAL REACTOR (CERAMIC) 


Description: 


Fuel: 


Coolant: 


Geometry: 


This problem is the same as sample problem No. 1 except for a change in the fuel/ai 
conditions and the inclusion of a convective flow of air over the outer surface of 
the cylindrical reactors 


Propane/air 

Equivalence ratio 
Mixture flowrate 
Preheat temperature 
Pressure 



0.8 

0.0278 g/s 
lOOOK 

1 X 105 Pa 


Air 

Bulk temperature 

Heat transfer coefficient 


Tc = lOOOK , 

h(, = 1.51 X 10 ^ cal/cm‘^-s-K 


Cell diameter 
Reactor length 
Void fraction 


Dfub 0*32 cm 
£ = 20 cm 
Vp = 0.87 



0.8 


SAMPLE PROBLEM BZ: COOUD MONOLITH TUBE, EQUIVALENCE RATIO • 


f) 2 

1 0 1 

n 







0.317S 

0.87 

0,3404 

0,0005 

0.6 





2 , 7fl 

innn , 

1. 







1 nno • 

1500. 








0 ,n 

.O') 

,10 

.20 

.40 

,60 


,60 

1.0 

1.5 

2.0 

2.5 

3 ,0 

3,5 

4.0 


4,5 

5.0 

fi.O 

^».5 

7,0 

6.0 

^.0 

10,0 


11,0 

12.0 

13 ,n 

14,0 

15,0 

16,0 

16,0 

20.0 




1700. 

1700. 

1700. 

1700. 

1700, 

1700. 


1700. 

1700, 

1 700. 

1700, 

1700. 

1700. 

1700. 

1700, 


1700, 

1700. 

1 700. 

1700, 

17U0, 

1700, 

17 0 0, 

1700, 


1 700. 

1700, 

1700 . 

1700, 

1700, 

1700, 

1700. 

1700, 


1700, 

1700, 

1 1 1 1 : 

11111 

11111 

11111 

11111 

1 1 1 

1 1 

11111 

1 1 1 1 1 

o.onoifii 

0,000151 

0.000151 

0,000151 

0,000151 

0,000151 

0,000151 

0,000151 

0,000151 

0,000151 

0,000151 

0 ,000151 

0.000151 

0,000151 

0.000151 

0,000151 

0,000151 

0,000151 

0.000151 

0,000151 

0,000151 

0.000151 

0,000151 

0,000151 

0,000151 

0,000151 

0.000151 

0,000151 

0,000151 

0.000151 

0,000151 

0 ,000151 

1000, 

0.0005 

,01145 







U2 

,7#^4 

.7401 

1.03475 






02 

• 203 

.0394 

1,10045 






CSHfl 

.0325 

3, E-051.2758 






H20 

1,0 E- 

10 ,126 

.844069 






CJ2 

1.0 t- 

10 ,0945 

1,27446 






2 1 
C3HB 02 


C3H8 

02 H20 

C02 

4.1 

E + 9 

40.61 

.25 1.04 




-.75 -3,' 

964.0 

3,0 



C3H8 


OP 

H20 COP 


1.1 

e+09 0. 

17.6 

1.0 




-5.0 4,0 

3,0 






SAMPLE PROBLEM B2: COOLED MONOLITH TUBE, EQUIVALENCE RATIO - 0.8 


♦ INTEGRAL PARAMETERS • 


IS I NUMtJER OP SPECIES I = 5 

riL I GRTO POINTS I = 30 

NIT { NUrtRER OF OVERALL ITERATIONS » = 2 

ILCSS ( heat loss OPTION » = 1 

IGEOM i problem geometry option ) = 0 

ICCEFF C HEAT TRANSFER INPUT OPTION ) = 1 

KH7 I PRINT OPTION ) = 0 


♦ NOnINTEGKAL PARAMETERS * 

REACTOR data *•* 


42 * 

or 


• RF.SERVOIR temperatures ♦♦♦ 

UPSTREAM = 1000.0 

noUNSlREAM =. 1500,0 


CHIn POINTS CHOSEN (SI 


.OUUO 

.5000-001 

• lOUO-fOOO 

.2000 

• 40un 

• 6Uon 

2,f^00 

5.000 

3.50U 

4 .uon 

4 .500 

5.000 

9.000 

10.00 

11.00 

12.00 

13.00 

14.00 


IfilTT Al 

1700, 

guessed wall 

1700, 

TEMPEKArURES (Tl» 
1700. 1700. 

1700. 

17U0. 

1700, 

1700. 

17U0, 

1700. 

1700. 

1 700, 

1700 

, 1700. 

170U, 

1700. 

1700. 

1700. 


monolith honeycomb catalytic combustor 


NijMREK OF tubes 
VOID fraction 
OIAMETER OF COMBUSTOR 
EFFECTIVE CELL OIAMETER 


inlet temperature 

PPESSURF 
MASS flow rate 
conductivity 
emissivitt 


S 1*00 
= *870 

= .3** (CM) 

s ,32 (CM) 


slOOO.OO (K) 
s 1,00 (ATM) 

= ,278^001 (G/S) 

= .500-003 (CAL/S-CM-K) 
= .60 


• 6000 
6.000 
15,00 


1700. 

1700, 

1700. 


1.000 

6.500 

i6.no 


1700. 

1700. 

1700. 


1.500 

7.000 

18.00 


1700, 

1700. 

1700. 


2.000 
6.000 
20. OU 


1700. 

1700. 



♦ COOLANT f’KOPCK IIF.S 


CnCLANT 


inlet TLMPERATirHL 

=ioon. 

on { 

K) 

PPESSURL 

= 

00 I 

ATM> 

MASS Fl(.)W hail 

= .000 


<G/S) 

APFROACii VELOCITY 

^ « 

no ( 

CM/S» 

MOLECULAR UiEIGhT 

“ • 

00 


FOR COOLANT tijhe: 

TlieE OIAMETER 

= 

00 1 

CM» 

Tl.HE ThICKNF SS 

= .0 

11 1 

CMI 

COAUUCTIVIY 

= .bOo 

-003 

» FCAL/S- 


42 * 


* THCRMOCHEMICAL OATa FOK coolant 


tff'pehatukl thermal conductivity 

.00 . n 0 u 

.00 .000 - 


HFAT capacity cPH = cl ♦ C2 • T ♦ C3 / T **2 

Cl = .non (CAi./o-hOLE - K) 

c? = .non (CAi/r,-MOir - k**2i 

C3 = *000 (CAL-K/G-MOLE) 


Viscosity parameters 

SIGMA = .non (A| 

r./K = ,000 (K) 


(CAL/G-MOLE “ K> 



SHr.ci;', Irni;>L ,UJL: FK ACT IONS . H-IST CULSSFC . aMii IUFFDSI.IN FacTOHS 


U r- 3 - O' ^ iT) 

< 3- -= /I » » 

u. *o s r» 7 

o 1-1 rvi X ru 

O ^ ^ 


3 *-• 3 »-» 

C S 3 C O 
= 1 X 030 
I I •»> I 
o 3 o X tn 
3 O' C (V 3 
^• f} in ^ tr 


c c ^ 0 % 

^0030 

o c o o c 

♦ ■♦■III 

3 O J*. X O 

J X cv. c o 

r- r\i fO »H 


T X ”V* 
.% rv C 

^ w u X u 
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• 7Hl KMOCHLMISTHY DATA • 


CUK>/F hlJ iJF- DATA IN FORM CP=RU tRC +«0/ I T *T > t CAL/ < MOLE »K » 


HF<CAL/MQL) 

m(cal/mdl) 

RD 

RC 

HD 

S(CAL/M0L/KI 


TU(k) 

NAME 

N2 

./vuuu^iino 

?H,nnn 

,613^1+001 

,16415-002 

.29727+005 

,64548+002 

300, 

1000. 

N2 

790UO+000 

.22i«n+00S 

.83049+001 

,23357-003 

-^72606+006 

.63774+002 

1000. 

2500, 

N2 

U2 

onouo 

52,000 
,25541 *005 

.67044+001 

.16795-002 

-,19919^005 

.689084002 

300 • 

1000. 

02 

oooou 

,2345(UU05 

.60725+001 

,50513-003 

-.21319+006 

.67976+002 

1000, 

2500. 

02 

m20 

b779fl+On5 

i«,ono 

,31716+UOb 

.68290+001 

,29575-002 

,27b43^005 

•69022+002 

30U • 

1000. 

H20 

577VM+005 

.30303+005 

,99698+001 

,12774-002 

-.14328+007 

.68463+002 

1000, 

2500. 

H20 

C02 

44,000 

.4l43fl^005 

,95778+001 

.36020-002 

-,16389^006 

.81997+002 

300. 

1000. 

C02 

^4054+005 

,36573+005 

,13972+002 

.38449-003 

-.13412+007 

.79867+002 

1000. 

2500. 

C02 

C5H3 

24820+UOS 

44.000 

•lU868fU06 

.32u9u*0U2 

•11928-001 

-,2o372+o07 

,159404003 

300. 

2500. 

C3H8 

24820+005 

,10888+006 

,32090+002 

,11928-001 

-, 203724007 

,15940+003 

2500. 

6000. 

C3h8 



TCTpL NUrlHt < OF KtACTItiNS 
MiMHEK OF OAS PHASE KEACTIONS 


z u o o 

o u o o 

O r<t O 

£ sO sC 

< N • • 

> -I o r» 

^ .H 

»- o 

U X 

< -■ 


1 . o o 

X o o 

UJ o o 

• • 

CL 

i: 

u 


X 

o 

^ so 

U (/> ^ tH 

< t o o 

u. « ♦ ♦ 

u o o 

CL I C C 

K U <H 

UJ -I ^ r 4 

o • • 

u r 

X 

0 . 


CNJ M 

o o 
u u 

o o 




II II 

• I 

• I 

X X 
X X 
O K> 

u u 


c o 


c>j rs< 
o o 

o o 

if> IT 





fM 


50 



OVF^ALL HFaT THArjSFER RtSUl TS ♦♦ 
coolant TrrMLRATURL IN =1UUO,00 (K) 

COOLANT TfMFEKATURt OUT =10UO,00 (K) 

TOTAL HFAT THANSFEHEO TO 

COOLANT = ,451+001 (CAL/S ) 


51 



-\jiIaL liISlAMCr ALONG MOHOLITH 


CYl INf.rr oil f»LATt 


S (CM) 





,')^U 

.mu 

• 200 

.4U0 

• 60(J 

.800 

i.uon 

1.500 

2.000 

2,500 

HULK 

T) MP 

(K ) 

1 nol.i 

1U06 

1U13 ,0 

1026,7 

1040.6 

1054 .6 

1066.6 

1102.6 

1136,4 

1169.2 

wall 

TLMP 

(K) 


?312,6 

2364 

2421.6 

24b6,.’i 

2473.1 

2462,6 

2490.4 

2493.4 

2494.5 


hOLt ff^ACTTOM 

HULK S^’FrlLS 


Tj2 

.764 53 

.76427 

.76416 

.76394 

.76373 

.76351 

,76330 

.76278 

,76226 

.76176 

02 

.20273 

,20235 

.20160 

.20012 

, 19864 

.19718 

.19573 

.19217 

•16668 

,16525 

C3H4 

,03244 

.03237 

.03222 

.03192 

.03163 

.03134 

,0310b 

.03035 

.02965 

.02897 

H20 

,00029 

,00058 

.00115 

.00230 

.00343 

.00455 

.00567 

.00641 

•011U9 

.01372 

CO? 

,(l 0 022 

.00043 

.00066 

.00172 

.00257 

.00341 

.00425 

.00630 

•00832 

.01029 

WALL SPLCIES 

N2 

.74031 

,74031 

.74031 

.74031 

.74U31 

,74031 

,74031 

,74031 

.74031 

.74031 

02 

,U39?5 

.03^25 

,03925 

.03925 

•03925 

.03925 

,03925 

.03925 

.03925 

.03925 

C3M(‘ 

•uoono 

.00000 

.00000 

.UOOOO 

•00400 

.00000 

.00000 

.00000 

.00000 

.00000 

H?0 

. I2b4b 

.12596 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

CU2 

.US447 

,09447 

,09447 

.09447 

.09447 

.09447 

.09447 

.09447 

.094*^7 

.09447 

LOUl.Ai'iT 

TLMP (K) 

1 0 U U , 0 •) 

1 U 0 0 . U 0 

11)00.00 

loun.uo 

inun.oo 

loon.oo 

lonn.uo 

1000.00 

1000.00 

lUUO.OO 

CUOI.ANT Tll"L 

rr.nP (I- 1 

.0 * 

.00 

.00 

• 00 

• 00 

.00 

.00 

.00 

• OU 

• 00 

ilLnl Lf SS TM 

COnUA-n t( (.1 / : ) 

'll? 

.4'uo-nn? 

.?UhO-noi 

•4?93*n0l 

.4398-001 

. *4449-0 Hi 

,44 77-ODl 

•1125 

• 1126 

.1126 




AX I AL 
3. uon 

HIST aMCE 
3,500 

ALorjG nn'ioi iTH , ct 
4.000 4.500 

LINUrH OH 
5 • UOO 

plate 

6.000 

S (CM 
6,500 

7.000 

8.000 

9.000 

HULK TE flP (K) 

1201, 2 

1232,5 

I2b3.1 

1293.0 

1322.4 

1378,6 

1406,3 

1433,4 

1465.6 

1536,3 

WALL TFMP (K) 

n 

2495,3 

2495.6 

2496,4 

2497.1 

2496.2 

2497.5 

2496.3 

2499.8 

2500.3 

MULL FKACTIO^t 

HULK Sf»K.ClES 

N2 

./6127 

.7607ft 

,76031 

.75964 

.75939 

.75851 

.75808 

,75765 

.75683 

.75604 

02 


.17861 

,17538 

.17221 

.16911 

.16314 

.16019 

.15730 

,151/2 

.14630 

C3HA 

,02831 

,02766 

,02701 

,02639 

,02577 

.02459 

.02400 

.02343 

.02232 

.02124 

H20 

.01630 

•01683 

.02131 

.02375 

•02613 

•03072 

.03299 

.03521 

• 03950 

•0h367 

C02 

.01223 

,01412 

.01598 

.01781 

.01960 

•02304 

.02474 

,02641 

.02982 

•03275 

WALL SPrclLS 

m2 

,74051 

,74031 

.74031 

.74031 

.74031 

.74031 

.74031 

.74031 

.74031 

.74051 

02 

,03925 

.03925 

.03925 

.03925 

,03925 

,03925 

,03925 

,03925 

.03925 

.03925 

CSHf* 

•ooono 

•00000 

.00000 

•OOQOU 

•OOUOO 

•00000 

•00000 

.00000 

• 00000 

•00000 

H20 

.12597 

.12597 

,12597 

.12597 

.12597 

.12597 

.12597 

,12597 

.12597 

,12597 

C02 

.09447 

.09447 

.09447 

.09447 

•09447 

.09447 

.09447 

.09447 

.09447 

,09447 

coolamt 

TLr"H (K) 

inon, O'l 

1000,00 

1000,00 

lOOU.OO 

1000. 00 

1000.00 

1 (I n U • 0 0 

lOOO.OO 

lOUO.OU 

1000,00 

COOLANT TliOt 
TLMH (K) 

,00 

.00 

,00 

• 00 

.no 

.00 

.00 

.00 

.00 

.00 

HEAT LOSS 10 

COOLANT (CAL/S) 

,1129 

.1129 

.1129 

.1130 

.1130 

• 2262 

• 1131 

• 1131 

.2265 

• 2266 



rtXiAU JJISTANCr AL0^6 MorMOLlTH . CYLINDCH OR HLaTL S C CM » 



lu.ouu 

11.000 

12.000 

13.000 

14.000 

15.000 

16.000 

18.000 

20.000 

HULK FFMP <K) 

1585. €. 

1634,0 

1481.7 

1729.0 

1776.3 

1823.9 

1872.0 

1963.5 

2055.1 

wall TfHP (KJ 

2500./ 

2501.2 

2501.6 

2502.0 

2502.6 

2503.3 

2504.2 

2505.0 

2479.3 

molL fraction 

3ULK SPFCICS 
NP 

.75526 

.75450 

.75374 

.75299 

.75225 

.75149 

.75073 

.74928 

.74779 

UP 

•iHlUl 

.13581 

.13068 

.12559 

.12049 

.11536 

.11017 

.10031 

.09015 

C3MP 

.02019 

.01916 

•01615 

.01713 

•01612 

.01510 

.01408 

.01212 

.01010 

H20 

.047/3 

.05173 

.05567 

.05959 

•U635I 

•06745 

.07144 

.07903 

•06664 

COP 

•03580 

.03680 

•0H175 

•04469 

•04763 

•q5059 

.U5S5B 

•05927 

.06513 

*^ALL SPrCIES 
N2 

.74031 

.74031 

.74031 

.74031 

•74031 

.74031 

.74031 

.74031 

.74031 

OP 

.03925 

.03925 

•03925 

.03925 

.03925 

.03925 

.03925 

•03925 

.03925 

C3mP 

.uooou 

•oonoo 

•00000 

.uuooo 

•OUOOO 

•00000 

.00000 

.00000 

.00000 

MPO 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

COP 

•09446 

.09448 

.09448 

.09448 

.09448 

.09446 

.09448 

.09448 

.09448 

COol-ANl 

TtMP CK) 

luoo.un 

1 non. 00 

1 tiUO.OO 

loon . 00 

moo. 00 

1000.00 

looo.uo 

1000.00 

1000.00 

COf'LANT Tunr 

1 L 1 K I 

.n 1 

.on 

.00 

• 00 

• 00 

.00 

.00 

• 00 

.00 

Mu*. 1 1 r*;s TM 

Ct)Ml.^iJT (CAI./Sl 


.2>r,7 

.2267 

.2268 

,2269 

.2270 

.2271 

.4545 

.4467 



SAMPLE PROBLEM NO. 3 

NONADIABATIC FLAT PLATE REACTOR (STAINLESS STEEL) 


Description: 


Fuel : 


Coolant: 


Geometry: 


This problem represents a semi -infinite stream of a fuel /air mixture passing over one 
side of a catalytic flat plate. A flow of coolant (air) flows over the other side. 
The mass flowrate is given in terms of the unit cross-sectional flow area. 


Propane/air 

Equivalence ratio 
Mixture flowrate 
Preheat temperature 
Pressure 


Air 

Bulk temperature 

Heat transfer coefficient 


Plate thickness 

Initial noncatalytic length 


^ = 0.8 2 
m = 0.3574 g/cm'^-s 
Tin = lOOOK 
P = 1 X 1Q5 Pa 


Tc = lOOOK 

h(, = 1.51 X 10 ^ cal/cm^-s-K 


t = 0.3175 cm 
Hq = 30.5 cm 



SAMPLE PROBLEM B2-FLAT PLATE EQUIVALENCE RATIO • 0.8 


5 16 15 

111 

0 






1.0 

0.87 

1.0 

0.055 

0.6 





1000. 

1. 






1000. 

1000. 







0.0 

10. 

20. 

30.5 

30.51 

30,52 

30.53 

30.54 

30.55 

30.75 

31. 

31.5 

32. 

32.5 

33. 

34. 

1600. 

1600. 

1600. 

1600. 

1600. 

16Q0. 

1600. 

1600. 

1800. 

1600. 

1600. 

1800. 

160U. 

1600. 

1600. 

1600. 

0 0 0 1 

1 1 1 1 1 

11111 

11111 

11111 

1 1 1 1 1 

Ill 


0.000151 

0.000151 

0.000151 

O.OOOlSl 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

0.000151 

.3175 








N2 

.764 

.7401 

1.03475 





02 

.203 

.0394 

1.10045 





C3H8 

.0325 

5. r- 

051.2758 





H20 

1.0 E- 

10 .126 

.844069 





C02 

1.0 E- 

10 .0945 

1.27446 





2 1 








C3h« 02 


02 

H20 C02 


1.7 E*13 

21.0 

1.0 1.0 




-4.0 4.0 

3.0 



C3H8 


02 

H20 C02 


1.1 0. 

17.0 

l.O 




-5.0 4.0 

3.0 





SAMPLE PROBLEM B2-FLAT PLATE EQUIVALENCE RATIO * 0.8 


• integral parameters ♦ 


IS « NUMHER OF species ) = 5 
NL I GRID POINTS ) = 16 

NIT I number of overall iterations I = 15 
ILOSS « meat loss OPTION ) = 1 
IGEOM ( problem geometry option ) s 1 
ICOEFF t HEAT TRANSFER INPUT OPTION ) - 1 
KR7 I PRINT OPTION ) = 0 


• NONIKTEGRAL PARAMETERS ♦ 


*<♦ REACTOR data •** 


FCR CYLINDER COOLED EXTERNALLY OR FLAT PLATE PROBLEM 

THICKNESS OF CYLINDER OR 
THICKNESS OF PLATE s .317 CCM) 


cn 


INLET temperature 

PRESSURE 

MASS FLOW RATE 

conductivity 

ehissivity 


slOOO.OO (K) 
r 1.00 (ATMI 
s . 357+00 CG/S» 

= .550-01 ICAL/S-CM-K) 

s .60 


reservoir temperatures 

UPSTREAM s 1000,0 

OCWNSTREAM = 1000.0 


GRID POINTS CHOSEN (S» 

.0000 lO.nn 20«00 SO.SO 30.5l io.52 

31.00 31,50 32.00 32.50 33,00 34.00 


INITIAL GUESSED HALL TEMPERATURES (TL) 

1800. 1600. 1600. 1600. 1600. 1600. 
1800. 1600. 1600. 1600. 1800. 


30.53 30.54 30.55 30.75 


1800 


1800 


1600 


1600 



••• COOLANT PKoPLKTIES ••• 


COOLANT 

t^LET TEMPERATURE = 
PRESSURE. = 

MASS FLOW PATE = 
approach VELOCITY = 
MOLECULAR WEIGHT s 

FOR COOLANT TUBE I 
TUBE diameter s 

TL8E THICKNESS = 
CCNOUCTIVIY = 


.00 (K1 
•00 (ATMl 
.000 (G/Sl 
.00 (CM/SI 

• on 


•00 (CMI 
.317 (CMI 

.550-01 ICAL/S-CM-K) 


CENTER TO CENTER DISTANCE I TUBE RANK) r ,0U JCM) 
distance between monolith and tube bank « .000 (CM) 


(J1 

00 


* therhochemical data for coolant 


TEMPERATURE 

.00 

.00 


THERMAL conductivity 
• 000 
.000 


heat capacity CPH 3 Cl ♦ C2 • T ♦ C3 / T**2 (CAL/G-MOLE • K) 


Cl 3 *000 
C2 3 ,000 
C3 3 ,000 


(CAL/G-MOLE - K) 
(CAL/G-MOLE • K«*2) 
(CAL-K/G-MOLE) 


VISCOSITY parameters 

SIGMA 3 .000 (A) 

E/K s .000 (K) 



♦ spfCiES initial molt fhacttons.first rulsses.anu diffusion FACTOKS • 


SPECIE 

ALPF 

ALPL 

DIFFUSION FACTOH 

N2 

, 764 + 00 

,740+00 

1.0347 

02 

,203+00 

,394-01 

1 , 1004 

C3H8 

,323-01 

,500-04 

l.275fl 

H20 

,100-09 

.126+00 

.6441 

C02 

.100-09 

.945-01 

1,2745 



♦ THEPWOCHEhlSTRY DATA ♦ 


CUHVr FIT OF DATA IN FORM CP=R0+RC*T>RO/ ( T • T » I C AL/ I MOLE*K » 


HF(CAL/MOL) 

H(CAL/M0L> 

RB 

RC 

Rn 

S(CAL/MOL/K> 

tuik » 

NAME 

N2 

.7V000+00 

28.000 

.23993405 

.61391401 

.16415-02 

.29727405 

.64546402 

300. 

1000. 

N? 

•79000^00 

•22188405 

.83049401 

.23357-03 

-.72806406 

.63774402 

1000. 

2500, 

N2 

02 

•00000 

32.000 

.25541405 

.67044401 

.16795-02 

-.19419405 

.68908402 

300. 

1000. 

02 

•ooono 

•23450405 

.80725401 

.50513-03 

-.21319406 

.67976402 

1000. 

2500. 

02 

H20 

..57796>05 

18.000 

•51716405 

•68290401 

.29575-02 

.27843405 

.69022402 

300. 

1000. 

H20 

►,57798f05 

.30303405 

.99698401 

.12774-02 

-.14328407 

.68463402 

1000. 

25U0. 

H20 

C02 

>•94054405 

44.000 

•41438405 

.95778401 

.36020-02 

-.16389406 

.81997402 

300. 

1000. 

C02 

>•94054405 

•36573405 

.13972402 

.38449-03 

-.13412407 

,79867402 

1000. 

2500. 

C02 

C3H8 

-•24820405 

44.000 

.12T85406 

.32090402 

.11928-01 

-.20372407 

.15567403 

300. 

2500. 

C3H8 

>•24820405 

.12785406 

.32090402 

.11928-01 

-.20372407 

,15567403 

2500. 

6000. 

C3H8 



• kinetic KLACTInN DATA ♦ 


TOTAL NUPIHEH OF REACTIONS 2 

NUfIRER OF GAS PHASE REACTIONS 1 


M 



REACTION 



PRF rXP FACTOR 
(MOLF-CM-S) 

TEMP EXP 

ACTIVATION 

(KCAL/MOLF.) 

1 

2 

GAS 

SURF 

S.O U2 
b.O 02 

♦ 1.0 C3Hfi--=H.O 

♦ l.O C3He--=4.0 

H20 

H20 

♦ 3.0 C02 
4^ 3.U C02 

.17004>m 

.IIOO+IO 

.000 
• QUO 

21.0000 

17.0000 


or 



ovehall heat transfer results •• 

COOLANT TEMPERATURE IN = .00 

COOLANT TEMPERATURE OUT = .00 IK> 

total heat transfereo to 

COOLANT = .5?5^0l ICAL/S) 


62 






AXIAL 

DISTANCE 

ALONG 

monolith , 

CYLINDEH 

OR PLATE 

S <ch» 







10.000 

20.000 

30. SCO 

30.510 

30,520 

30.530 

30.5**0 

30.550 

30.75U 

31,000 

BULK 

TEMP 

(K| 

1000.0 

1000*0 

1000.0 

1000.0 

1000.0 

lonn.o 

1000.0 

loon.o 

1000.0 

1000.0 

WALL 

TEMP 

(K) 

988.7 

988.7 

1006.6 

1763,5 

1762,2 

1758,2 

1751,6 

1742.3 

1646.5 

1516.4 


MOLE FRACTION 

BULK SPECIES 


N2 

,76400 

.76400 

.76400 

.76400 

.76400 

.76400 

.76400 

,76400 

.76400 

,76400 

02 

•20300 

.20300 

.20300 

.20300 

.20300 

.20300 

.20300 

,20300 

.20300 

.20300 

CSh8 

.03250 

.03250 

.03250 

.03250 

.03250 

.03250 

.03250 

.03250 

.03250 

.03250 

H20 

.00000 

.00000 

.00000 

.00000 

.00000 

,00000 

.00000 

•00000 

.00000 

,00000 

C02 

.00000 

.00000 

.00000 

.00000 

»00000 

.00000 

.00000 

.00000 

.00000 

,00000 

wall species 











N2 

•76436 

.76438 

.76438 

.74031 

.74031 

.74031 

.74031 

.74031 

.74031 

.74031 

02 

• 20310 

•20310 

•20310 

.03924 

.03924 

.03924 

.03924 

.03924 

.03924 

.03924 

C3H8 

.03252 

.03252 

.03252 

.00000 

.onooo 

.00000 

.00000 

.00000 

•00000 

.00000 

H20 

.00000 

•00000 

•00000 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

C02 

• 00000 

•00000 

.00000 

.09448 

.09448 

.09448 

.09446 

.09448 

.09448 

.09448 

coolant 

TEMP (Kl 

• 00 

.00 

.00 

• 00 

.00 

.00 

.00 

.00 

.00 

.00 

heat loss To 
coolant (CAL/S) 

1.493 

1.493 

1,596 

.2663-02 

•26bl-u2 

,2655-02 

.2645-02 

.2631-02 

,4972-01 

.5/24-01 



AXIAL OTSTANCr ALONG hONOLlTH 


CYLlNr)rK OK PLATL 


S ICM> 



31.500 

32.000 

32.500 

33.000 

34.000 

bulk temp ik» 

1000.0 

1000.0 

1000.0 

1000.0 

1000.0 

WALL TEMP <K) 

1366.3 

1254.1 

1188.6 

1150.6 

1120.0 

MOLE FRACTIOA 

Bulk species 

N2 

.76H00 

.76400 

.76400 

.76400 

.76400 

02 

•20300 

•20300 

•20300 

•20300 

•20300 

C5hB 

•03250 

•03250 

.03250 

•03250 

.03250 

H20 

•00000 

• 00000 

•OQOOO 

•00000 

•00000 

C02 

•00000 

• 00000 

•00000 

•00000 

.00000 

WALL SPECIES 

N2 

•7H031 

.74031 

.74031 

.74031 

•74031 

02 

.03^24 

.03924 

.03924 

.03924 

.0^925 

C3H8 

•OOOQO 

•00000 

.00000 

.00000 

.00000 

H20 

.12597 

.12597 

.12597 

.12597 

.12597 

C02 

.09448 

.09448 

.09448 

.09448 

.09448 

coolant 

TEMP CK) 

• 00 

• 00 

.00 

.00 

.00 

HEAT LOSS TO 

coolant (CAL/SI 

.1032 

.9468-01 

.8974-01 

.8687-01 

.1691 


iterations 


15 


TEMP error 


iflosa-o? 



SAMPLE PROBLEM NO. 4 
COANNULAR TUBES 


Description: 


Fuel : 


Coolant: 


Geometry: 


In this problem, a coolant passes axially through the ID of the inner tube and a fuel/ai 
mixture passes axially and in the same direction in the annulus formed by the coannul ar 
tubes. The outer surface of the larger tube is adiabatic and the inner surface is cata- 
lytic. Heat is transferred from the larger tube to the smaller tube by radiation and 
convection. 


Propane/air 

Equivalence ratio 
Mixture flowrate 
Preheat temperature 
Pressure 


Air 

Velocity 

Flowrate 

Initial temperature 
Pressure 

Molecular weight 


Coolant tube conductivity 
Coolant tube thickness 
Coolant tube outside dia. 

Reactor tube inside dia. 
Length 

Catalyst begins 
Catalyst ends 


(^ = 0.8 2 
m = 1.045 gm/cm'^-s 
Tin = 866K 

P = 1 X 10= Pa 


Vjj, = 3000 cm/s 
iTir = 2.046 g/s 
T,* = 294K 
P^ = 1 X 105 Pa 
M = 28.84 g/mole 


kj = 0.055 cal/s-cm-K 
Dp = 0.05 cm 
Dj = 0.925 cm 

DTub = 1*70 cm 
il = 61 cm 
= 5 cm 
8-2 “ 50 cm 



SAMPLE PROBLEM D-COANNUIAR TUBE EQUIVALENCE RATIO 

■ 0.8 




5 2B 5 

0 2 0 

0 






X.5875 

0.87 

1.702 

0.055 

0.6 




l.0«»5 

866. 

1.0 






666.0 

-1.0 







0.0 

2.5 

5.0 

5.1 

3.2 

5.3 

5.4 

5.5 

5.75 

6. 

8.5 

T. 

a. 

10. 

12. 

15. 

20. 

25. 

30. 

55. 

40. 

45. 

50. 

50,6 

51. 

55. 

60. 

61. 





X700, 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700. 

1700, 

1700. 

1700. 

1700. 

0 0 0 1 

11111 

11111 

11111 

11110 

0 0 0 

0 


3000. 

2.0H6 

294. 

1.0 

28.84 




0.055 

0.05 

0.925 

1 . 

1. 




373. 

573. 

0.000073 

0.000155 





3.711 

78.6 

6.9 

0.00092 

•18504 

. 



N2 

.76% 

• 7401 

1.03475 





02 

• 203 

.0394 

1.10045 





C3H6 

.0325 

5. E-051.2798 





H20 

1.0 c- 

10 .126 

,844069 





C02 

1*0 E- 

10 .0945 

1,27446 





2 1 








C3h6 02 


C3M8 

02 H20 

C02 

4.1 

C49 

40.61 

.25 1.0% 




•.75 -B.* 

96 4. 

3.0 


C3H6 


02 

H20 C02 


1.1 

E409 0. 

17.6 

1.0 




•5.0 4.0 

3.0 





SAMPLE PROBLEM D-COANNUALR TUBE EQUIVALENCE RATIO - 0.8 


* INTEGRAL PARAMETERS • 


IS ( NUMbEH OF SPECIES ) = 3 
NL i GRID POINTS ) a 28 
NIT ( number of overall ITERATIONS ) = 5 
ILOSS ( HEAT LOSS OPTION ) = 0 
I6EOM ( PROHLEM GEOMETRY OPTION I a 2 
ICOEFF i HEAT TRANSFER INPUT OPTION ) s 0 
KR7 i PRINT OPTION ) * 0 


• N0NI^TE6RAL PARAMETERS • 

*•« REACTOR DATA •*« 


MONOLITH HONEYCOMB CATALYTIC COMBUSTOR 
NUMBER OF TUBES a 1.00 

VOID FRACTION ■ .870 

DIAMETER OF COMBUSTOR a 1.70 (CM) 

EFFECTIVE CELL DIAMETER a 1.3S (CM) 


INLET temperature a 866.00 (K) 

CT» PRESSURE B 1.00 (ATM) 

MASS FLOW RATE a ,105^0l (6/S) 

CONDUCTIVITY a .SSO-OX (CAL/S«CM-K) 

EFISSIVITY a .80 


RESERVOIR TEMPERATURES 

UPSTREAM a 666.00 

DOWNSTREAM a 1.0000 


GRID POINTS CHOSEN (S) 

.0000 2.300 S.OOO 3.100 3.200 3, SOU 

6.300 7.000 8.000 10.00 12.00 13.00 

HO.OO 85.00 50.00 50.60 51.00 35.00 


INITIAL GUESSED WALL TEMPERATURES (TL) 

1700. 170U. 1700. 1700. 1700. 1700. 

1700. 1700. 1700. 1700. 1700. 1700. 

1700, 1700. 1700. 1700. 1700. 1700. 


5,800 3.500 3,730 6,000 

20.00 23*00 30.00 35.00 

60.00 61.00 


1700. 

1700. 

1700 


1700, 

1700 


1700, 

1700 


1700* 

1700 



COOLANT PROPEHTIES 


COOLANT 

IKLET TEMPERATURE » 294.00 (Kl 
PRESSURE 8 1.00 (aTH) 

MASS FLOW RATE s .2U5t01 (G/SI 
APPROACH velocity kSOOO.OO (CR/SI 
MOLECULAR WEIGHT 8 28.64 


FOR coolant tube* 
TLBE UIAMETER 
TUBE THICKNESS 
CCNDUCTIVIT 


•92 (CM) 

•050 lent 

.590-01 JCAL/S-CH-K) 


• THERMOCHEMICAL DATA FOR COOLANT 


temperature 

373.00 

973.00 


THERMAL CONDUCTIVITY 
.730-04 
.155-03 


HEAT CAPACITY CPH m Cl ♦ C2 ♦ T ♦ C3 / T**2 ICAL/G-MOLE - K) 


Cl a .690^01 ICAL/G-MOLE - K) 

C2 ■ .920-03 ICAL/6»H0LE - K**2> 

C3 a ••165405 ICAL-K/G-MOLE ) 

VISCOSITY PARAMETERS 

SIGMA 8 3*711 CA) 

E/K 8 76.600 IK) 



• species initial note fractions«first 6uesses«ANo diffusion factors • 


specie 

ALPF 

ALPe 

DIFFUSION FACTOR 

N2 

•76H40Q 

•THO^OO 

1^03%7 

02 

•203^00 

•39%. 01 

1^100% 

C5H8 

•325*01 

•500*0% 

1^27S6 

H20 

•100*09 

•126^00 

• 0%%! 

C02 

•100-09 

•9H5-01 

l,27%5 



• thermochemistry data • 


CURVE fit of data IN FORM CPsRB4RC*T^HD/ ( T#T I < CAL/ I nOLE«K I 


HF(CAL/fOL» 

H(CAL/MOLI 

RB 

RC 

KU 

SICAL/MOL/K) 

TUIK » 

NAME 

N2 

•79000400 

28,000 

•23993405 

•61391401 

,18915-02 

.2972^405 

•89598402 

300. 

1000. 

N2 

,79000400 

•22186405 

•83U99^U1 

.23357-03 

-• 72808408 

,83779402 

1000. 

2500, 

N2 

02 

,00000 

32,000 

•25591405 

•87099^01 

.18795-02 

-•19919405 

•68908402 

300. 

1000. 

02 

,00000 

.23950405 

•80725401 

.50513-03 

-,21319408 

.67978402 

1000. 

2500. 

02 

H20 

.57790405 

xa»ooo 

•31718405 

•88290401 

•29575-02 

•27893405 

.69022402 

300, 

1000, 

H20 

.57796405 

•303U3405 

•99898401 

•12779-02 

-•19328407 

•88983402 

1000, 

2500. 

H20 

C02 

,99059405 

99,000 

•91938405 

.95778401 

.36020-02 

-.18389406 

.81997402 

300. 

1000, 

C02 

,99059405 

•38573405 

.13972402 

•38999-03 

-.13912407 

.79887402 

1000, 

2500. 

C02 

C3HS 

•29820405 

99.000 

•10888408 

•32U904U2 

.11928-01 

-.20372407 

•15990403 

300, 

2500. 

C3H8 

.29820405 

•10888408 

•32U904U2 

•11928-01 

-.20372407 

.15990403 

2900, 

8000. 

C3M8 



• kingtic reaction data • 


TOTAL NUhBER OF REACTIONS 2 

NUHBER OF GAS PHASE REACTIONS 1 

activation 
(KCAL/R0LE> 


REACTION 


PRE EXP Factor 

(ROLE-CR-S) 


TEHP EXP 


1 

GAS 

5.0 02 

♦ 1.0 C3H6--=H.O 

H20 

♦ 3.0 

C02 

•41004010 

• 000 

40.4100 

2 

SURF 

5.0 OP 

♦ 1.0 C3H8 — =H.O 

M2o 

♦ 3.0 

C02 

•110Q401U 

• 000 

17.4000 





OVtKALt HEAT TRANSFER RESULTS *• 
COOLANT TEnPCKATURC IN s 294.00 (K) 
COOLANT TEMPERATURE OUT s SOO.S8 CK) 
total heat TRANSFEREU TO 


72 



AXIAL DISTANCE 


ALONG MONOLITH « CYLlNOEK OR PLATE S ICMI 



2.500 

5.000 

5.100 

5.200 

5.300 

5.400 

5.500 

5.750 

4.000 

6.500 

bulk temp (KI 

6HS.5 

844.4 

848.3 

652.6 

857.0 

841.5 

866.0 

877.1 

6B7.5 

904.4 

WALL TEMP (K) 

dio.e 

936,9 

1511.2 

1555.3 

1583.6 

1596.9 

1594.9 

1579.1 

1557.5 

1525.4 

MOLE FRACTIOK 

BULK SPECIES 

N2 

•7S43H 

.74434 

•76434 

.74421 

.76409 

.74397 

.76385 

.74356 

.74327 

,T42T3 

02 

•20279 

•20279 

.20279 

• 20195 

•20111 

•20029 

.19947 

.19749 

.19555 

,19U7 

C3H8 

•03244 

•03245 

.03245 

•03229 

.03212 

•03194 

.03180 

•03140 

.03102 

.09029 

H20 

•00024 

•00024 

.00024 

.00089 

.00153 

•00214 

.00279 

•00432 

.00580 

.00043 

C02 

•OOOIB 

•00018 

.00018 

.00047 

.00115 

,00142 

.00209 

.00324 

.00435 

.00400 

WALL SPECIES 

N2 

•74454 

.74434 

.74434 

.74031 

.74031 

.74031 

.74031 

.74031 

.74031 

.74031 

02 

•20279 

.20279 

•20279 

• 03925 

.03925 

•03925 

•03925 

.03925 

•03925 

.03925 

cshs 

•03244 

•03245 

.03245 

.00000 

.00000 

.00000 

• 00000 

.00000 

.00000 

.00000 

H20 

•00024 

.00024 

.00024 

•12597 

.12597 

.12597 

.12597 

.12597 

.12597 

.12597 

C02 

•oooia 

•00018 

• 00016 

•09448 

.09446 

.09448 

.09448 

.09448 

.09448 

.09448 

coolant 

TEMP IK) 

5ll^07 

328.73 

329.86 

331.06 

332.33 

333.64 

334.97 

338.54 

342.07 

349.23 

coolant TUBE 

TEMP IK) 

653,51 

496.66 

436.91 

985.20 

1U21.43 

1046.89 

1063.20 

1119.01 

1120.86 

1144.12 

HEAT LOSS TO 

coolant ICAl/S) 

0«H2H 

8.765 

• 5600 

• 6023 

• 6331 

• 6539 

.6661 

1.781 

1.766 

3.583 



AXIAL distance ALONG MONOLITH « CYL1NDE« OK PLATE 


S ICMI 



7.000 

0^000 

10.000 

12.000 

15,000 

20.000 

25,000 

30.000 

35.00U 

40,000 

BULK TEnp CKI 

^23*H 

991.9 

994,0 

1026.1 

1U59.2 

1102.4 

1120.5 

1143.1 

1149.3 

1148.9 

WALL TEMP (K) 

i**se*H 

1459.9 

1411,4 

1375.5 

1331,1 

1276.3 

1230.6 

1190,5 

1154.4 

1121.2 

MOLE FRACTION 

BULK SPECIES 

N2 

.76222 

.76127 

.75966 

.75825 

,75648 

.75416 

.75225 

.75066 


.74818 

02 

•18096 

.10194 

•17095 

•16136 

•14930 

,13349 

•12054 

•10970 

.10056 

.09281 

C9H6 


.02892 

•02614 

•02429 

•02164 

•01670 

• 01613 

•01998 

• 0W17 

.01069 

H20 

•01139 

•01627 

•02472 

• 09209 

•04136 

•05351 

•06347 

•07180 

,07.05 

.08479 

C02 

•00050 

• 01220 

•01654 

• 02407 

•09102 

.04014 

•04760 

•05985 

,05912 

.06959 

wall species 

N2 

.74031 

•74091 

.74091 

.74091 

•74031 

.74031 

,74031 

.74091 

•74091 

.74091 

02 

•09925 

•09925 

•09925 

•03925 

•09925 

.03925 

•09925 

.09925 

•09925 

.09925 

C5H8 

•00000 

•00000 

•00000 

.00000 

,00000 

.00000 

•00000 

•00000 

•00000 

.00000 

H20 

•12597 

.12597 

•12597 

.12597 

.12597 

.12597 

.12597 

•12597 

•12597 

.12597 

C02 

•09448 

•09448 

•09448 

.09448 

•09448 

•09446 

•09446 

.09448 

•09448 

.09448 

coolant 

temp ik» 

356^I0 

969.43 

994.18 

416,66 

446,92 

454,09 

460.96 

467.55 

479,08 

479.95 

coolant tube 

TEMP 

1129^H4 

1129.84 

1117.06 

1101.69 

1084,62 

1212.99 

1160.56 

1150.90 

1129.19 

1096.77 

HEAT LOSS TO 

COOLANT (CAL/Sl 

3,4H5 

6.700 

12.47 

11.99 

15.41 

3,670 

9.524 

3.386 

9.253 

3.129 



AXIAL OISTANCC ALONG MONOLITH • CYLINDER OR PLATE 


S CCM) 



ifS.OOO 

50.000 

50.600 

51.000 

59.000 

60.000 

61.000 

•ULK TCnp (Kl 

114S.5 

1133.6 

1130.9 

1130.0 

1103. U 

1069.1 

1062.6 

WALL Ttnp (K» 

1090.5 

1061.1 

1032.1 

998.1 

686.1 

852.1 

858.1 

MOLE FRACTION 

BULK SPECIES 

N2 

.71721 

.71610 

.71627 

.71627 

.71622 

.71617 

.71617 

02 

•0SS21 

•06067 

.07961 

.07979 

.07919 

.07916 

.07913 

C3Hd 

•0093S 

• 00622 

.00605 

.00805 

,00796 

,00792 

•00791 

H20 

.0B961 

.09112 

.09178 

.09180 

•09506 

.09528 

.09531 

C02 

•06736 

•07059 

•07109 

• 07110 

• oaso 

•07116 

•07118 

WALL SPECIES 

N2 

•71031 

•71031 

•71031 

.71627 

.71622 

.71617 

•71617 

02 

•03925 

•05929 

• 03925 

.07979 

.07915 

.07916 

•07913 

C3HS 

•00000 

• 00000 

.00000 

• 00005 

,00796 

•00792 

.00791 

M20 

.12997 

.12597 

.12997 

.09180 

.09506 

,09528 

•09531 

C02 

•09H1S 

•09116 

.09118 

.07110 

.07130 

,07116 

.07118 

coolant 

TEMP (Kl 

165.77 

191.33 

192.16 

192.35 

195*62 

199.61 

900.56 

coolant tube 

TEHP (Kl 

1071.62 

1046.91 

1009.07 

959.61 

921.58 

892.92 

676,93 

MEAT LOSS TO 

coolant (CAL/SI 

2.996 

2.668 

• 1306 

.9716*01 

1.790 

2.062 

• 3971 



SAMPLE PROBLEM NO. 5 

CYLINDRICAL REACTOR (CERAMIC), ASSIGNED BED TEMPERATURE 


Description: This problem represents a fuel air mixture flow in a catalytic ceramic tube with no 

heat loss to the environment, and the catalyst bed temperature is assigned at 1200K. 
This problem is also the same as the catalytic combustion in a honeycomb reactor. 


Fuel : 


Geometry: 

•VI 

cn 


Propane/air 

Equivalence ratio 
Mixture flowrate 
Preheat temperature 
Pressure 


Cell diameter 
Reactor length 
Void fraction 


(}. = 0.2 

rh = 0.0394 g/s 
Tin = book 
P = 3 X 105 Pa 


DTub “ 0.32 cm 
I = 2.0 cm 
Vp = 0.87 



f’Lr l>«(h l.Lii - ASSlGr-LO IlMf'. UMIun 


in 

«o 

o 


• C O i-t 

(\J 3 o 
oj ry 
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IfjTf.TKrtL MArtAhtTtKS 


r) tT rj o 3 o 


li II II II II II II 


-H Z 0. 

o o 

<I — 


J U. O C «/; 

^ y, D <. O 

n V ’J' 1 ^ 

X Z C -J lJ C. 




i- 13 ^ u ^ 

1 / _ — t 2 L. <x 

— Z z -N i-i P-« X 


^ 4 

tH fO 
O O 

o o o o o 

O C I I c 

• • ^ o • 

O fO > O 

o -o r 


r •»■>*- 

xJ uJ 3 
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SPVLil.S luniAL ^0L^ J WALT GUE bSJ S • ANt t niFMjSION FACTORS 


X ^ in 

X C ./I T 3> 

.-o 3 r- 3 - 

o 1-1 rsj X rvj 


o 3 ,H 

o c o o o 

0 3 0 3 0 

^ 4> 4- I t t 

X f- n o -• ® 

_J r- vO o V) ^ 

c r- ry »o ca 


c cr rv- T> J' 

c o c o o 

0 3 0 0 0 

11 . 4 - > I I I 

X «0 ® -O 3 O 

_» ® 3 *r 3 o 

h- X «-l 


X 3 ry 
r, X 3 
^ J U X o 
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• THLRriOCHL.’lISTRY DATA ♦ 


CUHVr fit of UATA in form CP=Rb^RC*T4HO/(T*T) »CAL/(MOLE*K » 


HF iCAL/f'CL » 

HICAL/MOL » 

HB 

RC 

HO 

SICAL/MOL/KI 

tuiki 

NAME 

L2 

,7vuuo+ono 

28.000 

.239934005 

.613914001 

,16415-002 

.297274005 

.645484002 

300. 

1000. 

N2 

TSUUOfUOO 

.221884005 

,630494001 

.23357-003 

-,7-:B064U06 

.637744002 

1000. 

2500. 

n2 

02 

,oouoo 

32.000 

.255414005 

,670444001 

•16795-002 

-.194194005 

.069084002 

300. 

1000. 

02 

,oouoo 

.234504005 

.807254001 

.50513-003 

-.213194006 

.679764002 

1000. 

2500. 

02 

H20 

,b77‘?ft^U05 

18.000 

,317184005 

,682904001 

.29575-002 

,278434005 

.090224002 

300, 

1000. 

H20 

,577Sfl400b 

.303034005 

.996984001 

.12774-002 

-.143284007 

•684634002 

lOOO. 

2500. 

H20 

CO? 

,940544005 

44,000 

,414384005 

,957784001 

.36020-002 

-.163894006 

.819974002 

300. 

1000. 

C02 

,940544005 

.365734005 

.139724002 

.38449-003 

-.134124007 

.798674002 

1000. 

2500. 

C02 

C3Hft 

,24820+005 

44.000 

.108884006 

.320904002 

•11928-001 

-.203724007 

.159404003 

300. 

2500. 

C3H8 

,248204005 

.108884006 

,320904002 

•11928-001 

-.203724007 

• 159404003 

2500. 

6000. 

C3H8 



KINtLTIC KEACTION DATA 


Z LJ O O 

O -I O O 

i-t O O 

»- Z 

< N. • • 

> -j o r- 

i-i < » f-t 

t- u 
U 

< — 


a. o o 

X o o 

UJ o o 

• • 

0 . 

c 

UJ 



M <M 

o o 

U U 

CD CD 

• • 

fO fO 




i: o 
r) 

2 X 


o o 
:m cm 
X X 

0 o 

^ 3 - 
It II 
t I 

1 ■ 

X s 
X X 
fO K) 

u u 


o o 



CM CM 

o o 

o o 


-IX rfi m 

« r 

t- D 

K U. U. 

a CT X 

>- -c — 

I- a (/3 
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Axial 

• 05U 

lUSTANCE 
• 100 

ALONG monolith • CYLINDEH OK 
•150 .200 .250 

PLATE 

.300 

S ICM) 
.350 

.400 

• 450 

• 500 

liUL^ TEMP IK) 


823.2 

631 ,3 

838.4 

644.8 

850.7 

856.2 

861.4 

866.3 

870.9 

WALL TEMP IK) 

1200.0 

1200,0 

1200.0 

1200.0 

1200.0 

1200.0 

1200.0 

1200.0 

1200.0 

1200.0 

MOLL FKACTION 

hulk species 

N2 

.Tail** 

•7fl299 

.78287 

.78276 

.78266 

.78257 

.78249 

.78241 

.78233 

.78226 

02 

,20b98 

•20800 

•20518 

.20446 

.20382 

.20322 

.20266 

.20214 

.20164 

.20116 

C3H8 

•00807 

,00788 

•00772 

.00758 

.00746 

.00734 

.00723 

.00713 

• 007U3 

.00694 

H20 

•00104 

.00179 

•00242 

•00297 

.00347 

•00393 

.00435 

.00476 

• 00514 

.00550 

C02 

•00078 

•00134 

•00181 

.00223 

.00260 

.00294 

.00326 

.00357 

•00385 

.00413 

WALL SPECIES 

N2 

•77687 

•77687 

•77687 

.77687 

.77687 

.77687 

.77687 

.77687 

.77687 

.77687 


• 16!j31 

•16531 

•16531 

•18530 

•16530 

.16530 

•16530 

•16530 

.16530 

.16530 

CJHfl 

• u n 0 Q u 

,00000 

•00000 

.00000 

•00000 

•00000 

•00000 

.oooou 

.00000 

.00000 

H20 

•03^04 

•03304 

•03304 

•03304 

•03304 

.03304 

•03304 

.03304 

•033U4 

.03304 

C02 

, 02478 

.02476 

.02478 

•02478 

.02478 

.02478 

•02478 

.02478 

•02478 

.02478 



00 

OJ 





axial 

Instance 

ALONG 

monolith , CYLINDEK oh plate 




• 750 

1 .000 

1.500 

1.750 

HULf' 

n MP 

IK ) 

890. A 

907.2 

935.3 

948.4 

WALL 

TEMP 

IK) 

1200.0 

1200.0 

1200.0 

1200.0 


MOLt FHACnON 

HULK SPECirS 


N2 

.7819G 

.78170 

.78127 

.78106 

02 

•19919 

.19744 

.19453 

.19317 

C3H8 

•00&56 

.00622 

• 00566 

•00540 

H?0 

•007U2 

• 00636 

• 01059 

•01164 

C02 

.00527 

•00627 

.00795 

• 00873 

wall species 

N2 

.77687 

• 77687 

.77687 

.77687 

02 

•16530 

.16530 

•16530 

.16530 

C3H8 

•00000 

• OOOOO 

.00000 

• 00000 

H20 

•03304 

•03304 

.03304 

.03304 

C02 

.02478 

•02476 

.02478 

•02478 


S (Ch) 



SAMPLE PROBLEM NO. 6 

CATALYTIC REACTOR WITH UP- AND DOWNSTREAM HEAT EXCHANGER 


Description: 


Fuel : 


Coolant: 


Geometry: 


In this problem, two banks of heat exchanger tubes are located up- and downstream of 
a monolith catalytic reactor to remove radiant heat transmitted from the reactor. 
Additionally, convective heat transfer between the reactive gas and the coolant is 
also permitted when the reactive gas flows pass the tube banks. 


Propane/air 


Equivalence ratio 
Mixture flowrate 
Preheat temperature 
Pressure 


^ = 0.8 
m = 5.75 g/s 
Tin = 644K 
P = 1 X 105 Pa 


Hydrogen 

Inlet temperature 

Pressure 

Flowrate 


Reactor diameter 
Reactor length 
Void fraction 
Cell diameter 

Coolant tube diameter 
Coolant tube thickness 
Coolant tube bank center 
to center distance 


To = 700K 
P = 100 X 105 Pa 
iric = 1.19 g/s 


D = 5.00 cm 
a = 5.75 cm 
Vf = 0.87 
DTub ~ 9-32 cm 

dc = 0.70 cm 
tc = 0.06 cm 

Ad = 0.80 cm 



SaMPLL 

PKOBLLM E- 

KONOLITH HLfULLN HLAT 

LXCHAfiGEKS 



b 

3 0 3 

0 0 






n.3l7b 

0.87 

5.0 

0.0005 

0.6 




b,75 

b^^ . 

1.0 






7-JO. 

inoo. 







0.0 

0,02 

0.04 

0.06 

0.06 

o.in 

0.15 

0.20 

0,25 

0,30 

0,40 

0.50 

0.70 

1.0 

1.5 

2.0 

2,5 

3. 

3.5 

4. 

4.5 

5.25 




1800, 

1600. 

1600, 

1800. 

1600. 

1800. 

1600. 

1«00, 

1600, 

1600. 

1600. 

1600. 

1600. 

1800, 

1600. 

l^lOO, 

1800. 

1600, 

1800. 

1600. 

1000. 

1600. 

1600. 

1 1 1 

11111 

111111 

11111 

11111 

1 1 1 

1111111 

till 

1333, 

1.165 

700. 

100. 

2. 




n .055 

0,065 

0.7 

0.8 

. 0.1 




fef3. 

1073, 

0.000756 

0,00106 





2.627 

59.7 

6.52 

0.00078 

12000.0 




N2 

.764 

.7401 

1.03475 





02 

.203 

.0394 

1,10045 





C3HB 

.0325 

5. t*0bl.2756 





H20 

1.0 E 

-10 .126 

.644069 





C02 

1.0 L 

-10 .0945 

1.27446 





2 1 
C3HB 02 


C3H6 

02 H20 

C02 

4.1 

£♦9 

40,61 

.25 1, 

04 



-.75 -3.964.0 

3.0 


C3H6 


02 

H20 C02 


1.1 

£♦09 0. 

17,6 

1.0 




•5.0 4.0 

3.0 




00 

cn 



sample PHOHLF-rt K- monolith htTWFEN HEAT LXCHANGEKS 


• INTEGHAL PAHAMLTEHS * 


IS I NUMBER OF SPECIES I = 5 
NL ( GHIO POINTS » = 22 
IHT { NUnUFH OF OVERALL ITERATIONS I = 3 
ll.OSS { HEAT LOSS OPTION » = 0 
iGEOn I PROBLEM GEOMETRY OPTION I = 3 
ICCEFF C MEAT TRANSFER INPUT OPTION ) = 0 
KR7 I PRINT OPTION ) = Q 


• NONINTEGRAL PARAMETERS • 

•** REACTOR data •** 


MONOLITH honeycomb CATALYTIC COMBUSTOR 
NUMBER OF TUBES s 54.94 

VOID FRACTION s ,B70 

diameter of COMBUSTOR s 5.00 ICM> 

EFFECTIVE CELL DIAMETER s .32 CC«) 

00 

or INLET TEMPERATURE 

PRESSURE 
mass flow RATE 
CONDUCTIVITY 
EMISSIVITY 


••• RESERVOIR temperatures 

UPSTREAM = 700,00 

OOUnSTHEAM = 1000,0 


GRID POINTS CHOSEN IS) 


, 0 U P 0 

,4U0n 

4,5no 

,20UU-U01 

,500U 

5,250 

,4000-001 
• 7UUU 

,6000-001 ,6000-001 
1,000 1,500 

,1000<fUOO 

2.000 

• 1500 
2,500 

,2000 

3,000 

• 2500 
3,500 

.3000 

4,000 

INITIAL 
1800 . 
1600, 

GUESSED WALL 
1600 • 
1600, 

TEMPERATURES CTL) 

1600, 1600, IBOO. 

1600, 16UU, 1800, 

1600. 

1600, 

180U. 

18U0. 

1800. 

1600, 

1600, 

1600, 

1600. 

1600. 


IROn , 


s 644,00 IK) 
s 1.00 (ATM) 
s .5754001 IG/S) 

= .500-003 (CAL/S-CM-K) 
: .60 



♦♦♦ COOLANT PHOHEKIIES 


= 700.no (K) 

= 100.00 lATM) 
z ,119^001 (G/SI 
= 1333.00 iCfl/SI 
= 2.00 


COCLANT 

irULT TLMPEHATliHf 
PHESSUhK 
MOSS Fl.UW hate 
APhHOACH \/LLOCnY 
rinLLLULAK Ut IGMT 

FOM coolant inhe: 

HjPE: niAMETLH 
Tubt THICKNESS 
CONUIICTIVIY 


.70 (CM) 

.Obb ICM) 

.550-001 (CAL/S-CM-K) 


CFMCH TO CFNtEH niSTANCE ( TORE BANK) = .BO (CM) 

niSTANCE BETWEEN MONOLITH ANO TUBE BANK = .100 (CM) 


00 


• THERMOCHEMICAL data for coolant 


temperature 

673.00 

1073.00 


Thermal conductivity 
•756-003 
•106-002 


heat capacity CPH = Cl + C2 ♦ T ♦ C3 / T*«2 


Cl = •652f001 (CAL/G-MOLE - K) 

C2 = ./mO-003 (CAL/G-HOLE - K**2> 

C3 = . 120 + 005 (CAL-'K/G-MOLE) 

viscosity parameters 

sigma = 2.027 (A) 

E/K = 59.700 (K) 


(CAL/G-MOLE - K) 



• SPECIES initial WOLC FM ACT IONS ,F IRST GUESSES«ANO DIFFUSION FACTORS • 


SPECIE 

ALPF 

ALPE 

DIFFUSION FAC 

K2 

.764+000 

•740«000 

1.0347 

C2 

.205^000 

.394-001 

1.1004 

C4M8 

,32b-U01 

• SOO-OO** 

1,2758 

h^O 

.100-009 

•126«000 

.8441 

C02 

•lOU-009 

•94S-001 

1.2745 


C30 

CD 



• rHEHMOCHtmSTHY UATA 


CUHVE fit of data 

IN FORM CP 

=RB+RC*T+RU/«T*TI ICAL/(M0LE 

• K» 




HF « CAL/MOL ) 

HiCAL/MOL ) 

RH 

RC 

KO S 

tCAL/MOL/KI 


TU(KI 

NAMC 

7<^UUU4-000 

26 ,uon 
, 23'393>005 

.61391+00 1 

. 16415-002 

.29^274005 

.645464002 

3Q0, 

1000, 

N2 

,7^UUU4^00U 

,221664005 

.030494001 

,23357-003 

-.726064006 

.637744002 

lOOO, 

2500. 

N2 

02 

,uoooo 

32,000 

,255414005 

,670444001 

,16795-002 

-.194194005 

*669064002 

300, 

1000. 

02 

,UOQOO 

,234504005 

,607254001 

,50513-003 

-.213194006 

.679764002 

1000. 

2500. 

02 

H20 

,5779rt*005 

16,000 

,317164005 

,682904001 

.29575-002 

.276434005 

.690224002 

300, 

1000. 

H20 

,5779fl4.005 

,303034005 

.996964001 

,12774-002 

-.143284007 

,664634002 

1000, 

25U0. 

H20 

C02 

,^^40544003 

44,000 

,414364005 

.957784001 

.36020-002 

-.163694006 

•619974002 

300. 

1000. 

C02 

,'#40544003 

,365734005 

.139724002 

.38449-003 

-.139124007 

,798674002 

1000, 

2500. 

C02 

C3H8 

, 246204^003 

44.000 

.108664006 

.320904002 

.11926-001 

-.203724007 

,A59404003 

300, 

2500. 

C3H6 

246204.005 

,108664006 

,320904002 

.11926-001 

-.203724007 

•159404003 

2500, 

6000, 

C3H6 



total NUriHtH OF HCACTIONS 
^UhHEK OF GAS PHASE REACTIONS 


z u o o 

0-1 o o 

o ^ o 

H> C ^ ^0 

< X • • 

> -I o r* 

M «X » ^ 


^ ^ 3 3 

u i/> ^ ^ 

< • o o 

u. c ♦ ♦ 

t_) o o 

(LI CO 

M kJ •-« 

UJ _J 

o • • 

UJ IT 
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AXIAL 

.02U 

nisTANCE 

*040 

ALONG MONOLITH , CyLINHEK OR 
,060 ,060 *100 

PLATE 

.150 

S (CM) 
,200 

.250 

• 300 

,400 

BULK TEMP «K) 

665.0 

677.7 

689,9 

700.7 

710.5 

731.1 

749.4 

766,1 

781.6 

009,3 

wall temp (Kj 

1814,6 

1885.5 

1698.0 

1904,4 

1912.7 

1926,8 

1949.5 

1969,7 

1987,4 

2009.2 

MOLL fraction 

BULK SPECitS 

N2 

,76403 

,76378 

,76356 

.76338 

.76321 

,76266 

.76255 

.76227 

.76202 

,76157 

02 

.20072 

,19896 

.19753 

.19625 

.19511 

,19272 

.19063 

.18875 

.18702 

• 18396 

C3h6 

*03204 

,03170 

,03141 

• 03116 

.03093 

•03046 

•03004 

.02967 

.02932 

• 02672 

H2o 

*00183 

,00317 

.00429 

*00526 

•00615 

.00798 

,00958 

,01103 

.01237 

•01470 

C02 

*00137 

,00236 

,00321 

,00395 

•00461 

.U0596 

,00719 

,00828 

.00927 

• 01102 

WALL SPECIES 

N2 

,74031 

.74031 

.74031 

,74031 

.74031 

,74031 

.74031 

.74031 

.74031 

.74031 

02 

.03925 

.03925 

.03925 

,03925 

.03925 

•03925 

.03925 

.03925 

•03925 

• 03925 

C3HA 

*00000 

,00000 

•oonoo 

,00000 

•00000 

• 00000 

•00000 

.00000 

•00000 

.00000 

H20 

.12596 

.12597 

.12597 

.12597 

.12597 

.12597 

,12597 

,12597 

.12597 

.12597 

C02 

,09447 

,09447 

,09447 

,09447 

.09448 

,09448 

,09446 

,09448 

.09448 

,09448 



I 



AXI'\L 

niSTANCf 

ALONG monolith , lylindek ok 

PLATE 

S (CM» 





• Son 

.700 

1.000 

1.500 

2.000 

2.500 

3.000 

3.500 

4.000 

4.500 

HUL*^ irMM <K» 

834, b 

878.5 

935.1 

1013.4 

1081.3 

1141.8 

1196.8 

1247.4 

1294.1 

1336.2 

WALL TfflP (K) 
MOLc fraction 

2032,1 

2057.5 

2U83.7 

2104.2 

2112.7 

2114.3 

2111.9 

2102.7 

2079.9 

2026.0 

HULK SPECieS 

N2 

.76117 

.76048 

.75961 

.75841 

.75736 

.75641 

.75555 

.75474 

.75397 

.75323 

02 

.18125 

.17656 

.17062 

•16242 

.15527 

.14885 

.14296 

.13746 

.13223 

.12718 

C3hA 

•02816 

.02725 

•02607 

•02444 

*02302 

•02175 

.02058 

.01949 

•01845 

.01745 

H20 

•01680 

.02040 

.02497 

.03128 

.03677 

.04170 

.04623 

•05046 

.05448 

.05837 

C02 

•01260 

.01530 

.01673 

.02346 

.02758 

.03128 

.03467 

.03785 

•04086 

.04378 

WALL SPECIES 

N2 

•74031 

.74031 

.74031 

,74031 

.74031 

.74031 

.74031 

.74031 

.74031 

.74031 

02 

•03S26 

.03925 

.03924 

.03924 

.03924 

.03924 

.03924 

.03924 

•03924 

.03924 

C3h^ 

.oonoo 

.00000 

•00000 

.00000 

•oonoo 

•00000 

.00000 

.00000 

•00000 

.00000 

H20 

.125^17 

.12597 

.12597 

.12597 

.12597 

.12597 

•12597 

.12597 

.12597 

.12597 

C02 

.09448 

.09448 

.09448 

.09448 

•09H48 

•09446 

.09448 

•09448 

•09448 

.09448 



AXIAL UISTaNCE ALONG MONOLITH , CYLINOLK OH PLATE S (CM) 


5.250 


BULK TEMP (K) 

1385.1 

wall temp (K) 

1878.7 

MOLL fraction 


BULK SPECIES 


N2 

•75217 

02 

.12000 

C3HB 

.01603 

H20 

.06389 

C02 

.04792 

WALL SPECIES 


N2 

•74031 

02 

•03924 

C3Ha 

•00000 

H20 

•12597 

C02 

.09448 



• ♦ HAT LXCHAiJUrK PNOHLffl OPTION 2 


cnCLANT templhatuhe; in <K» 
COOLANT TEmPEHaTUMF out IK» 
TOTAL HEAT TRANS^LRELl TO 
COCLANT (CAL/S) 

Trr'PtRATUHf OF FULL/AIR (KH 
UPSTREAM OF HEAT EXCHANGER 

downstream of heat exchanger 


exchanger 1 

700.00 
7oy.51 

•4004002 

G37.93 

649.00 


exchanger 2 

/C9.S1 

/40.i4 

. 131*003 

1385.05 

1325.34 
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SECTION 5 


PROGRAM AND SUBROUTINES 

Brief descriptions of the main program and subroutines are given in 
this section. The code consists of the main program, MAIN!, and eighteen 
subroutines. A cross reference between these routines is given in Table 5-1 
and a simplified flow chart is shown in Figure 5-1. 

5.1 MAINl MAIN PROGRAM 

The main program, MAINl, has a number of functions and serves as the 
driver for the entire program. Primarily, MAINl initializes some constants 
and reads a variety of input parameters. Subroutines READIN, OBTAIN and 
KINKIN are called by MAINl to read in initial species concentration, thermo- 
chemical and kinetic reaction data, respectively. Following this operation, 
the main program loop on grid solution iterations is entered. In this loop, 
space step quantities and terms necessary for the governing equations are 
constructed. The FLAME subroutine is then called to calculate the change in 
concentrations and temperatures due to chemical reactions. Following this 
loop, the solution is checked for convergency, and the loop is repeated until 
the solution is converged. 

BEES Subroutine 

This subroutine sets up the initial terms for solving the 2-D heat 
conduction problem when the traverse heat loss option is invoked. 

COEF Subroutine 

This subroutine calculates the reactor film transfer coefficient based 
on specified catalytic reactor geometries and flowrates. 
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TABLE 5-1. SUBROUTINE CROSS REFERENCE 


Routine/ 

Subroutine 

Calls Routine/ 
Subroutine 

Called by 
Routine/ 
Subroutine 

BEES 

CPROP 

MAINl 

COEF 

— 

MAINl 

C0ND2 

— 

FLAME 

CPROP 


BEES 

FLAME 

0PT2 

FLAME 

CPROP 

C0ND2 

RERAY 

MAINl 

GETDAT 

— 

OBTAIN 

KINKIN 


MAINl 

MAINl 

READIN 

OBTAIN 

KINKIN 

VIEWl 

VIEW 

BEES 

FLAME 

TRIDM 

PROP 

COEF 

STNPRT 

OPT2 


OBTAIN 

GETDAT 

MAINl 

OPT2 

CPROP 

PROP 

MAINl 

PROP 

— 

MAINl 

READIN 

— 

MAINl 

RERAY 

— 

FLAME 

STNPRT 

— 

MAINl 

TRIDM 

TRID 

MAINl 

VIEW 

— 

MAINl 

VIEW! 

-- 

MAINl 
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Figure 5-1 • Program flow chart 








C0ND2 Subroutine 


This subroutine solves the heat conduction term in the surface energy 
balance based on the information obtained in the TRIDM subroutine. 

CPROP Subroutine 

This subroutine calculates the coolant heat transfer coefficient. 

FLAME Subroutine 

This is the program's key subroutine which solves the species and 
energy equations including the effects of chemical kinetics. The routine 
begins by setting up the individual species thermochemical properties. Util- 
izing current solution estimates, errors and their derivatives with respect 
to the variables are formed for the enthalpy, total pressure and species mass 
balance equations. The resulting equation errors are then checked for con- 
vergence. If convergence is not achieved the matrix of error derivatives is 
inverted and multiplied by the errors to obtain the corrections to the vari- 
ables needed to drive the equation errors to zero. If the corrections exceed 
a certain multiple of the variables, all of the corrections are uniformly 
damped so as not to overcorrect the variables. Corrections are then made 
and the solution procedure cycles to the top of the subroutine to repeat the 
above process until either a converged solution is found or the allowed number 
of iterations is exceeded. 

READIN Subroutine 

This subroutine reads the species names, initial species concentrations 
and diffusion factors. 

OBTAIN Subroutine 

This subroutine scans the curve fit thermochemical data file for the 
species of interest. Once the proper species thermochemical data file is 
found, subroutine GETDAT is called to extract and store the data. 
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GETDAT Subroutine 

This subroutine extracts curve fit thermochemical data from input/ 
output unit 11 and stores the data for use during a calculation. 

KINKIN Subroutine 

This subroutine reads in the kinetic reactions and their associated 
rates. It then develops the stoichiometric coefficients for the reactions 
and checks to see if the reactions balance. The stoichiometric coefficients 
are applied in the FLAME subroutine. 

RERAY Subroutine 

This subroutine is a generalized matrix inversion routine which either 
gives a set of solution vectors or the solution vectors plus the full matrix 
inversion. 

0PT2 Subroutine 

This subroutine sets up the computer code to perform the problem with 
a catalytic reactor radiating to up- and downstream heat exchangers. 

PROP Subroutine 

This subroutine calculates the transport properties for the reacting 
gas given the gas composition, temperature, and pressure. 

STNPRT Subroutine 

This subroutine sets up the computer code to print out the results as 
well as the input information of the specified problem. 

TERMS Subroutine 

This subroutine further reduces the terms set up by the BEES subroutine 
into a smaller set. 
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TRID Subroutine 


This subroutine is a general routine which solves a set of tridiagonal 
linear equations. 

TRIDM Subroutine 

This subroutine uses the terms set up for the 2-D heat conduction 
problem by the TERMS subroutine and rearranges them into a set of tri diagonal 
linear equations. It then calls the TRID subroutine to solve the equations 
and uses the subsequent result to calculate reactor wall temperatures. 

VIEW Subroutine 

This subroutine calculates the view factors for monolith reactors given 
the monolith cell geometry. 

VIEWl Subroutine 

This subroutine calculates the view factors for reactors with concentric 
tube configurations. 
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SECTION 6 


SYMBOLS AND COMMON BLOCKS 

Variables are contained in COMMON blocks and INCLUDE statements which 
represent a set of COMMON blocks. The COMMON'S in each INCLUDE statement are 
shown in Table 6-1 and the COMMON blocks and INCLUDE statements in each rou- 
tine are shown in Table 6-2. A list of symbols in these COMMONS blocks and 
INCLUDE statements is given in Table 6-3 and local variables in the routines 
are given in Table 6-4. 
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TABLE 6-1. COMMONS BLOCKS IN INCLUDE STATEMENTS 


INCLUDE 

COMMON Blocks 

BASIC 


ENERGY 

ENERGY, VFA, SAVE 

HEAT 

HEAT, TRANS, PROPS, CHEAT, 
F2, AMB, TWS, BBB 

HEATl 

HEATL, MEAN 

PROCl 

FI, F3 

PR0C2 

F4, F5 

PR0C3 

F6 

PROF 

CARDl, lUNIT, LIN, WANTS 
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TABLE 6-2. COMMON BLOCK CROSS REFERENCE 



INCLUDE 

COMMON 


BASIC 

ENERGY 

H- 

<C 

LU 

3: 

HEATl 

PROCl 

PR0C2 

PR0C3 

PROF 

BLKl 

BLK2 

BLK6 



BLK9 
BLKll 
BLKl 9 

BLK21 

s 

CQ 



BEES 

B 

B 

X 

X 

X 




X 




X 


X 



COEF 

X 


X 

X X 










X 




C0ND2 

D 

■ 

X 

X 

X 








X 





CPROP 

D 

I 

X 

X 





X 



X 

X X 





FLAME 

B 

B 

X 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

GETDAT 

X 



X 


X 

X 











KINKIN 

X 



X 

X 

X 

X 











MAINl 

X 

X 

X 

X X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


X 


X 

OBTAIN 

X 



X 



X 











0PT2 

X 

X 

X 

X 


X 



X 

X 


X 






PROP 

X 


X 

X 

X 

X 

X 











READIN 

X 



X 

X 


X 











RERAY 

X 

















STNPRT 

X 


X 

X 

X 

X 



X 

X 






X 

X 

TRIDM 

X 


X 

X 

X 













VIEW 

X 

















VIEWl 

X 

X 

X 





X 
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TABLE 6-3. SYMBOLS IN COMMON STATEMENTS 


Variable 

Common 

Description 

A"( , ) 

F6 

Basic matrix of coefficients set up in FLAME, 
fully inverted on last iteration in FLAME, and 
used to set up linearized corrections in ACEF. 
Equivalenced to DAA for storage economy only. 

ADSM 

F5 

Composite term used for integration 

AEMl 

F2 

Iteration variable 

AEM2 

F2 

Iteration variable 

AH 

F5 

Not used 

AHH 

F5 

Composite term used for integration 

AK 

F2 

Thermal conductivity of reactor (cal/s/cm-K) 

AKT 

BLK2 

Thermal conductivity for coolant tube 
(cal/s/cm-K) 

ALPF( ) 

F3 

Assigned initial or frozen values of the a. 's 
which are input to the code. These values^ are 
assigned to station one and are never varied. 

ALPHI( ) 

F6 

The current set of a-j's at a given station as 
generated by the FLAME subroutine. Entering 
FLAME these variables contain the terms a. 
introduced in Eq. 44. 

AM 

ENERGY 

Mass flowrate of reactants (gm/s) 

AMH 

BLK2 

Coolant flowrate (gm/s) 

A1 

TRANS 

Interation variable 

A2 

TRANS 

Iteration variable 

A3 

TRANS 

Not used 

A4 

TRANS 

Not used 
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TABLE 6-3. CONTINUED 


Varible 

Common Description 

B ■ 

BLK8 

4 

Radiation variable oT 

BF( ) 

F3 

Not used 

BIN(JS) 

TRANS 

Not used 

BP3(J) 

BBB 


BP5(J) 

BBB 


BP6(J) 

BBB 

Variables defined in Subroutine BEES 

k 

BP7(J) 

BBB 1 


B1(J) 

BBB ' 


B16(J) 

BBB 

i 


C 

y 

F5 

Perimeter of reactor cells (cm) 

CC 

BLK8 

Distance between monolith and tube bank (cm) 

CH 

HEAT 

Nondimensional heat transfer coefficient 

CHH 

F5 

Composite term used for integration 

CHM 

HEAT 

2 

Mass transfer coefficient (gm/cm -s) 

CHT 

BLK21 

Same as CH 

CM 

HEAT 

Nondimensional mass transfer coefificient 

CMH 

HEAT 

LeO.667 

CMSTG 

TRANS 

Not used 

CMT 

BLK21 

Same as CM 

CPF( ) 

F6 

Molal heat capacity of each species (cal/mole/K) 
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TABLE 6-3. CONTINUED 



Variable 

Common 

Description 

CPG( ) 

F6 

Heat capacity of gas mixture in cal/gm-K. 
Within FLAME, CPG is heat capacity times PM, 

CPPH 

BLK2 

Specific heat capacity (cal/gmole-K) . 

Cl 

BLK9 

Specific heat coefficients (Eq. 50) 

C2 

BLK9 

Specific heat coefficients (Eq. 50) 

C3 

BLK9 

Specific heat coefficients (Eq. 50) 

DBTW 

TRANS 

Not used 

DC 

BLK8 

Diameter of combustor (cm) 

DIAM 

F5 

Reactor diameter (cm) 

DQCON 

BLK19 

Conduction term defined in Subroutine C0ND2 

DR 

BLK2 

Coolant tube thickness (cm) 

DS(I) 

F2 

One-half of the distance S(I)-S(I-2) 

DT 

BLK2 

Coolant tube diameter (cm) 

DY 

F2 

Same as DR 

eak( ) 

F4 

Activation energy (Kcal/gmole-k) 

EK 

BLK9 

Viscosity parameter e/k (Eq. 46-48) (K) 

EMIV 

VFA 

Emissivity of reactor 

EMIVT 

BLK2 

Coolant tube emissivity 

EXK( ) 

F4 

Temperature exponent 

F{ . ) 

VFA 

View factor 

FC 

F5 

Not used 


106 




TABLE 6-3. CONTINUED 


Variable 

Common 

Description 

FKF{ 

) 

F4 

Pre-exponential factor (mole, cm, s) 

FL( , 

) 

BLKl 

View factor 

FT 


FI 

Not used 

6KT 


BLK3 

Thermal conductivity (cal/cm-s-K) 

GKTl 


BLK9 

Coolant thermal conductivity at TKl (cal/cm-s-K) 

GKT2 


BLK9 

Coolant thermal conductivity at TK2 (cal/cm-s-K) 

H 


FI 

Enthalpy of gas mixture as input to FLAME. In 
FLAME solution the parameter includes other input 
terms (h°). 

HAMB(J) 

AMB 

External heat transfer coefficient (cal/cm -s-K) 

HEW( 

) 

SAVE 

Enthalpy based on wall composition and bulk 

HHO( 

) 

BK4 

Stored local value of HO 

HI{ ) 


F6 

Molal enthalpy of each species 

HN 


BLK2 

Heat transfer coefficient on the coolant side 
(cal/cm^-s-K) 

HNN( 

) 

BKl 

Convective heat transfer coefficient for external 
cooling (cal/cm^-s-K) 

HO 


BLK2 

Heat transfer coefficient on the hot gas side 
(cal/cm^-s-K) 

HOS( 

) 

F6 

Product of molal enthalpy and partial pressure 
of each species; modified to include heat 
transfer terms for computational efficiency 

HW 


ENERGY 

Enthalpy based on wall composition and temperature 

ICC 


FI 

Value of ITC as carried by program with regard 
to heat loss operations (see Section 5.1) 

ICOEF 

BKl 

Input external heat transfer coefficient flag 
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TABLE 6-3. CONTINUED 


Variable 

Common 

Description 

ICON 

FI 

Flag indicating convergence (=0) or nonconver- 
gence (=1) of FLAME routine 

IGEOM 

FI 

Reactor geometry option 

IKIN( ) 

BLK7 

Noncatalytic (0) or catalytic (1) node 

ILOSS 

FI 

Heat loss option 

INMAS 

lUNIT 

Mass storge unit 

IS 

FI 

Number of species 

ISP 

FI 

IS+1 

ISPECI(2) 

CARDl 

Species name 

ISS 

F5 

Number of species x 2 

ITER 

FI 

Iteration count in FLAME 

ITURB 

HEATl 

Not used 

IWANT(3, ) 

WANTS 

Flat (IWANT ( , ) and species name (IWANT 
( )) of desired species. 

KR 

FI 

Kinetics flag in FLAME 

KR7 

FI 

Diagnostic output option 

L 

FI 

Axial station number 

MGAS 

FI 

Number cF gas phase reactions 

MT 

FI 

Total number of reactions 

N 

FI 

Number of species 

NAMA( ) 

F6 

First four characters of species name 

NAMB( ) 

FI 

Second four characters of species name 
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TABLE 6-3. CONTINUED 


Variable 

Conmon 

Description 

NL 

F2 

Number of grid points 

NWANT 

WANTS 

Number of species in input list. Must equal IS. 

P 

FI 

Pressure in atmospheres 

PHI( , ) 

F4 

Not used 

PLN 

FI 

Log of pressure 

PLP( ) 

F6 

Log of partial pressure of each species 

PM 

F6 

Product of pressure and molecular weight 

PMSl 

SAVE 

Saved value for pressure molecular weight 

PMS2 

SAVE 

Saved value for pressure molecular weight 

PMU( , ) 

F4 

Stoichiometric coefficient of each product 
species in each reaction 

PP( ) 

F6 

Partial pressure of each species 

PPM( ) 

BLK6 

Local value of PM(1 ) 

PR 

PROPS 

Prandtl number 

PS( , ) 

BLK6 

PP(I)/PM(1) 

PT 

BLK8 

Coolant pressure (atm) 

PX 

BLK3 

Not used 

QCHEM 

ENERGY 

Chemical energy term in Eq. (6) 

QCON 

BLKll 

Conduction term defined in Subroutine C0ND2 

QCOND 

ENERGY 

Conduction energy term in Eq. (5) 

QCONV 

ENERGY 

Convective energy term 

QLOSS 

ENERGY 

Total heat loss 
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TABLE 6-3. CONTINUED 


Variable Common 

Description 

QRAD 

ENERGY 

Radiation loss, Eq. ( 8 ) 

QRADIN 

ENERGY 

Radiation from up- and downstream reservoirs 

QR1( ) 

BLKl 

Radiation terms defined in Subroutine VIEWl 

QR2( ) 

BLKl 

Radiation terms defined in Subroutine VIEWl 

RB( ) ^ 

RC( ) 1 

RD( ) 

m ) j 

RF( ) 

RDS( ) ) 

> F 6 ^ 

1 

Curve fit constants for enthalpy, entropy and Cp 
1 of each species. See page 5.14. RB = Fi + 

F 2 , RC = F 3 , RD = F 4 , RE = Fg and RF = Fg 

RMU( . ) 

F4 

Stoichiometric coefficient of each reactant 
species in the reaction 

RU 

PROPS 

Variable defined in MAINl 

RUD 

PROPS 

Variable defined in MAINl 

RT(2) 

F3 

Product of gas constant and temperature 

S(JSl) 

F2 

Grid point locations (cm) 

SB( ) 

F 6 

Entropy/R of each species at pressure (not 
standard state). Later modified to be log Kp 
of formation log p (free energy of formation 
at pressure. 

SG 

BLK9 

Viscosity parameter (Eq. 46-48) 

SNSM 

F5 

Not used 

SPLP( ) 

SAVE 

Saved value for log partial pressure 

SS 

BLKll 

Mid-point between S^ and S^._-j 


no 






TABLE 6-3. CONTINUED 


Variable 

Common 

Description 

ST' 

PROP 

Not used 

SUMH 

F5 

Summation of products of partial pressure and 
enthalpy 

SUMP (2) 

F3 

Summation of partial pressures 

T( ) 

F6 

Temperature (K) 

TAMBI 

AMB 

Not used 

TAMBO 

AMB 

Not used 

TAU 

VFA 

Stefan Boltzman constant 

TCOND 

PROPS 

Thermal conductivity (cal/s-cm-K) 

THMU ( , ) 

F4 

Not used 

TI 

F2 

Inlet temperature (K) 

TKl 

BLK9 

Lower coolant reference temperature (K) 

CVJ 

BLK9 

Upper coolant reference temperature (K) 

1 TL(JSl) 

F2 

Wall temperature (K) 

i TRES(2) 

F2 

Upstream (1) and downstream (2) reservoir 
temperature 

TSQ(t) 

F3 

Temperature to the square power 

TSl 

SAVE 

Saved value for gas temperature 

TS2 

SAVE 

Saved value for wall temperature 

TT( ) 

BK2 

Coolant temperature (K) 

TTI 

BLK2 

Coolant inlet temperature on the coolant side (K) 

TTW( ) 

BLK4 

Stored local value of TTWl 
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TABLE 6-3. CONTINUED 


Variable 

Common 

Description 

TU< ) 

F6 

Temperature at which switch is made from low 
temperature to high temperature thermodynamic 
data curve fits (see page 5-14). 

TW(J) 

TWS 

Iterative value for temperatures 

TWT 

F5 

Not used 

TWl(J) 

TWS 

Iterative value for temperatures 

TW2(J) 

TWS 

Iterative value for temperatures 

TZ 

BLK6 

Average temperature for property calculation 

VG 

CHEAT 

Void fraction of monolith 

Vise 

PROPS 

Viscosity (gm/cm-s) 

VL 

BLK8 

Coolant inlet velocity (cm/sec) 

VNU( , ) 

FI 

Stoichiometric coefficients for equilibrium 
reactions 

W 

BLK6 

Constant =1.0 

WM( ) 

UNIT 

Molecular weight of each species as input 

WT 

BLK9 

Coolant molecular weight (gm/mole) 

X(JS) 

PROPS 

Distance normalized to tube diameter 

XALPF( ) 

WANTS 

ALPF( ) as read 

XHRF 

CARDl 

Heat of formation (cal/gmole) 

XMW 

CARDl 

Molecular weight 

XNU 

PROPS 

Not used 
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TABLE 6-3. CONCLUDED 


Variable 

Common 

Description 

XNUB 

MEAN 

Heat transfer correlation variable 

XNUI 

CHEAT 

Heat transfer correlation variable 

XP 

CHEAT 

Defined in Subroutine COEF 

XPPS 

MEAN 

Not used 

XTP 

CHEAT 

Not used 
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TABLE 6-4. SYMBOLS IN ROUTINES/SUBROUTINES 


Variable 

Subroutine 


Description 

AA 

FLAME 

Inverse of A loaded into AA and corrected for 



rearrangements 

in kinetics package of FLAME 

AAK 

BEES 

Temporary value for reactor thermal conductivity 

AAN( ) 

VIEWl 

Locally defined variable 

AAl 

FLAME 

Locally defined variables for Newton-Raphson 



coefficients 


AA2 

FLAME 



AA3 

FLAME 



AC 

0PT2 

Locally defined constant 

1 

ADSM 

MAINl 

1 

! 

Locally defined variable 1 

AEMN 

MAINl 

Locally defined variable 

AFF( ) 

FLAME 

Affinity of given kinetic reaction (= TAS) j 

i 

AHM( ) 

0PT2 

Same as AMH 

i 

i 

i 

) 

AKK 

MAINl 

j 

i 

Temporary value of AK 

ALFA 

FLAME 

Normalize coolant tube conductance h^Ar/k 

ALPSV( ) 

MAINl 

Saved value of ALPHI( ) 

AM 

MAINl 

Mass flowrate (gm/s) for IGEOM ^ 1 



Mass flux (gm/cm^-s) for IGEOM = 1 

AMU{ ) 

KININ 

y. for reactants 

AN 

VIEWl 

Locally defined variables 

ANK 

CPROP 

Exponential n in Eq. (49) 

ANP 

VIEWl 

Locally defined variables 

AR 

OPT2 

Locally defined constant 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

ARAT 

VIEWl 

Locally defined variable 

AREA 

MAINl 

O 

Cross sectional area (cm ) 

ARG( ) 

TERMS 

Call list function 

AS 

MAINl 

Locally defined variable 

AV 

MAINl 

Surface area to volume ratio (1/cm) 

AVS 

CPROP 

Eq. (46) 

AX 

OPT2 

Locally defined constant 

B 

TRID 

-Tw(I+l) X B6(I+1) - B7(I+1) 

BB 

VIEWl 

Locally defined variable 

BB( ) 

0PT2 

Locally defined variable 

BIG( ) 

FLAME 

Largest positive contribution to each mass 
balance. After kinetics package becomes 
largest absolute contribution. 

BMU( ) 

KININ 

u-j for products 

BRl 

MAINl 

Locally defined variable 

BR2 

MAINl 

Locally defined variable 

BUST 

FLAME 

Value assigned to BUMP 

C 

OPT2 

Locally defined variable 

Cl 

VIEWl 

Locally defined variable 

C2 

VIEWl 

Locally defined variable 

C3 

VIEWl 

Locally defined variable 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

C4 

VIEWl 

Locally defined variable 

C5. 

VIEWl 

Locally defined variable 

CONST 

FLAME 

Locally defined iteration variable 

CPDF 

FLAME 

Local summation 

CPDW 

FLAME 

Local summation 

CPEW 

FLAME 

Local summation 

CPEWW( ) 

OPT2 

Specific heat capacity of edge gases at 
wall 

CPH( ) 

0PT2 

Same as CPPH 

CPP 

OPT2 

Average gas specific heat capacity 
(cal/gmole-K) 

CPTIL 

PROP 

Locally defined variables for thermodynamics 
properties calculations 

CT 

FLAME 

IT times coolant tube diameter 

CT 

0PT2 

Pi times coolant tube diameter 

D( , ) 

VIEW 

Locally defined variable 

DAMP 

fLame 

Damping factor 

DCHM 

FLAME 

Zero 

DCTW 

FLAME 

Zero 

DENS 

CPROP 

Density (gm/cc) 

DERR 

FLAME 

Derivation of error 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

DIA6( ) 

TRID 

-B4 (I+l) 

DIAM 

VIEW 

Same as DTUB, effective cell diameter (cm) 

DIAMC 

VIEWl 

Same as DTUB, effective cell diameter (cm) 

DIAMS 

MAINl 

Diameter/height of reactor (cm) 

DKPT( ) 

FLAME 

The partial of log K with respect to log T for 
each kinetic reaction 

DLNTW 

FLAME 

Temperature correction 

DLTX 

MAINl 

Temporary value for inlet temperature (K) 

DM1 

FLAME 

Locally defined variable 

DM2 

FLAME 

Locally defined variable 

DNE 

FLAME 

Number of reactants in a kinetic reaction ! 

DOM 

FLAME 

Log of mole fraction of a species 

DS 

VIEW 

Locally defined variable 

DSS 

. VIEW 

Locally defined variable 

DTT 

CPROP 

Same as DT, coolant tube diameter (cm) 

DTTWI 

FLAME 

Normalized temperature correction , 

DTUB 

MAINl 

Effective cell diameter (cm) j 

DTW 

COND2 

Local variable 

DTWDT 

FLAME 

Temporary value for hot order wall temperature; 
derivatives for iteration corrections 

DTWDTH 

FLAME 


DTWDTW 

FLAME 


DTW I 

FLAME 

Locally defined variable ; 

DTWII 

FLAME 

Locally defined variable ; 

. . . _ _ . - _ J 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

DUM 

FLAME 

Locally defined temporary values 

DUMl 

FLAME 



DUM2 

FLAME 



DUMN 

FLAME 



DUMQRA 

FLAME 



DY2 

BEES 

2 

0.25 X square of coolant tube thickness (cm ) 

E 

1 

1 

FLAME 

Prior to inversion in FLAME, equal to error 
of mass balance and equilibrium equations; after 
inversion equal to corrections in log partial 
pressures of species 

E 

OPT2 

Locally defined variable 

! 

EMAX 

1 

1 

FLAME 

Maximum correction in species partial pressure 
on a given iteration 

1 

' EMIVV 

MAINl 

Temporary value of emissivity 

1 EMXP 

FLAME 

EMAX from prior iteration 

! 

ERR 

FLAME 

Iteration error 


ESUM 

FLAME 

Locally defined variable 

ETW 

FLAME 

Not used 


ETWW 

FLAME 

Not used 


FCC( ) 

FLAME 

Temporary storage for PP( ) 

FF 

OPT2 

Locally defined variable 

FFl 

0PT2 

Locally defined variable 

FIB 

MAINl 

Locally defined variable 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

FM{ ) 

VIEWl 

Locally defined variable 

FP. 

OPT2 

Local radiation variable 

FP( ) 

VIEWl 

Locally defined variable 

FQ( ) 

VIEWl 

Locally defined variable 

FZ 

FLAME 

Locally defined variable 

6MAX 

OPT2 

Mass flux (gm/cm -s) 

HH 

BEES 

Locally defined variable 

HH( ) 

0PT2 

Same as HN 

H0( ) 

0PT2 

2 

Heat transfer coefficient (cal /cm -s-K) 

HSV 

MAINl 

Locally defined variable 

HT 

FLAME 

Total conductance 

IB 

MAINl 

Iteration control integer 

IBGN 

STNPRT 

Locally defined output counters 

ICB 

FLAME 

Count of iterations with BUMP activated 

ICP 

KININ 

Local counters 

ICR 

KININ 

Local counters 

IDENT 

KININ 

Output alphanumeric 

IDUM 

OBTAIN 

Dummy integer 

lEND 

STNPRT 

Locally defined output counters 

lEOF 

OBTAIN 

End-of-species-file indicator 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

I FLAG 

KININ 

Stop program if there is a problem with kinetic 
data 

IFRZ 

FLAME 

Local control integer 

IHDG 

OBTAIN 

Dumn\y integer 

IMXP 

FLAME 

Index for maximum iterations 

ISL 

FLAME 

Local control integer 

ISPEC( ) 

READIN 

Same as NAMA, NAMB 

ISS 

FLAME 

Number of species plus 2 

ISS 

READIN 

Number of species times 2 

ISW 

KININ 

Flag on whether a species name in a reaction can 
be matched with the species input names 

ITDR 

MAINl 

Iteration counter 

ITEMS 

STNPRT 

Locally defined output counters 

ITERZ 

FLAME 

Interaction counter 

ITMX 

FLAME 

Maximum number of iterations allowed 

IX 

FLAME 

Reaction equation index 

IX 

KININ 

Local index 

J 

READIN 

Locally defined integer 

UMAX 

READIN 

Same as PARAMETER JP (see BASIC PROG) 

KIN 

KININ 

Input unit 

KNY 

READIN 

Unused integer 

KOUT 

KININ 

Output integer 

KR7S 

MAINl 

Stored value of KR7 


120 





TABLE 6-4. CONTINUED 


Variable 

Subroutine 


Description 

KS 

FLAME 

1 or IS+1 

JSS 

FLAME 

IS or 2+IS 

M 

FLAME 

Local index 

MATCH 

OBTAIN 

Local counter 

M6P 

FLAME 

Identity of first surface reaction 

NA( ) 

KININ 

Reactant names 

NASYM( ) 

KININ 

Output alphanumeric 

NB( ) 

KININ 

Product names 

NBSYM( ) 

KININ 

Output alphanumeric 

NEQ 

FLAME 

Total number of equations used in corrector 
step. 

NEQU 

TRIDM 

Number of grid points minus 1 

NI 

0PT2 

Temporary index 

NIT 

OPT2 

Iteration counter 

Nil 

FLARE 

Reray control integers 

NI2 




NI3 




NLL 

MAINl 

Number of grid points minus 1 

NLNGTH 

OBTAIN 

Same as NLNGTH 

NLP 

MAINl 

Number of grid points plus 1 

NMA 

KININ 

Temporary alphanumeric 

NP( ) 

KININ 

Output alphanumeric 

NR( ) 

KININ 

Output alphanumeric 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 


Description 

NT 

OPT2 

Locally defined integer 

NTUB 

MAINl 

Locally defined variable 

NUMNA 

KININ 

Maximum value of ICR for each reaction 

NUMNB 

KININ 

Maximum value of ICP for each reaction 

OMG 

CPROP 

Eq. (47) 


OK 

OBTAIN 

Local counter 

PAGES 

STNPRT 

Locally defined output counter 

PCM 

FLAME 

Locally defined variable 

PI 

VIEWl 

Constant 3. 

1416 

PKP( ) 

FLAME 

Reverse rate of kinetic reaction 

PKPE( ) 

FLAME 

Natural log of PKP 

PKR( ) 

FLAME 

Forward rate of kinetic reaction 

PKRE( ) 

FLAME 

Natural log of PKR 

PLIM 

FLAME 

Maximum allowed change in log partial pressures 

PLN 

MAINl 

An P 


PLPSV( ) 

MAINl 

Saved value 

of log partial pressure 

PMR( ) 

FLAME 

Net forward 

rate of each reaction as used 



in FLAME 


PMSVl 

MAINl 

Locally defined variable 

PMSV2 

MAINl 

Locally defined variable 

PMVl 

PROP 

Locally defined variables for thermodynamics 

PMV2 

PROP 

properties < 

:alculations 

PMV3 

PROP 



PMV5 

PROP 



PMV6 

PROP 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

k 

Description 

i . . . ' . . . 

PP 

VIEWl 

Locally defined variable 

PPSV( ) 

MAINl 

Saved value of partial pressure 

PRMU( , ) 

FLAME 

Net stoichiometric product coefficient for 
each species in each reaction subsequently 
multiplied by CKIN (JE) for computational 
convenience 

PRT 

CPROP 

Prandtl number 

PSUM 

KININ 

Used in reaction balance check 

QC( ) 

OPT2 

Total heat transferred to coolant (cal/s) 

QDIFF 

FLAME 

Diffusion energy 

QQ( ) 

MAINl 

Locally defined variable 

QQ( ) 

STNPRT 

Variable dfined in MAINl 

QR( ) 

FLAME 

Local Radiation variable 

QRAT 

FLAME 

Local radiation variable 

QRR 

FLAME 

Iterative value of radiation loss 

QRR( ) 

0PT2 

Not used 

QRSAVE 

FLAME 

Saved value of radiation flux 

QTOT 

OPT2 

Not used 

QTOTAL 

MAINl 

Total energy transferred to coolant (cal/s) 

R 

VIEW 

Effective cell radius 

RAT( ) 

FLAME 

Maximum of PKP (JX) and PKR (JX). Basis of 
ordering reactions 

RATIO 

TRID 

Locally defined variable 

RE 

COEF 

Reynolds number 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 


Description 

RET 

CPROP 

Reynolds number 

REY 

0PT2 

Reynolds number 

RH 

MAINl 

Locally defined variable 

RR 

CPROP 

Reynolds number per unit length (1/cm) 

RR 

VIEWl 

Ratio of effective cell diameter to coolant 
tube diameter 

RSUM 

KININ 

Used in reaction balance check 

SIGN 

FLAME 

Locally defined variable 

SL 

VIEWl 

Locally defined variable 

SR 

VIEWl 

Locally defined variable 

STEP 

MAINl 

Iteration ste 

p size 

STP( ) 

KININ 

Net stoichiometric coefficient for products 

STPR( ) 

KININ 

Difference between stoichiometric coefficients 
for each species in each reaction equation 

STR( ) 

KININ 

Net stoichiometric coefficient for reactants 

SUB( ) 

TRID 

B3(I+1) 

1 

1 

SUM( ) 

OPT2 

Local summation 

SUMD 

FLAME 

Locally defined summation variables 

SUMHEW 

FLAME 



SUMK 

FLAME 



SUMHW 

FLAME 



SUMR 

FLAME 



SUP 

TRID 

B5(I+1) 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

Til 

MAINl 

Temporary value of T(l) 

TAV 

0PT2 

Temporary value of TAVR 

TAVG( ) 

OPT2 

Locally defined variable 

TC( ) 

FLAME 

-HI (JP)/RT = -d(log Kp)/d(log T) 

TDIF 

MAINl 

Temperature difference control value 

TFZ( ) 

FLAME 

Fixed values of temperature (data statements) 

TG( ) 

OPT2 

Temperature upstream of heat exchanger (®K) 

TG( ) 

STNPRT 

Variable defined in MAINl 

TG( ) 

MAINl 

Locally defined variable 

TGAS 

FLAME 

Temporary value for gas temperature 

TGI 

0PT2 

Temperature of fuel/air mixture (°F) 

TGNLP 

OPT2 

Temperature downstream of heat exchanger (®K) 

TII 

0PT2 

Iterative temperature 

TKE 

CPROP 

00 

• 

cr 

UJ 

TPMU 

PROP 

Locally defined variables for thermodynamics 
properties calculations 

TRESl 

MAINl 

Stored value of upstream reservoir temperature (K) 

TRES2 

MAINl 

Stored value of downstream reservoir temperature (K) 

TSV 

MAINl 

Temporary value for inlet temperature (K) 

TTF( ) 

OPT2 

Coolant temperature out (°K) 
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TABLE 6-4. CONTINUED 


Variable 

Subroutine 

Description 

TTP 


FLAME 

T(l) 

TTRA( 

) 

0PT2 

Locally defined variable 

TTWl 


FLAME 

Temporary value for coolant wall temperature 

TTW2 


FLAME 

Temporary value for hot side wall temperature 

TWSI 


FLAME 

Locally defined variable 

TWSII 


FLAME 

Locally defined variable 

TX 


BEES 

) 

Same as TTI, coolant inlet temperature (K) 1 

i 

f 

i 

1 

U 


0PT2 

1 

i 

1 

! 

V( ) 


FLAME 

! 

Not used i 

j 

VA,VB, 

VD,VE 

VC 

FLAME 

Locally defined variables used in thermo- ! 

dynamics calculations j 

1 

VF 


MAINl 

Void fraction of monolith • 

1 

VFF 


FLAME 

Locally defined variables used in thermo- 
dynamics calculations 

VLK 


FLAME 

Log of (see Eq. (12)) j 

VMU 


PROP 

1 

Locally defined variables for thermodynamics i 

properties calculations ' 

j 

WMM 


MAINl 

[ 

1 

Locally defined variable 

WMl 


MAINl 

Molecular weight of gas at inlet (gm/gmole) 

WM2 


MAINl 

Molecular weight of gas at outlet (gm/gmole) 
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TABLE 6-4. CONCLUDED 


Variable 

Subroutine 

Description 

XA 

OPT2 

Locally defined variable 

XB 

0PT2 

Locally defined variable 

XBP 

READ IN 

Diffusion factors 

XC 

0PT2 

Loca*'ly defined variable 

XF 

READ IN 

Species initial mole fraction 

XE 

READIN 

Species first guesses, mole fractions 

XLE 

COEF 

Lewis number, Pr/Sc 

XNUBM 

MAINl 

Local summation value 

XSUM 

MAINl 

Local summation value 

XX 

FLAME 

Saved values of FKF (JX) prior to pressure and 
third body modifications 

Y 

OPT2 

Locally defined radiation variable 

Z 

OPT2 

Locally defined variable 

Z1 ,Z2,Z3 

FLAME 

Locally defined iteration variables 

ZDS 

VIEWl 

Locally defined variables 

11 

FLAME 

Unused variable 
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SECTION 7 


CODE LISTING 
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SUBROUTINE BEES 
REAL HH 
INCLUDE BASIC 
INCLUDE HEAT 
INCLUDE PRDCl 
INCLUDE PROC2 

C0MM0N/BLK2/TT( JSl) .AKTiDRtDT*ENIVT«AMHiCPPH»TTI.HN«HO 

COMMDN/BLKll/SS 

COMMON/BKl/HNN( JS1» tICOEF 


TX=TTI — 

DT2 = DY • DY / 4.0 
DO 100 I = 2. NL 
SSsCS(I)-Sll-l) )/2.*S(I-U 
IF(ICOEF.NE.l) CALL CPROP(TX) 

AAK=AK 

HNN(I)sHN - 

AK=AKT 
HHSHN 

IFncOEF.EO.DGO TO 1 
GO TO 2 

1 HH-HAIiB(I) 

HNN(1)=HH‘ - 

2 SSN=S{I)-S(1-1) 

C 

Bid) X OY2 / (DSII) « SSN) 

B2(I> X 0Y2 / COSd^’l) * SSN) 

B3d) X AK • OY / 4.0 / OS(I> 

B4d) X SSN • (HH ♦ 2.XAK/0Y) AK*OY/4.*U./DS(I+1)*1./OSII)J 

B5d) X AK • DY / 4.0 / OSd + 1) 

B6d) X 2.*AK / DY • SSN 
B7d) X HH * TX • SSN 
BBd) X 2.*AK * C / DY • SSN 
B9d) X AK • C * DY / (4 . «DS d4-l ) ) 

BIO(I) X AK * C * OT / T4.*0Sdn ' 

BP7d) X B7(I) / B4d) 

BP3d) X B3(I) / B4d) - ■ 

BP5(I) X B5(I) / B4d) 

BP6(I) X B6(I) / B4(I) 

811(1) X 2.0 *■ Bid) ■» B2d) - BPGd) 

BM = 1. - BP3d) 


B12d ) 


BP7d ) 

/ 

BP1 

B13(I ) 

= 

BP5d ) 

/ 

BH 

B14(I ) 

= 

BPGd ) 

/ 

BM - 2. - B2(I) 

Bl5d) 

= 

BP3d) 

/ 

BM 

B16(I ) 

= 

BPGd ) 

/ 

0M - 2. - Bl(l| 


C 

100 CONTINUE 
AK=AAK 
C 

RETURN 

END 
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subroutine coef 
INCLUDE BASIC 
INCLUDE HEAT 
INCLUDE HEATl 
INCLUDE PROCl 
COBHON/BLKai/CHl «CHT 


RE=RUD/VISC 
XLE = PR / SC 
XP=RE*PR/XIL) 

XNUI=.81*PR**.«»*S0RTIRE*«»t/(l»/VF-l.H 
IP(XNUI«LE.3.66I XNUIs 3«66 
CMH=tXLE)**.667 
IF (IGEOM .NE. 1) GO TO 10 

c flat plate laminar flou correlation — 

XNUD s *332 » pR «* *333 • SORllRCI 

GO TO 40 ■ 

10 IF (RE .GT. 10000. » 60 TO 20 

XNUB = 3.66 • (1.0 ♦ .095 • XP • RUD / RU> •• .45 
GO TO 30 

20 XNUB = .023 • RE •• .6 * PR ** .4 
30 IF (L .GT. 2J 60 TO 40 

IF (XNUB .GT. KNUD -XNUB“s XNUl ■ 

40 CONTINUE 

CH=RU*XNUB/PR/RE 

CH=CH*CMH 

CHTsCH 

CMT=CM 

90 FORHAT(ex.l0E10.5) 

RETURN 

END 
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SUBROUTINE C0ND2JNNI _ 

INCLUDE BASIC ® 

INCLUDE HEAT 
INCLUDE PROCl 
include PR0C2 
COMHON/BLK19/OCON,DOCON 
C— — — 

IF(NN .NE. 2» GO TO 100 

TW(NNi s IB12(NN) * B13(NN) » TU2(NN4-1) ♦ TUKNNI * 

% B2(NN) • TWINN^lll / (•B14(NNII • 

“ - TU(NN-1V = TUINN) 

TUl(NN-l)sIWl<NN> 

TU2«NN-1»=TW2(NN> 

OTU=oTWl(NN)/Bl>MNN) 

GO TO 2S0 
100 CONTINUE 

c ■■■■■ ■ ■■■ ■■ - ■ 

IP(NN .NE. NU GO TO 200 

TW(NNI s |B12(NN) * BIS(NN) • TU2(NN-1| * TWl(NN) ♦ 

S Bl(NN) * TU(NN>1) I / <>B16(NNH 

TW(NN-»1) s TW(NNt 
TW1(NN*1)=TW1 (nN) 

TU2(NN+1)=TW2(NN) ■ ■ 

0TWs-TWKNN»/B16<NN» 

GO TO 250 
200 CONTINUE 

C 

TU(NN) : (BPT(NN) *■ TWllNN) ♦ TW(NN-1» • BKNN) 4 
$ TUINN+l) • B2(NN) ♦ TU2(NN-1) • BP3(NN) ♦ 

S TW2(NN'»1) * BP5(NNM / BIKNN) 

0TW=TU1(NN)/B11(NNI 
250 CONTINUE 

C 

TU2(NN) s TW(NN) • <2.0 BKNN) B2<NN)I - TWl(NN) • 

“ » BlINN) * TU<NN-D - B2(NN) • TW(NN4l) 

C 

QC0N=-(Be<NN)'»D9(MN)'fB10(NN) >*TH1<NN)'»BS(NN)*TU(NN)'»B9(NN)«TU1(NN'*’ 

$1)«B10(NN)«TU1(NN-1) ~ 

OOCONs-<B&<NN)-fB9INN)-»B10INN) )*TU1 (NNI'fBOlNNI^DTU 
IF ( NN . EO . 2 ) OOCON=DOCON-»Bl 0 1 NN * *TWl < NN ) 

IF(NN.EO.NL) D0C0N=DQC0N>B9(NN)*TW1<NN) 

C 

RETURN 

END 


131 



on o o no ! n n no 


SUBROUTINE CPROPJTZJ 
INCLUDE BASIC 
INCLUDE HEAT 
INCLUDE PROCl 
COMMON/BLKH/SS 

C0HM0N/BLK2/TT IJSl ) » AKT«OR *0T «EHI VT « AhH*CPPH« TTI «HN«hO 
C0MM0N/BLK6/VL « PT » WT , B » DC « CC # 6K T 
C0MM0N/BLK9/TK1 «TR2*GKT1 •GKT2«»G.EK«C1 tC2tC3 
C*»« VISCOSITT 

TKE=T2/EK 

OHGs1.16/TKE**0.1«»8»0,529/EXPI0.77S*TRE»*2.162/EXP|2,'»38*TKE>- " 
AVS=26.69«SeRT(WT«TZ)/IOn6«SG»S6)*0. 000001 

THERMAL CONUUCTlVITt 

ANKrALOG <GKT2/GKT1 ) /ALOGI TK2/TK1 » 

IFUNK.LT. 0.0)60 TO 7 

GKTsGKT1*JTZ/TK1)**ANK - 

GO TO 8 

7 ANK=-ANK 

6KTsGKTl/JTZ/TKl)**ANK 

heat CAPACITT 

8 CPPHsC1+C2*TZ*C3/TZ**2 “ 

CPPHsCPPH/HT 

•• DENSITT 

OENS=PT/TZ*MT/82,OH 

PRANDTL AND REYNOLDS NUMBER CALCULATIONS 

PRT=AVS*CPPH/6KT 

RET=VL*DENS/AVS 

RR=RET 

RET=DT*RR 

IF(IGE0M.EQ.1)RET=SS*RR _ 

•• HEAT transfer COEFFICIENT CORRELATIONS 
DTTsDT 

ifiigeom.gt.dgo to I 

IFIIGEOM.EO.DGO TO 2 
HN=1.1*PRT**0.333*GKT/DT 
IF(RET.GT.500. )GO TO 3 
HN=HN»(0,48*SORT(BET>+0,«»3) 

GO TO 99 

3 HN=HN*C0.46*SORT(RET)+0.00128*RET) 

GO TO 99 

1 IF(RET,GT.10000)GO TO 9 
HN=3.66*GKT/DTT 

GO TO 99 

4 HN=0.023*GKT/DTT*RET**0.6*PRT**0.4 
GO TO 99 

2 IFIRET.GT.50000IGO TO S 

HN = .332 • GKT / SS * PRT *♦ .333 * SQRT(RET) 

GO TO 99 

5 HN=0.0296*RET**0.8*PRT**0,333*GKT/SS 
99 RETURN 

END 


i 

! 

I 
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SUBROUTINE FLAME 
INCLUDE BASIC 
INCLUDE HEAT 
include ENERGY 
INCLUDE PROCl - ■ - 

INCLUDE PR0C2 
INCLUDE PROCS 

COMMON/BLKl/ FLIJS1« JSI) tORlI JSII (OR?* JSII 

C 0 MM 0 N/BLK 2 /TT»JS 1 > .AKT.OR«OT.EHivT.AHHtCPPH.TTItHNtHO 
C0MM0N/BLK6/PS( JSt JP2I «PPM< JSX» tW.TZ 

■ COMMON/BLKT/IKINt JSl) 

C0MM0N/BLKe/VL<PT«UT«B«0CfCC«6KT 
COMMON/BLK21/CHT »CHT 
C0MH0N/BLK19/0C0N,D0C0N 
COMHON/BKl/HNN( JSl ) « ICOEF 
COMMON/BK3/PX( JS«JP2I 

COMMON/BK*/HHO(JS1) .TTWJJSir 

DIMENSION ORIJSl) 

DIMENSION TC < JP2 ) , E < JP2 » » FCC J JP2 1 , XX ( JP2 ) t PRMU ( JP2 • JX I « AFF « JP2 » 
$ t DKPT ( JP2 ) « PKP ( JP2 ) • PKR I JP2 ) • PMH ( JP2 1 « RAT I JP2 ) • BIG I JP2 1 
StPKPE( JP2) .PKRE<JP2» 

■ DIMENSION VtJYI *AA(JT.Jt) ' ' 

DIMENSION TF2I9I 

DATA TFZ/1900. « 1300. t 1200.* 1100., 1000. « 900. «800*«700.« 600./ 
DATA NUL/0/«ICOLD/0/ 


NSSsISS 

DO 1 J=1,ISS 

PP(J)=EXP<PLP(J>) 

FCC(J)=PP(J> 

1 SPLPU)=PLPtJ) 

PMS1=PH(1) 

PMS2=PH(2) ■ ■ 

TS1=T(1) 

TS2=T»2) 

IFRZ=0 

ICB=0 

BUST=100. 

ITHX=90 ■ 

IHXP=1 

PLIM=2.302b65l 

TTP=T(1) 

LMXP=1.E+10 

ITER=1 

DAMP=1. 

IC0N=1 
ISL = 0 
N=ISS+4 

IFCIGEOH .EO. 2» tJ = N ♦ 1 
NEGlsN 

30 DO 2 K=l,2 
CRG(K)=0. 

«/A=ALOG(T(K )/3000. )/l .9B65 
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VP=T<K»-3000, 

VC=IT«K>+3000,)*.5 
vn=T«K )*3000, 

Vt=VC/J VO*VD) 

VFF=VD/1.98bb 

RT(K)=1.9665*T(K) 

TSC«K)=T(K)*T(K) 

SUhP(Kt=0. 

KS=1-MK-1 )*1S 
KSS=IS+(K-1I*IS 
DO 2 I=KStKSS 

H=I-CK-1I*1S 

J=2 '■ ■■ 

IF(T(K).LT.TU(H«X>) J=1 

CPF ( I J=RC (H* J)+T(K ) •RDCHt J>+RE(Ht J)/TSQ(K I 
HK I)=RBCHt J)+VB*(RC(H« J)+RD(H i J)*VC+RECH t J»/VD» 

SBC I )=rF(H» JI+RCIH t J1*VA+VFF*(R0|M« J)*RE(H» JI*VEI"PLP( II 
CPGCK)=CPGCK)-fPP(I)*CPFCll 

PS(L*I)=PPCI)/PM(1)' ■ 

PXCLtI)=PP(II/P 

PPH(L)=PM(1I 

SUHP(K)=SUHP(K)*PP(11 

TCCII=-HI(I)/RT«K) 

ECll=TCm+SB(II_ 

“ SB(I)=Em 

2 OCHMsO. 

CHM=CHt*CH 
DCTW=0.0 
DOll 1 = 1. N 
V(II=0. 

DOll J=liN 

11 A(I.J)=0. 

SUnHW=0. 

SURH=0. 

SHHEU=0. 

DO 3 1=1. IS _ 

HOS(I)=HIlTI*PP(l) 

SUMH=SUNH+HOS ( H 
HOSCI+IS)=Hl(I+ISI*PP(I*ISI 
SUHHW=SUMHU^HOS ( I+IS I 
MEUCI >=H1(1+ISI*PPCI) 

smhem=smhew+hew ( 1 1 
IF (KR7 .UE. 3) GOTO '3' 
hRITECb. 10330) L 

10330 FORMATI//20X.7HSTATION .13) 

UPITEC6. 10000 ) HOSm .HOS(I*IS) .HEW(I) tSUMH.SUHHU.SMHEU 

3 continue 

10000 FORHATCZBX.AEIO.B) 

IF (L .EO. 1) RETURN 

lFClLOSS,Etl.O.AND.I6EOM.EO.O)GO TO 60 

IFCICOEF.EO.DGO TO 201 

TZ=TTILI 

IFIITER.EQ.I) call CPR0PIT2) 

60 CONTINUE 

201 continue 

IF CL .EO. 2) TGAS = TI 
CPDU=0. 
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CPDF=0, 

CPtW=0. 

DO 4 I=1«1S 

CPEW=CPLW*PP(n*CPF«I^IS» 

CPDU=CP0W*PPCI*1S»*CPF(1+1S» 

4 CPDF=CPOF+PP(I)*CPF(l> 

RHS ( 1 ) =-Pn 1 1 ) CHR* ( SURH>SHHCH>CRH* ( SRHEU-SURHU*PR 1 1 1 /PR ( a )l 

• ) 

A ( 1 1 1 > s- 1 1 . ♦CHM ) *T ( 1 ) *CPOF 

A(l.l)sAn«l)>( SUHH>SRHCU4^CRH« ( SRHEW«SUHHU*PR ( 1 1 /PR ( 2 ) )) «OCHR 
All t2» = -CHM*T<2»*( (CRH-1. »*CPEW-CRH*CP0W*PRm/PR<2) ) 

♦ -OCTH* « SURH-SMHEU+CRH* C SRHEW-SURHW*PH ( 1 »/PR 12 » » » 

Aa,3>=PM»l)*H ♦CHH*SUHHH*PRai/PR|2» 

A(1.4»=-PRm /PR ( 2 ) *CHR*SURHW 
00 5 1=1. IS 

A < 1 . 1 +4 ) =-HOS ( n -CHR* ( HOS ( 1 1 -HEH ( n * ( 1 . -CRH ) ) 

5 A(1.I>IS'«-4)=CHR*H0S(I + IS»»PR(1»/PR(2J 
DURQKA=ERIV*(l.-F(L«Ln*TAU*T(^»**4 
0RAT=IQf(A01N'f0UH0RA)*AHH*C ' " 

IFCILOSS .LT. 21 60 TO 32 

OCON=0.0 

DOCONsO.O 
00 31 I=1«N 

31 A(2tlls0.0 

A(2.2)=1.0 

RHS(2)=0.0 
60 TO 3^ 

32 IF (1LOSS.6T.OJ GO TO 33 
IF(L-2) 7.6.7 

6 TWSIs CHH*(TL(L-»1)-T(2J J/OS(L+U 
TUSI I=AER1* ( T ( 2 ) **4 . -TRES ( 1 ) **4 . ) 

DTWI=-CHH*T(2)/0S(L*1I 

0TUH=4.*AER1*T(2)**4. 

GO TO 10 

7 1F(L-NLJ 9.(5. 8 

6 TWSIsAER 2*(T(2)**4.-TRES(2)**4.I 
TUS1I=CHH*(T(2)-TL(L-1 J l/OS(LJ 
0TUr=4.*AEM2*T(2)**4. 

0TWII=CHH*T(2)/0S(L) 

GO TO 10 

9 TUSI= CHH*(TL(L+1)-T(2) )/OS(L+l) 

TWSII=CHH»(T(2)-TL(L-1 ) )/OS(LI 
DTWI=-CHH»T(2)/DS(L+1) 

DTWII=CHh*T(2)/0S(L) 

10 GC0N=TWSI-THSII 
DOCON=OTWI - DTWII 
GO TO 34 

33 TWKL) = T(2) 

CALL C0N02(L) . - 

34 KhS ( 2 ) =PM ( 1 1 •H-SURH+PR ( 1 ) •QCON-PR ( 1 1 *0RAT 
A ( 2 . 2 ) =-PM ( 1 ) • ( 00C0N-AHH*C»4 . *UURO>«A ) 

A ( 2 • 3 ) =-PH ( 1 ) * ( H4^0C0N-QRAT I 
A(2,1)=T(1)*CP0F 
39 continue 

DO 410 1 = 1. IS -- -- - 

A(2.4+I) = PP(I)*HI(I) 

A(3.4-.I)=PP(1) 
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410 Aj4f4*IS+I)=PP{I^IS) 

411 RHS(3)=P-SUMP(1) 

RHS(4)=P-SUMP(2) 

DO 720 IsliIS 

At4 + I,l)=A(4 + Itl) + «PP(n-PPII + IS)*PM(U/PN(2) »*DCHH*CMH 
A(4+^I.2)=ACt+It2) + (PP<H-PPn + IS)*PM(l)/PM(2» )*DCTW*CMH 
A(4*It3)=-PH(l)*ALPHUn-CH«*PPI I'fiS)*PH(l)/PH(2)*CMH 
A(I+4,4)=CHM*PP<I+lS>*Pf1(l)/PH(2»*CMH 
A«4 + I.4 + ll=PPtn+CHM*PP(I)*CMH 

A(4'Ht4+IS+I)=-CHH*CMH *PP ( I+IS ) *PM J 1 > /PM< 2) 

720 RHS(4-H)s*PH|l)*ALPHnl»-PP«I)-CHM*tPP(I»-PP(I + IS)*PM(l)/PM(2))* 
•CHH _ . 

DO 420 Isl.lS 

A(4+IS^I t3)=-PH(l )/PH(2l*PP( I*IS» 
A(4+1S+I.4I=PN(1)/PM(2)*PP<I+1S) 

A(4+IS+I»4+I)SPP(I) 

A(4+1S+1 «4+lS+l )="PP( 1+IS)*PH(1 »/PH(2) 

420 RHS < 44 IS+I ) =-PP ( I » *PP ( 1 ♦! S ) *PN ( 1 > /PH ( 2 ) 


IFdGCOn .NC. 2) 60 TO 422 

CT=3. 1416*01 

7T<1)=TTX 

IFdTER.EO.l) HO=CHT*CPG«l)/PH(l) 

ALFA=HN*DR/IAKT> 

HT=HO-AKT/DR*(l./(l+AtFA7-T.T~ * 

Zl=CT*HO/AM*PM(l)*IS(L)“S(L“l) > 

22=AKT»ALFA/(DR*(1,+ALFA) ) 

23=CT*HN/AHH/CPPH*(S(U-S(fin 

CRRsO.O 

DO 104 Ks2»NL 

IFIK.EO.DGO To 105' “ 

0RH=QRR*FL(K.L)*TL(K)**4 
60 TO 104 

105 aRR=URR+FL(K.U*T(2)**4 
104 CONTINUE 
EHlVTsENlV 

ClR(L)s-CMIVT*TAU*<QRR+QRl(U*TRESd)**4+0R2<U*TRES(2»**4t 

C0NST=H0*T(1)+Z2*TT(L)-0R(L» 

TTU1=TT(L) 

ITERZ=1 

107 ERR=tMIVT*TAU*TTUl**4+HT*TTUl-C0NST 
DEPR=4.*EHIVT»TAU*TTW1**3+HT 

DLNTW=-ERR/DERR 

TTU1=TTW1*DLNTW 

DTTW1=DLNTW/TTW1 

IF(ABS(DTTWl) .LE..002) GO TO 109 
ITERZ=ITERZ+1 
IFdTERZ.LE.SO) GO TO 107 
WRITE(6«10ei TRRtDLNTW 

IDfl F OHMAT ( //lOX «• COOL ANT TUBE ENERGY BALANCE EXCEEDS ITERATION LIMIT 
»ERkOK= • .ElO.bdX, ‘OLNlWr »fE10.5) 

STOP 

109 CONTINUE 

C0NST=4 .*EMI VT*TAw*TTU1**3+HT 
TTW2=( ALFA*TT (L I ♦TTWl ) / ( 1 . +ALFA J 
UTUllT=hn*T ( D/CONST 

L)TWDTiJ = 4.»t.NI\/T*TAl)*FL(LtL)*T(2)**4/C0rpST 
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OTWDTH=Za*TT (L » /CONST 
AA1=Z1*(TI1)-DTW0T» 

AA2=21*DTWUTW 

AA3=21*IT(1>-TTU1) 

A(1,1)=A<1.1)-AA1 
A(l«2)=A(lt2)4^AA2 
A(l«3)=A(lt3)>AA3 ' 

A(l,N»=21*ttTWDTH 

HHSU)=RHSm+AA3 

A(2.1)=A(2«l)<fAAl 

A(2«2)=A(2«2)>AA2 

A(2«3)=A(2«3)4^AA3 

A(2,N)=-A(1.N) “ 

RHS(2)=RHS«2»-AA3 

A(N,1)=-Z3*0TUDT/(1.+ALFA» 

A(N.2)=-Z3*0TW0TW/(l,+ALFAt 

A(N.N)=TTJU*Z3/(1.+ALFA»*(TT<L»-0TW0TH» 

RHS<N)=-TT(L>*TT<L-1»-23*<TT<L)-TTW2» 

HHO(L)sHO - — 

TTW(L»=TTW1 

ORSAVC=OR(l.> 

0H2=AMH*CPPH*CTHL)*TTJL»1» > 

«*?2 CONTINUE 

IF CKR7 .LE. 31 BO'TO 295 

klRITE(6<502l 
D0901 I=1«N 

501 WRITt(6«500) (A( I • J) • J=1«N> < RHS(I) 

500 FORMAT ( IX. 12E10. H/ « 2X. 12E1 0. ) 

502 FORMAT(1X.20HA(I.J) «KHS( I) .BEFORE > 

255 IFtKR-5) 33*» ,333 .33<* ' 

333 00 332 Msl.MT 
IX s 0. 

IFCM.GT.HGAS) 1X=IS 

ONEsO.O 

DO 327 1 = 1. IS 

327 ONE=DNE-RMU(IX+I.H) " 

FKF(n)=FKF(H)*e2.05G«»ONE 
EXK(M>=EXKIM)-»DNE 
xx(n)=FKF(ni 
332 CONTINUE 
KR=<» 

MGPrMGAS+1 

33H DO 337 MsHGP.HT 
337 FKFCM)=XX«M»*PM«1I/CH 
IFCICB.LE.l . BUHP=1. 

DO 340 M=1.HT 

suno=o. 

SUHK=0. 

SUHR=0. 

K=1 

IF(N.GT.HGAS) K=2 
Do 315 1=1. ISS 
PRMU(I.M)=PmMl .M)-RRU( I.H) 
SUHK=SUnK-fPRMU « I . H ) *SB (II 
SUMR=SUMR-.RhU ( I .H » »PLP ( 1 > 

SUnD=SUHU+PHMU< I .H ) *HI ( 1 » 
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PRHU(I,H)=PRMU<ItM)*BUHP 
315 CONTINUE 

AFF<HI=SUrK*RT(K> 

DKPT(M)=SUMD/RT(K) 

VLK=AtOG(FKF(M) ) +EXK ( M » *ALOG( T 1 K ) I -EAK ( M ) /RT ( K I 

IF(K.EO.l) VLK=ALOG(FKF<M) ) ♦EXM h ) *ALOG « TGAS » -EAK ( H> /1.986b/T6AS 
pkpecm)=amini(suhr-sumk+vlk«bo. » 

PKRE<H)=AHlNl(SUnR+VLK*eO, ) 

PKP(M)=EXP(PKPE<M) » 

PKR(m=EXP(PKRr(H) J 
PHRCrDsPKRCm-PKPCHl 
RAT(M)=AMAXl(PKP(f1) *PKR<Hn 
IF (KR7 .LT. 3) GO TO 3A0" 

URITE(6«903) 

WRITE(6i 902) (PRMU( ItM) «1=1*ISS> 

WRITE (6. 341) M.FKFdl) .SUNK iSUHR «DKPT (H) iVLKtEXK(H) ,EAK(M) f 
• PKR)M) «PKP(H) «PMR(H) 

902 F0RMAT(1X«12E10.3| _ _ „ 

903 format ( lHOt49HPRMU.HVFKF.SU«K.SUHRiOKPTiVLK,EXK,EAKtPKR.PKPiPHR /) 
341 F0RHAT(I3«2X«11E11.5> 

340 CONTINUE 

DO 379 Msl.MT 
SIGN=IKIN(L) 

IF(M.LE.MGAS) SIGNs-PM(X)*ADSM 
- DO 375 Isl.ISS 

SUMD=RMU(I.H)*PKRtM)-PMU(I«M)*PKP<M) 

DO 375 Jsl.ISS 

375 A(J*4,I+4) = A(J-»4,l + 4)+SUMD»PRMU(J.M)*SIGN 

SUf1D=PMR(M)*(EXK(M)4EAK(M)/RT(2) )4 PKP«M)*DkPT(H) 

DO 360 Isl«lS 

IF (M.GT.MGAS) GO TO 328 

DUM1=PM ( 1 ) •ADSH*PMR (Ml* PRMU( I *M ) 

A( 1+4,1 l = A( I+4.1 )-DUM1*PKP(M)*DKPT(M) 

A(I+4,3)=A(I+4,3I-0UH1 
RHS(I+4)=RHS(I+4)+DUMl 
GO TO 380 . . 

328 DUHISPMR (M) »PRMU( I+IS,H)*1KIN(L • 

A ( I + IS+4,1 )=A (I + IS+4,1 )-DUf11*DCHH/CHM 

A<I+lS+4,2l=A(I+lS+4,2l+SUMU*PRMU(I+IS.M)*IKIN(U-DUni*DCTW/CHM 
A( 1+1S+4,3)=A( I+IS+4,3l+DUMX 
PHS( I+IS+4 )=RHS( I+IS+4I-DUMX 
380 continue 
379 CONTINUE 

IFdGEOM ,Nt. II GO TO 498 
FLAT PLATE PROHLEH 

RhS( 2) = RHS(l) ♦ RHS(2) 

DO 50 1=1, NEO 

A(2,I) = A(2,I) + A(X,I) 

A(1,I ) = 0.0 
A(3,I)=0.0 
A(I,1)=0. 

A(I ,3)=0. 
fSO CONTINUE 

RHS(l) = o.n 

AC1,1) = 1 .0 

RHS(3)=0.0 

A(3,3)=1.0 
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00 S2 1=1. IS 

no M J=1«NFQ 
AII44.J) s 0.0 
A( J.I-f4) = 0. 

51 CONTINUE 

A(I+4,I*4) = 1,0 

RHS(I«4l s 0.0 

52 CONTINUE 
C 

490 continue 

c 

IF (KR7 ,LE. 3) GO TO 905“ 

HRITE<6,499) 

499 FORMATC1H0.93HA(I,J>«RHS(1). AFTER KINET AND BEFORE INVERT) 
DO 904 1 = 1. N 

904 URITE(6.S00> lA(I.J) .J=1«N) .RHS( 1) 

905 CONTINUE 
FZ=1.0 

IF(IFRZ.EQ.O) 60 TO 711** 

ETWW=RHS«2) 

RHS(2)=0. 

FZ=100. 

IFUFRZ.lt. 0) FZ=1.0 

7114 TF<BUMP.GT.l.) GO TO 101 
DUM=T<1)/A(1.D 
0U)U=T(1)/A)2,1) 

IF (ICOEF .EO. 1) OUni s T <2) / A(2.2) 
IFCaBS(RHSC1)*0UH).GT..2) go to 101 
1F(ABS(RHS«2)*0UM1) .GT. .1*FZ) GO TO 101 
IF(ABS<RHS(3) )/P.6T..00001*FZ) GO TO 101 
1F(ABS<RHS(4))/P.GT..00001*FZ) go to 101 
C 

lF<IGEOn .NE, 2) GO TO 322 
DU«N=TT(L)/A(N.N) 

IF(ABS(OUNN«RHSIN)I.GT.0.2)GO to 101 
322 CONTINUE 
ESUM=0. 

DO 301 1=1. ISS 
PCH=PHU)*ALPHI(I) 

IFU ,GT. IS) PCM=AMAXUPPJ1-IS).PMC1)*PPU)/PM(2) ) 

BIG(1) = AhAX1(PCH*PPU)) 

301 ESUH=AHAX1(ESUM.ABS(RHSU'»4) )/(BIGU)'*’l.E>30) ) 
lF(ABS(ESUn).GT..00001*FZ) GO TO 101 
ICON=0 

IFUFRZ.LE.O) go TO 7115 

ICON=l 

IFRZ=0 

ITER=0 

7115 DO 102 1=1. N 
DO 102 J=1.N 

102 AAU.J)=A(I.J) 

101 NE0=ISS-»4 

IFUGEon .EQ. 2) NEO = NEO ♦ 1 
IF<IFRZ.EU.0) GO TO 7113 
no 7111 1=1. N 

J=2 
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7111 A(J«I)=0. 

J=2 

A(J«J)sl.O 

7112 RHS<J)=0. 

7113 CONTINUE 
NIlsO 

N12=0 

NI3=0 

CALL RERAY(NEO«A«Nll*RHS«l«N12tNI3> 

ETW=-1./«A(2.2)*T|2)I 
00 190 Is5«N 
190 Ea-4»=RHSa) 

ISSsNSS 

ir(KR7.LC*3> GO TO 90 
URITE(6.S09I 
DO 503 lel.N 

503 WRITE(6«500)( A(l« J) f JsltN) «RHS(1 > 

509 F0RnAT(lH0tl9HA(l«J)»RHS(II«AFTCR ) 

90 DAHP=1. ' 

EMAXs 1.E«30 
DO 15 Isl.ISS 

lFJABS(EHAX».6T.ABS«Em»»G0 TO 92 
ERAXsEd) 

IMAX=I 

92 OOHsPLPCD-PCW 
IF(E(I»» 20.15,18 

18 0UH=(9.+9,*00m/C3,-00M» 

GO TO 13 — 

20 0UM=(3.*00M»9, )/l9,+D0M> 

IF(OUM) 13,15,15 

13 CONTINUE 

19 oARPsAHINKOARP, (OUM-OOM)/eU > > 

15 continue 

IF (EdMXPJ/EHXP.LT. -0,331 PLIRsPLIR/2, 

IF (ROD (ITER, 2 » ,EO.O » ERXP=EHXP/tXP(PLIH» 

IF(ABS(1.-E(IRxP>/EHXP) ,LT. .25) PLIRsARINl |PLIR*2. ,9.2103909 ) 

ERXP=ERAX 

IHXPsIRAX 

UARP=ARIN1(DARP,PLIR/( ABS(ERXP) I ) 

0R1=ARAX1(ABS(RHS(1 I ) ,ABS(RHS(2) I ) 

IFdGEOR .EG. 2) OHl = AHAXl (DRl , ABS(RHS(N) I I 
DH2=ARAX1(ABS(RHS(3) ) ,ABS(RHS(9I ) ) 

DARPrAMINl (OAHP,0.9/ARAX1 (DR1,DH2) I 
IF(DARP.GT..99> GO TO 21 
DO 27 1=1, ISS 
27 F(I)=DARP*E(I) 

C HAKE CORRECTIONS 

21 DO 28 1=1, ISS 
PLPd)=PLPd)+E(I) 

PP(I)=EXP(PLP(1)) 

2fl CONTINUE 

IFdGEOR .EQ. 2) TT(L) = TT(L) / d . 0-RHS ( N) •DAHPI 
DO bUO K=l,2 

T(K)=T(KI/(1.-RHS(K)*DAHP) 
bOO PH(K>=PM(K)*(1.+RHS(K+2)*DAHP) 

14(j1 0C0NV=AM*(SUMH/PR(1 )-U) 

UDIFF=r.HM*Art*(SUMH-SRHE(,)/PR(l ) 
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OCHFM=CHM*CMH*AM* J SMHEW/PH C 1 1 -SUHHW/PH « 2 ) ) 

OCONO=AH*awSI-TWSin 

0HA0=»0RAT*AM 

IF (KR7 .LT. 3) GO TO «*8 

URITE.(G*50SI 0AnP«TUI«Phllt«T(2l*Pn(2) 

!>05 FORHATdH .GHOAHP= Ts ttlO.St SH PMs ,E10.!>«HHT2s lElO.S 

•.5HPH2= tCIO.SI 

uRlTE(6t506l (PLPII) «Isl«lSS) 

506 FORHAT(5X.10E10.5» 

WR ITE ( 6 « 510 ) OCONV • QDIFF < OCHEMt UCONU • ORAO 
510 F0HMAT(8 X. ‘EMErGY TERHS»/10X, •CONVECTION a*EX0.5t2Xt 'DIFFUSION s» 
•E10.5.2X,»OChEH =•« 

•E10,3t2Xt»BE0 CONDUCTION =* «E10.5«2X« «WALL RADIATION s'lElO.SV 
WRITE (6.512) OH2.TT(U .TTW1.TTW2.AMH.CPPH.H0.HN 
WKITE(6.512)DT.0R.AKT.HT.ALFA.CT.QRSAVE 
512 FORMAT(5X.8E10.5) 

48 IF (ICON .GT. 0) GO TO 45 

XF(XFKZ.LT.0.AN0.L.NE.2) go to 1721 

GO TO 1708 

1721 T(2) = TF2(-1FR2) 

IFRZ s IFRZ + 1 
ITER =0 
ICON = 1 
GO TO 30 

1708 DO 55 Isi.IS ' — ' 

ALPHI(I*IS)=PP(I*IS)/PM(2) 

55 ALPHHI>=PP«I»/PH(1» 

1414 continue 

H=SUHH/PH(1 I 
HW=SUMHW/PM(2) 

0L0SS=0L0SS^0C0N0+0RAD 

TGAS=T(1) 

RETURN 

45 ITERsiTER+1 

C 

IF(ITER.GT.ITMX-lO) KR7=4 
IF(ITER.LT.ITMX) GO TO 30 
WRITE(6.47) 

47 FORMAT(1H0.27HALLOWEO ITERATIONS EXCEEDED) 

4700 IF(IFHZ.GT.O) GO TO 46 
IFRZ=1 
ITER=0 
T(1)=TS1 
T(2)=TL(L) 

IF(ISL.GT.2) GO TO 46 
IF(ISL.EO.l) T(2)=TL(NL) 

IF(ISL.EQ,2) T(2)=TL(2) 

ISL=ISL*1 
PM(l)=PliSl 
HH(?)=PMS2 
UO 1413 1=1. IS 
PP(I)=FCC(I) 

PP( I+IS)=PP( I ) 

PLP(I + IS)=SPI.P(I) 

1413 PLPd ) = SPLP«I ) 

60 TO 30 

46 CONTINUE 
RETURN 
END 
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SUBROUTINE GETDAHJ) 

INCLUDE BASIC 
INCLUDE PROF 
INCLUDE PKOCl 

INCLUDE PROCS 

C 

DIMENSION 22(JP»2) 

C 

NAMA(J) = ISPECKl) 

NAMB(J) s ISPECK2) 

ALPFIJ) = XALPF(J) 

ALPF(J4IS> = XALPFIJ'flS) 

C READ CURVE FIT DATA 

WM(J)=XMW 

DO 220 Ksli2 

READ (ZNMAS«20220I RAt R9(J«K)i RC(J«K)« RO(JtK)< RE(JiK) 

• RF(J«K). Z2(JiK)i TUld.KI 

20220 FORMAT (6E9.6i 2F6.0) 

220 continue 

WRITE (6«20m0) NAMA(J)t NAMB ( J 1 1 WM ( J) 

20140 FORMAT ( lOX t 2A4 i2 (2X «F10.3 > ) 

DO 240 Kzl«2 

WRITE (S«202RO) ftAi RBIJiKli RC(J<K)f RD(JfKt< RE(J«K) 

* iRFIJtK)* 22(J«K)< TU(J«Kl« NAMA( J) i NAMB(J) 

20240 FORMATfSXt 6(El2»5«lX)« 2(F6«0«lX>t 2Xt 2A4) 

RB(J«K) = RB(J«K) 4 RA 
RF(J«K) s RF(J«K) / 1.9665 
240 continue 

300 RETURN 

END 
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on o on on on 


SUBROUTINE KININ 

integer blank 

INCLUDE BASIC 
INCLUDE PROF _ 

INCLUDE PROCI 

INCLUDE PR0C2 
INCLUDE PR0C3 
C 

DIMENSION Na(S)« NB(S)« AMUCS)* BNUISI 
dimension NASYMlSIt NBSTM(S> 

DIMENSION NR<5».NP<5) *STR(5) «STP<5)tSTPRJJPI 
C 

DATA BLANK /4H / 

DATA NASYM / «*H « 4*«»H ♦ / 

DATA NBSYM / 4h— = ♦ *»*4H f 

DATA NGAS/4HGAS / 

" ■ DATA NSUR/4HSURF/~ 


KIN = 5 ‘ 

KOUT = 6 

read total no. of reactions'and no. of gas phase reactions 

100 HEAd(KIN«20100| mt. mgas 

WRITECKOUT. 201101 MT. MGAS 
URITEJKOUT. 20120) 

IFLAG s 0 

DO FOR EACH REACTION 

DO 400 M = 1. MT 
ISU = 0 
RSUM = 0. 

PSUM = 0. 

' rEAD(KIN. 20400) NA«NB.FKF(M)iEXK|N)«EAK|M) 

REAOtKIN. 20410) AMU. BMU 


IX = 0 

IF(M .GT. MGAS) IX = IS 

c 

DO 300 I = 1, IS 
NMA = NAMA(I) 

THMUJI.M) = 0. 

RMU (I.M) = 0. 

PMU (I.M) a 0. 

■ RMUd+IS.MI =“0. 

PMU(I-.IS.H) = 0. 

C 

STPR(I)=0. 

DO 200 J = 1. 5 
C 

IF(NA(J).NE.NMA) go TO 150 
RMU(I-dX.M) = AMU(J) 
NA(J) = BLANK 
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150 


IF(NB( J1 .NC.NMA) 60 TO 200 
s BMU(J> 

NB(J> = BLANK 
200 CONTINUE 

C 

STPR (I »=PMU(I+IX.M)-RMU(I+IX,N) 

C 

RSUH = RHUd'l'IXtM) « HHd) * RSUH 
PSUH = PMUd + IXiM) « UHd) ♦ PSUM 
C 

300 CONTINUE 

1CR = 0 

ICPsO 

DO 310 J=l«5 
STR( J)=0. 

STPtJ)=0. 

NR( J)sBLANK 
310 NP(J)sbLANK 

00 340 t s 1, IS - 

IF (STPRd)-O.) 320i340t330 
320 ICRsICR*! 

NRdCR)sNAHAd) 

STRCICR)=-STPR(I ) 

60 TO 340 

330 ICPsICP^l — - -- - - - 

NPdCP)sNAHAd) 

STP(ICP)=STPRd) 

340 CONTINUE 
NUHNA=ICR 
NUMNBsICP 

IDENTsNGAS - 

IF(H.GT.HGAS) lOENTsNSUR 

URITE(KOUT<20420) H«IDENT« ( NASYHI J) . STR IJ ) • NR ( J ) tJsltNUHNA) 

Si (NBSTHt J) iSTP( J) iNP( J)« J sliNUHNS) 

URITE(KOUTi2042S) FKF(H) iEXK(M) iEAK(H) 

EAK(H)=1000.*EAK(H) 

c 

IF(ABS(RSUH-PSUH> .LEi 1.E.04 * RSUN) GO TO 400 
URITE(KOUTi20300> H 
IFLAG =1 
STOP 
C 

400 CONTINUE 

C END DO 
C 

20100 F0RHAT(2I3) 

20110 FORHATdHli 30Xi 25H* KINETIC REACTION DATA • 

Si// 15Xi 30HTOTAL NUMBER OF REACTIONS i 13 

Si/ 15Xi 30H NUMBER OF GAS PHASE REACTIONSi 13 i/) 

20120 FORMAT(7XilHMi4Xi4HTYPEiT28iaHREACTIONiT60il4HPRE EXP FACTOR 
S T96i 8HTEMP EXPiTllOi lOHACTlVATION 

Si/ T81i IIH(MOLE-CM-S) iTllOi 1 IH ( KC AL/MOLE » /) 

20300 F0RHAT(66H SUM OF MOLECULAR WEIGHTS INDICATES IMBALANCE IN KINETIC 
S EQUATION iI3l 

20400 FORHAT(10(A4ilXl I ' 3E10.4) 

20410 FOKMAT( SFb.Oi T41i 5F5.0 ) 

20420 F0RMAT(5XiI3i4XiA4i6(A3iF3.1i1XiA4) ) 

20425 FORMAT(lH+i T82i E10.4i T96i F8.3i TllOi F10.4 » 

C 

RETURN 

END 
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o o 


maini 

INCLUOE BASIC 

include: heat 

INCLUDE HEATl 
INCLUDE ENERGY 
INCLUDE PROF " - 

INCLUDE PBOCl 
INCLUDE PROC2 
INCLUDE PROC3 
REAL NTUB 

COMMON/BLKl/ FLI JSI t JSI I t OKI ( JSI I t QR2 ( JSl ) 

C0H.10N/BLK2/TTI JSl) . AKT tOR •OT«EMIVT» AHH«CPPHt TTI *HN«HO 
C0MM0N/BLK6/PS( JS. JP2) *PPH( JSH «Mi <2 
COI".HON/BLK7/IKIN( JSl) 

COBMON/BLKe/VLtPTtWTtBtOCtCC.GRT 
C0MR0N/BLK9/TKltTK2tGKTltGKT2.SG,EK«Cl*C2iC3 
COHMON/Bk 1/HNN( JSl) tICOEF 

cohhdn/bk**/hho(Jsi) 

dimension AMF( JS) ,TG( JS) lOOl JSl) 
dimension ALPSV*UY)i PPSW(JY), PLPSV(JY) 
dimension TITLEC20) 


TAU=1.356E-12 

KR=5 

120 READEStIGOTiTLE 

WRITE (6,470) TITLE 

REA0(5,500) IS,NL,NIT,IL0SS,IGE0H,IC0EF,KR7 

IF(NIT .EO« 0) NIT = 20 
KR7S = KR7 

URITE(6,520) ISt NLt NIT. ILOSS, Il>EOM, ICOEF, KR7 
520 FORHAT(/» T10.23H* INTEGRAL PARAMETERS • 


1, //, T15.40HIS ( number OF SPt-CIES ) =• 13 

S. /, T15.40HNL ( GRID POINTS ) =t 13 

S, /, T15.40HNIT ( NUMBER UF OVERALL ITERATIONS ) 13 

S, /, T15.40HILOSS ( heat LOSS OPTION ) =. 13 

S, /, T15.40HIGEOM t PROBLEM GEOMETRY OPTION ) =« 13 

S, /, T15.40H1COEFF I HEAT TRANSFER INPUT OPTION ) =• 13 

S, /, T15.40HKR7 ( PRINT OPTION ) »« 13 ) 


READ 15.660) OTUB, VF. OIAMS. AK, LMIV 
AS=(l,-VF)*3.14l6*DIAMS*DIAMS/4. 

NTUU = VF • DIAMS • DIAHS / DTUB / DTUB 
IF (I6E0M .EO. 0 .and. NTUB .GT. X) ILOSS s 0 
IFdGEOM ,NE. 3) GO TO 34 ““ 
AS=(1.-VF)»D1AMS 

NTUB = VF • DIAHS / ( 3. l416*DTUB*DTUB/4. ) 
34 CONTINUE 

READ (5,660) AMt TI< P 
READ 15.660) TRES 
C=3.1416*0TUB*NTUB 
AV=4./DTUB 
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READ (5*730» C S « I ) « 1=1 «NU I 
HEAO<b.700MTL(n «I=2«NU 
READ«5.e30MIKlNCl) «I=2*NLI 

IF (ILOSS .EO. 0 .AND. 16E0H ,EU. 0>G0 TO 32 

IFdCOEF.EO.noO TO 20 

REAU<5.30)VL«AHH*TTltPT»WT 

REA0C5.30)AKT«0RfDT.0CtCC 

REAO<5.30)TKl,TK2.6KTltGKT2 

REAO(5«30)SG,EK«C1*C2«C3 

OY=DR 

CONTINUE 

FORHAT(8F10.5) 

IFdCOEF.Nt.DGO TO 32 

REAO(5«30MHAMB(I ) •Is2«NL) 

REAO(5«SO)TTI.AKT«OR 

Dt=OR 

continue 


WRITE (6*680 > 

680 FORHAT( ///* T10«26H* NONINTEGRAL PARAMETERS * 

$• //, T45*20H««* reactor DATA ••• > 

IFdGEOM ,NE. 1» MRITE(6«681) NTUB* VF. OIAMS. DTUB 

681 FORMAT(//« T15,38HH0N0LITH HONETCOMB CATALYTIC COMBUSTOR 


St 

/• 

T15.26HNUMBEH OF TUBES 

8« 

P1.2 

St 

"7% 

T15.26HV0ID FRACTION 


F7t5 

St 

/% 

T15.26H0IAMETEK OF COMBUSTOR 

8f 

F7t2 

St 

S) 

/% 

T15.26HEFFECTIVE CELL DIAMETER 

= t 

F7t2 


SH (CM) 
5H (CMI 


IFdLOSS.EO.O.OK* I6E0M. £0.0)60 TO 40 

WRITE(6.690) _ , 

3 FORMAT(//* T13.52HF0R CYLINDER COOLED EXTERNALLY OR FLAT PLATE PRO 
SBLEH ) 

URITE(6*692I OR 

> FORMAT( /* T15,26HTHICKNESS OF CYLINDER OR 

S« /, T15.26HTH1CKNESS OF PLATE =.F5,3. 5H (CM) ) 

3 continue 

WRITE (6*694 )TI«P» AH. AK.EMIV 

» FORMAT(//. T15,26HINLET TEMPERATURE =.F7.2,4H (K) 

S. /* T15.26HPRESSUKE =*F7.2* 6H (ATM) 

S« /, T1S«26HHASS FLOW RATE =*E9,3* 6H (G/S) 

S, /, T15.26HC0N0UCTIVITY =*E9,3*13H (CAL/S-CM-K) 

S. /* T15.26HEMISSIVITY a* F7.2 ) 


=.F7.2,4H (KI 

S« /* T15.26HPRESSUKE =*F7.2* 6H (ATM) 

S« /, T1S«26HHASS FLOW RATE =*E9,3* 6H (G/S) 

S, /, T15.26HC0N0UCTIVITY =*E9,3*13H (CAL/S-CM-K) 

S. /* T15.26HEMISSIVITY a* F7.2 ) 

DUHl = ABS(TRES(in 
DUH2 = ABS(TRES(2)) 

WRITE(6.670) DUMl* DUH2 
WRITE(6*725) (Sd)*I = l*NL) 

WRITE(6*715) (TL(I)*I=2*NL) 

5 FORHAT(//*T10*dNITIAL GUESSED WALL TEMPERATURES (TL)» 

S, /, 110*10610. 4*/*T11*10G10.4*/*T12*10G10.4*/*T13.10G10, 4 ) 

3 FORNAT(//»T10* *6010 POINTS CHOSEN (S)* 

S, /,T1()*10G10.4,/*T11*10G10.4*/*T12*10G10.4,/*T13*10G10.4 ) 

IFdGEOM. EQ.O. AND. ILOSS. EO.0)GO TO 45 
WRITE(6*20b00) TTI* PT * AHH* V(-* WT * DT * DR.AKT 
IF (I6E0M .EO. 3) WRITE (6 .20505 ) DC* CC 
3 FORHATdHl*//. T45t26H*** COOLANT PROPERTIES 
$, //, T15.19HC00LANT 

S, /, T15.19HINLET TEMPERATURE =* F7.2. 4H (K) 


a* F7.2 
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s • / • 

4 « / % 

S. /« 

s * / « 

s, //, 

s « / • 

$. /« 

* • / • 


S) 

20505 FORMATC 

$ //♦ 
$. 5H (CHI 


Tlb.l'JHPRrSSURE =• F7.2. 6H 

T15,1'*HMASS FLOW RAIL =t E9.3. 6H 
T15tX9HAPPR0ACM VELWCITY =. F7.2* 7H 
n5,19HH0LECULAR WEIGHT =. F7.2 
T10.21H FOR COOLANT TUBE: 
T15.19HTUBE DIAHETEH =. F7.2t 5H 

T15.19HTUBE THICKNESS s, F7»S. 5H 
Tl5, I9HC0NDUCTIVIT a, E9.3tl3H 


(ATHI 

(G/S) 

(CH/SI 


(CH) 

(CN) 

(CAL/S-CH-KI 


T15.41HCENTER TO CENTER DISTANCE (TUBE BANK! a« 


F7.2 


%, /, T15«41hDISTANCE between HONOLITH and tube bank a* F7.3 

4< SH (CHI I 

WRITE(6,20G00I TKl. GKT1« TK2« 6KT2« C1« C2< C3« S6« EK 
20600 FORMAT!//. T10*35H * THERHOCHEHICAL DATA FOR COOLANT 

S, //.TI5.37HTEMPERATURE THERHAL CONDUCTIVITY 

4. 2(/.T17.F7,2. T37.EI0.3I 

$. //.T15.65HHEAT CAPACITY CPH a Cl ♦ C2 * T ♦ CS / T**2 

S(CAL/G-HOLE - K» 

i. //.T20. 4HC1 a. E10.3.17H (CAL/6-HOLE - Kl 

$, /.T20. 4HC2 a. E10.3.20H (CAL/G-HOLE - K**2l 

5. /.T20. 4HC3 a. E10.3.15H (CAL-K/G-HOLE I 

S. //.T15.20HVISCOSITY PARAHETtRS 

$. /.T20. 7HSIGHA a, f7,3* 4H <AI 

». /.T20. 7HE/K a, F7.3. 4H (K| » 

45 CONTINUE 
C 


CALL READIN 


C 

CALL OBTAIN 
C 

CALL KININ 
C 

DO 190 lal.NL 
190 X(IlaS(II/DTUB 

IF (IGEOH - 21 196.197. 198" 

196 IF (IGEOH .EO. 01 GO TO 198 
C a 1.0 

RU a AH 
AS a DY 

DO 192 I a 1, JSI 
DO 191 J a 1. JSI 
1«1 Fd.JI a 0.0 
192 F (1.11 a 1.0 
GO TO 195 

197 AREA a 3.1416 / 4.0 • (DTUB»OTUH - DT ♦ DT I 
RU r AM / AREA 

PHD a HU • (UTIIB - DTI 
CALI. VIEWKUTUR.DTI 
GO TO 195 

198 ARLA a 3.1416 / 4.0 • DTUB ♦ DTUB 
RU = AM / AREA / NTUB 

RUD a RU ♦ OTUB 
CALL VIEW(DTUb.NL.S.FI 
195 CONTINUE 


t 
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nll=nl-i 

nlp=nl*i 

BHl = TAU*TRESm**4 

IF (TRESdj.GT.O.) GO TO 210 

BR1=0. 

TRESl = TRES(I) 

TRES(l) = - TRES(l) 

DO 200 J=2«NL 
200 F(J.2)=F(Jt2)+F( Jtll 

210 BR2=TAU*TRES(2)**4 

IF (TRES(2» .GT.O. » GO TO 230 

BR2=0. 

TRES2 = TRES(2) 

TRES(2> = - TRES(2) 

DO 220 J=2tNL 

220 F(JtNL)=FtJ.NL)+F(JtNLP) 

230 AEHN=AS*EHIV*TAU/AM 
EHIVV=EMIV 

— akk=ak ■ 

DO 240 LsStNL 

240 DS(L)=(S(L)-S(L-2) 1/2. 

DS(2) = DS«3) 

DS(NL>1I = DS(NU 

C 

IFdLOSS .EQ. II CALL BEES^ 

C 

IFdLOSS, EQ.O. AND. IGEOH.EQ.O)GO TO 245 
DO 241 L=1,NLP 

241 TT(LI = TTI 
245 CONTINUE 

TLd)=TL(2) 

DO 255 L=2.NL 
TW(L)= TL(L) 

255 CONTINUE 

TWd) = TW (2) 

TW(NL+1) = TW(NL) 

260 PLN=AL0G(P) 

DO 270 J=liISS 
PP( J)=ALPF( J»*P 
270 PLP( JI=ALOG(PP( J) I 
PM( 1 1=0. 

DO 260 1 = 1. IS 

260 PH<l)=PH(l)+ALPF(I)*WHd) 
PM(l)=P*PHd) 

PH(2)=PHd ) 

TG(l) = TI 

Td)=TI 

T(2»=TL(1) 

L = 1 

CALL FLAME 

IFdLOSS ,Ea. n CALL TRIDM 
TSV=TI 

DO 290 1=1. ISS 

ALPHI d )=PP(I )/PM(l ) 

PPSVd )=PPd) 

PLPSVd )=PLP( I ) 

290 ALPSU d »=ALPHl (I ) 
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n n 


PMSV1=PM(1) 

PHSV2=PM(2) 

HSV=SUMH/PM<1» 

H=HSV 

UMM=PM(1)/P 

RH=PM(l)/(e2,056*T(lH 

itdrso 

300 DLT>f = 0. 

AEHlsAEMN 
AEM2=-AEK1 
EMIu=EMI VV 
AK=AKK 

TF (TRFSl .LT. 0.» AEHlsO. 

IF (TRES2 .LT. 0.) AEM2=0. 

HsHSV 
T(1 l=TSV 
T(21=TL(2) 

PM(1)=PMSV1 

PH(?)=PHSV2 

CLOSSsO. 

Do 310 isltlSS 

ALPHHI)*ALPSVm 
PLP(1 )sPLPSV{l) 

310 PP(nsPPSVd) 

STEP=1. '■ ■ 

IB=1 

DO 420 L:2tNL 
320 KR7=KR7S 

IF( ILOSS.EQ.O. and. IGEO n.EQ. 0)60 TO 41 
TT(L»=TT(L-1) 

41 Continue 

IF (ITOR.GT.O) T(1)=TG(L) 

T11=T (1 ) 

lF(IGEOM.Ea.l)RUD=RU*S(U 
IF(IGE0M.E0.1)T(1)=(T(II+T(2) )/2, 

CALL PROP 
CALL COEF 
T(1 )=T11 
FIH=IB-1. 

DUP=1.-FIB*STEP 
AHH=(S(Ll-S(L-l ) )/AM*DUM 
ADSM = AHH * ah / RU 

CHH=AS*AK/AH*DUH " 

CH=AHH*C 

A1=A1*C»AHH 

A2=A2*C*AHH 

0RADIN=-EMIV*(F(L«1 )*8R1+F(L«NLP)*BR2) 
DO 330 1=2. NL 

IF (I.EO.L) GO TO 330 - 

B = TAU*TL( n**4 
ORADIN=QRADIN-F«L.I )*B*EHIV 
330 CONTINUE 

CALL flame 

IF (ICON.EO.O) GO TO 350 
GO TO 450 ■ “ ' 
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350 IF (IB.CO.l) GO TO 360 

IF (H0D(IB*2)«C^0.1) GO TO 320 

IB=IB/2 

STEPsSTEP*2. 

60 TO 320 

360 TDIF=ABS(1.-T«2»/THU) 

IF a01F.GT.DLTX> DLTX=TDIF 
IF (L.GT.2) GO TO 3S0 
PMSV1=PB«1> 

PHSV2=PH«2) 

TSV=TU> 

00 370 1=1. ISS 

PLPSV(I>=PLP(1> 

370 PPSV(I)=PP|I> 

380 CONTINUE 

TG<L>=T(1> 

TL«L>sT«2) 

IF(IUOSS.EO.O.AN0.1GEO«.EO.O>GO TO H2 
00 «L > s ( TT (L ) -TT < L-1 n •AMH*CPPH 

IFdGEOII .EO. 1 .OR. (IGEOM .EO. 0 .ANO. ILOSS .EQ. 1) 
% OOIL)=HNN(L)*<TJ2I-TTI )*>S«L)-Sa-l>) 

42 continue 

XSUHsO. _ 

00 390 1=2. f ■ 

390 XSUM=XS0B*XNUBXP<1> 

XNUBH=XSUM/XPP 
UM1=PH(1I/P 
WH2=PH(2>/P 
00 too J=1«ISS 
400 AMF(J|=PP(J)/P 

RH=PfMl)/C62.056*T«l) I 
RH2=Pf1(2)/«82.056*T(2> I 
410 CONTINUE 

IF (ICON.GT.O) GO TO tSO 

420 CONTINUE 

c 

IFIOLTX.LT. .001) GO TO 52 
IF<1T0R.GE.N1T)G0 TO 52 
GO TO 425 
52 CONTINUE 

IFdGEOB .EQ. 3 .OR. IlGEOM .EQ. 0 .AND. ILOSS. EQ. 0) 

QTOTAL = 0.0 

UO 51 JJ = 2. ML 
51 QTOTAL = QTOTAL ♦ OQIJJ) 

WRITEI6. 20700) TTI. TT(NL). QTOTAL 

421 CONTINUE 

20700 FORflAT ClHl .//.T15.36H** OVERALL HEAT TRANSFER RESULTS 
S. /.T20.26HC0OLANT TEMPERATURE IN =, F7.2 

S, /.T20.26HCOOLANT TEMPERATURE OUT =. FT. 2 

$, /.T20.26HTOTAL HEA» TRANSFERED TO 

S, /.T20.26HCOOLANT =. E10.3 

S) 

C 

CALL STHPRT(TG.QQ) 

IFdGE0M.E0.3)TZ=TTI 

IFCIGEOH .EQ. 3) CALL 0PT2 • 01 AMS.N1 UB . TG . AS.OTUB ) 


GO TO 421 


«« 

4H <K) 

4H «K) 

7H(CAL/S) 
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IF(KR7.LT.3)G0 to 425 
IFJ IGE0«.E0.3)G0 TO 425 

IF( IbEOH.EO.O.AND.ILOSS.EO.O>60 TO 425 
WKITE(6«54) (HNN(L) •L=2«NU 
WHITE (6 >55) (HHO(L) •L=2«NL) 

54 F0KHAT(/.5X.52HHEAT TRANSFEK COEFFICIENT COOLANT SIDE (CAL/S-K-CM2 
$) ,/,5x.8»2x.E10.3,2X) ./) 

55 FORMAT*/, 5X.52HHEAT TRANSFER COEFFICIENT FUELAIR SIDE (CAL/S-K-CN2 
S) ./.5Xf ei2X,E10.3,2X) ,/) 

425 CONTINUE 

IF(DLTX.LT..001) GO TO 110 
IFtlTDR.GE.NITjGO TO 110 

ITDR = ITDR + 1 — 

GO TO 300 

110 WR1TE(6,610)ITOR,DLTX 
450 STOP 

460 FORMAT (20A4» 

470 FORMAT<1H1,//,T15«20A4,/ ) 

500 FORMAT (1213) 

510 FORMAT (8F10.0) 

530 FORMAT tlOX,6HALPF , 8E10 .4/21X , 8E10 .4 ) 

540 FORMAT (E10.5) 

550 format (1H,I4,7E14.5,2A4) 

560 FORMAT ( 6E9 .6 ,6X ,F6 , 0 ,6X , 2A4 ) 

570 FORMAT tl3l 

580 format (10X.21HKINETC REACTION OaT A//15X , 20HNUMBER OF REACTIONS=I3 
* /) 

590 format (20F4.1) 

600 format (10X,20F4.1) 

610 FORMAT ( 5 ( A4 , IX ) , 5X , 3E10 .4 ) 

620 format (5X, I3,3X, A4»1H+*A4»1h+, A4,3H--=,4X, A4,1h+* A4,6X,E10,4,6X,F 
S 8.3,6x,F10.4) 

630 FORMAT (8E10.3» 

640 FORMAT (lOX.BHDlFF FAC.10E10.5) 

650 FORMAT ( 1 5, 5X , 7FX0 . 0 ) 

660 FORMAT (8F10.0) 

670 FORMAT*//, TIO, RESERVOIR TEMPERATURES 
4. //,T15, 'UPSTREAM =',G12.5 

4. /,T15, ‘DOWNSTREAM =• ,G12. 5 ) 

700 FORMAT *8E10,5J 

710 FORMAT *10X,10F10.4/15X,10F10,4/20X,10F10,4/25X,10F10.4//I 
7?n format *10X,10F10.4/15X,10F10,4/20X,10F10.4/25x,10F10.4//) 

730 FORMAT *8E10.5» 

740 FORMAT * lOX , 9HENTHALPY=, ElO . 5 ,8X ,6HM0LWT=, ElO . 5 ,8X , 4HRH0=, ElO ,5 
4,/ 4HALP=,10E10.5) 

750 FORMAT *1X,'XNUB = • ,E10 ,5 , lOX • • XPP = ‘,E1U,5) 

760 format *1X,‘MEAN NUSSELT number X •,E10«5) 

770 FORMAT *1H0,34HBULK AND WALL PROPERTIES ALONG BED) 

780 FORMAT * 1 OX , 9HOIST ANCE= , ElO . 5 , 5X , 5HNODE = , I4 , 5X , 5HDLTXX , El 0 . 5 ) 

790 FORMAT *10X,6HBULK ,3HTx ,El0.5,4HMUx ,E10.5,3HHx ,E10.5,5HRH0= , 
4 E10.5/10E10.5) 

800 FORMAT *10X.6HWALL ,3HTx ,El0.5,4HMUx ,E10.5,4HHWx , ElO . 5 ,5HRH0« 

4, E1U.5/10E10.5) 

810 FORMAT*///. T15,25(1H*) 

5, //,T15, ‘ITERATIONS =‘,I5, • TEMP ERROR x», G13.5 

4, //,T15,25*1H*) ) 

830 FORMAT(40I2) 

END 
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SUHROUTIME OBTAIN 
INCLUDE BASIC 
INCLUDE PROF 
INCLUDE PROCI 

DIMENSION LNGTHI2) 

DATA LNGTH/ 2t 30 / 

DATA IENO/3HEND/ 



OK s 0 

DO 9 IsliNWANT 

IWANT(1.I)=0 

9 CONTINUE 

C CURVE FIT DATA ASSIGNED TO UNIT 11 
C TABULAR DATA ASSIGNED TO UNIT 12 

ICC = -I 
INKAS =12 

IF(ICC.LT.O) 1NMAS=11 

rewind INMAS 

c originally COMPAR 

NLNGTH = LNGTH(2> 

1F«ICC,LT.0) NLNGTHsLNGTH(I) 

HATCH = 0 
lEOF = 0 

UR1TE(6* 300) 

300 FORMAT(1h1*30X«24h* ThERHOCHEMISTRY DATA • «/) 

IF(lNMAS.EQ.ll) WRITE (6.10) 

10 FORHAT(1H0.15X,59HCURVE FIT OF DATA IN FORM CP=RB+RC*T+RD/ ( T*T ) (C 

»AL/«M0LE*K)//TX*1?H HF(CAL/M0L) «2X« 10HH(CAL/MOL) . 5X . 2HRB * IIX t 2HRC 
».11X.2HRD.7X.12HS(CAL/M0L/K).6X«9HTU(K).4X.4HNAME /) 

IF(1NMAS.E0.12| WRITE (6. ID 

11 FORMAT! lH0tlbX«12HTABULAR OATA/12Xi 4HNAHE* 8X* 2HMW«6X« IIHHFCCAL/ 

♦M0LI/15X. 76HCP(CAL/(M0L*K) ),H(CAL/MOL) .FlCAL/(MOL*K) GIVEN EVER 
•t lOOK FROM ion TO 5000K ) 

C ORIGINALLY GtTSPt(IEOF) 

HEAD(INMAS.4) ihdg 
I* FORMAT I A2! 

3 READ! INMaS«1) ISPECI «XHW • XHRF 
1 FORMAT!10X.2A4.2X,2F10.0) 

IF llSPECI!!) .EO. lEND) lEOF = 1 

IF !1E0F .NE. 0) GO TO 200 

150 DO 160 I=1*NWANT _ 

IF llWANTIltl) ,NE.~0) GO TO 160 
DO 120 J=li2 

lF!IWANTlJ+ltl).NE.lSPLCI(J)) GO TO 160 
120 CONTINUE 

MATCH = MATCH 1 
IWAMT!1.D = match 
CALL GETDAT!!) 

IF IMATCH .EQ. NWANTt GO TO 200 
GO TO 3 
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160 continue 

C ORIGINALLY GETOUM(NLNGTh) 

00 fl I=1.NLNGTH 

READ (INMAStSOlOOl lOUN 
20100 FORMAT (Al) 

a continue 

GO TO 3 

C ORIGINALLY CHECKtOKi 

200 DO 5 I=1«NWANT 

IFdUANTdtD.NE.O) GO TO 5 

WRITE(6t6MlWANTIJ«I) tUsStS) 

6 FORMATdHOt 39H*** NO DATA IN SPECIES MASTER FILE FOR « 2A4» 
S 4H 

OK s OK ♦ 1 
5 CONTINUE 


if(ok.eQ>oi go to 7 
write (6.201201 

■20120 FORMAT (IHO, 47H*** PROGRAM STOPPED IN OBTAIN - UNKNOWN SPECIES ' “ 
S. 14H REOUESTEO •»* 

%) 

STOP 

7 RETURN 
END 
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SUBROUTINE 0PT2 | 01 AMS tNTUBt TG t AS lUTUB I 

REAL NTt NTUB 

INCLUDE BASIC 

INCLUDE HEAT 

INCLUDE ENERGY 

INCLUDE PROCl 

INCLUDE PR0C3 

C0MH0N/BLK6/PS( JS,JP2) iPPH(JSl ) iHtT^ 

C0BM0N/BLK2/TT(JSl)tAKTiDRiDT.tMlVTiAMH.CPPH.TTI.HNtH0 

C0MB0N/BLK8/VLfPT,WTiB*DCtCCiGKT 

DIMENSION TG( JSl) 

DIMENSION TII (2) •AHM(2) tHH(2) «CPH( 2 ) t HO ( 2 ) • C ( 2 ) 

DIMENSION U<2) ,E(2).SUH<2) .*A(2I iXb(2),XC«2)tIAVG(2> tTTWA<2) 
DIMENSION BB(2) «TTFC2> «0C(2»tCPEWW(2) tORR«2) 

W = 1. 

AREA=DIAMS*W 

BsOIAMS 

CT=3.1<416*DT 

Y = CC / B 

2=DC/DT 

NT = 1, / DC 

TAV=200. _ 

TFsSOO. 

TII(1)=TT1 

TII(2)=TTI 

FF = SQRT(1.+Y*Y> - Y 
1F(2.GT.2, »G0 TO 2 
FFl=1.01+0.049*2“0.059*2**2 
GO TO 3 

2 FF1s1,26-0.22H*2*0.0145'»2**2 

3 FF=FF»FFI 

FP=(FF*EMIV*TAU*NTUB)/(NT*B*CT) 

NITsl 

DO 99 K=l«2 _ _ 

GHAX=AM/B*DC/CC 

T2=TII(K) 

3*f CALL CPR0P<T2> 

PM(I)=PPMU) 

T(1)=TI 

IF(K.E0.2)TI1 l=TG(NL) 

IF(K.E0.2)PM(1)=PPM(NL) 

TAVG(K)=T(1) 

32 J = 2 

IF(TAVG(K).LT.TU(M.l) ) J=1 
CPLWW (K )=0.U 
DO 6 M=1.IS 

CPP=RC(M, J)4TAVG(K)*RD«H, J)*RE(M, J)/TAVG(K)**2 
NI = NL 

IF(K.E0.1)NIr2 

G CPEWW (K )=CPEWW<K )+PS(NI «M)*CPP 
CPGd )=CPEWW(K)*PM(1I 
CALL PROP 
REy=riT*GMAX/VISC 
U (RLY.GT.3U0. )G0 TO 30 
Ah = U 
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AX = t),5 
AC = 1.5 
GO TO 31 

30 AR=0.3 
AXsO.4 
AC=0.65 

IF(REY.LT.3000. )AC=1.0 

31 HOCK )=AR/REY*»AX*CPEWW(K)*GMAX/PR**0.667*AC 
AhH(K)=AMH 

HH(K)=HN 

CPH(K)=CPPH 

U(K )=!./( I./HHCK) +DR/AKT^l,/HO(K ) ) 

E(K)=U(K)*CT/(AHM(K)*CPH(K1 > 

C**» CALCULATE RADIATION TERM 
SUM(K)=0.0 
N=1 

IF(K.EQ>2)N=NL'fl 
DO 1 Is2<NL 

1 SUM(K)=SUM(K) + (S(I)-StI-l) > *F ( I *N) *TL ( 1 1 **4 
IF(K.EQ.1)N=3 

SUH(K)=(NTUB*3,1416»DTUB*SUM(K»+AS*TL(N-1)**4»*FP/NTUB 
C(K)=CT/{AHH(K)*CPH(KJ ) 

XA(K) = (B-UtK)*C(K)/E(K»*(B“«l.-EXP(-E(K)*B) )/E(K» > ) 

XB(KI=(B-XA(K ) )*SUH(K)/U(K)+TII (K)*tl.«EXP(-E(K)*B) l/ECK) 
XC(K)=2.*AM*CPEWW(K)/tCT*NT*U(K)) 

IF(K.E0»2>60 TO 4 

TAVG(K) = (XB(K)-XCCK)*T1 l/tXA(K)“XC(K) ) 

TGI=TAVG(K)*2.-TI 

TTWA«K>=TAVG(K»“AM*CPEWH(K>*(TGI-T1)/(H0(KI*B*CT*NT) 

GO TO 20 

4 TAVG{K)={XB{K)*XCCK)*TBCNL) »/<XA(K)*XC(K) > 

T6NLP=TAVG(K)*2.-TG(NL) 

TTWA(K»=TAVG(k)-AH*CPEWW(k>*(TG(NLI-TGNLP)/(HO(K)*B*CT*NT) 

20 IF(ABS(TAVG(K)-TAV).LT,1,0)GO to 55 
TAV=TAVG(K» 

NIT=NIT+1 

GO TO 32 ■' 

55 BR(K)=(SUH(K)*U(K)/H0(K)+U(K)*TAVG(K) » *CT/ ( AHM ( K J *CPH ( K ) 1 
C •* CALC, working FLUID EXIT TEMP AND HEAT REMOVAL RATE 

C •• 

TTF(K> = (BB(K)*(1.-EXP(-E(K)«B) »/E(K)+TII(KI*EXP(-E(K)*B) ) 
TII(2)=TTF(1) 

IF(ABS«TTF(K)-TF) .LT. 1.0)60 TO 33 
TFrTTFcK) 

T2=(TTF(K)+T1I (Kl )/ 2 . 

GO TO 34 

33 QHR(K)=SUM(K)*CT»B 

0C(K)=AHHCK)*CPH(K)*(TTF(K)-TII (K) ) 

99 CONTINUE 
C 

OTOT = CTTFC 1 )-TII ( 1 ) )»AHH(1 )»CHH( 1 )♦( TTF(2)-TII (2) ) * AHM ( 2 ) *CPH C 2 ) 
WRITECG. 20600) TII(l). TII(2). TTF ( 1 ) » TTFC2). QCCl). OCI 2 > 

TGI. TG(NL). TI. TGNLP 
IF(KH7.LT.3)G0 TU 50 
WRITE (b.51 )HH( 1 ) .HH(2) .HOCl ) ,H0(2) 

51 F0IIMAT(//5X.7HHH(1 )= . El 0 . 3 . 3X . 7HHH ( 2 ) = . El 0 . 3 . 3X . 7HHO ( 1 ) = »E10.3. 
$3X,7HMO(2)= ./) 
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so CONTINUE 

20800 FORMAT! IHXt//. no. 3«H** FOR HEAT EXCHANGER PROBLEM OPTION 2 
S. //.T*«6.2bHEXCHANGER 1 EXCHANGER 2 

%, /.T15.30HCODLANT TEMPERATURE IN (K> s, T47 . 2 ! FT ,2 ,8X ) 

S. /.Tlb.SOHCOOLANT TEMPERATURE OUT (Kl =.T47 .2 ( F7.2 .6X ) 

S, /.T15.30HTOTAL HEAT TRANSFERED TO 

S« /,T15,30HCOOLANT (LAL/S> = ,T47 .2 ( E1».3 .6X ) 

S. /.Tlb.30HTEMPEHATURE OF FUEL/AIR (K)l 

$. /.T1&.30HUPSTREAH OF HEAT EXCHANGER s, T47 .2 ( F7.2 , 0X ) 

$. /.T1S.30H00WNSTHEAM OF HEAT EXCHANGER s,T47 . 2 ( F7.2 , 8X ) 

S) 

RETURN 


END 
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PROCS 

BASIC 


PKOC 

parameter 

PARAMETER 

parameter 

PARAMETER 

PARAMETER 

PARAMETER 

PARAMETER 


JP = 12 
JP2 = 2 
JD = JP 
JS = 40 
JSl = JS ♦ 1 
JY = 2 * JP 
JX = 5 


JP 

1 


+ 4 


end 

HEAT 


PROC 

REAL HAMB 

COMMON/HEAT /CHtCHMtCHHf CM 

COMMON/TRANS/ A 1 t A2 « A3 » A4 • DBTW t ChSTG *0 IN ( JS ) 
COMMON/PROPS/X( JS) tPRfSCtTCONOf ViSCtXNUfSTtRUDtRU 
COMMON/CHEAT/XPtXTPf VFfXNUI 


/ F2 / 


COMMON 

Sf 

s « 

COMMON/ AMb 
COMMON/ TWS 
COMMON / BBB 

St 

St 

St 


AEMlt AEM2 
TRES(2) 

S(JSl)t AK, OY 
HANR(JSl)f TAMB0« 
TWCJSDt TWKJSlIf 
^ B1 (JS1)« B2 USDf 
B6 (JS1)« 87 (JSl)f 
BlllJSDf 0l2«JSl)f 
Bl6(JSl)t BP3(JS1)« 


0S(JSl)t NLt TZt TL(JSl) 


TAMBI 

TU2(JS1) 

B3 (JS1)« 
B8 (JS1)« 
ei3( JSl) t 
BP5(JSl)f 


B4 (jsn« 
B9 (JSl)f 
B14IJS1) « 
BP6(JS1) « 


B5 (JSl) 
B10( JSl) 
B15(JS1) 
BP7( JSl) 


END 

ENERGY PROC 

C0MM0N/ENER6Y/QC0NVt0C0N0fQCHEMf0RAD«0RA0IN«0L0SSt AMtHU 
COMMON/VFA/EMIVtTAUtF( JSlf JSl ) 

C0MM0N/SAVE/HEW(JP2) tSPLPi JP2) tPMSl «PMS2tTSl*TS2 


END 

HEATl 


PROC 

COMMON/HEATI/ITURB 
COMMON/MEAN/ 


XNUBXP(IOO) tXNUBtXPPtXPPS 


END 

PROF 


PROC 

COMMON /CARDl/ 

COMMON /lUNIT/ 
COMMON/LIN/ WM(JP) 
COMMON /WANTS/ I WANT ( 


ISPECI(2)f 

INMAS 


XMWt xhrf 


3t JP) f NWANT. XALPF( JP2) 


END 

PROCl 


PROC 

COMMON 

S f 
S I 
S t 

COMMON 

S f 


/ Fl / FT, Ht iCONt ISt ISPt ITER* 
Lt MTt MGASt Ni P# PLN 
ICCt NAMB(JP)t VNU(JP*JP) 
ILOSSt IbEOM 

/ F3 / ALPF( JP2 ) # BM JP) 

RT(2)t SUMP(2)t TSQ(2) 


KRt KR7 


END 

PR0C2 PROC 

COMMON 

S I 

COMMON 
:i> « 

S t 


END 

P'(0C3 


PROC 

COMMON 

St 

St 

St 

St 


/ F4 / EAK(jX)t EXK(JX)t FKF(JX)t PMU(JP2tJX) 
RMU( JP;>f JX) , PHK JP?» JX) t THMU( JDt JX) 

/ F5 / AHt AHHf Ct CHHt UlAMt FCt ISSt SNSM 
bUMH* Twr 
ADSM 


/ F6 / A(JYtJY)t ALPHI(JY)t CPF(JP2)t CPGI2) 
Hl(JP2)t H0SCJP2)t NAMA(JP)t PLP(JP2) 
PM<2)t PP(JP2)t KB(jP*2)t RC(JP*2) 
R0(jPt2)« RE(JPt2)t RF(JPf2)f RHS(JY) 
SB( JP2) * T(2) • TU( JP,2) 


end 
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SUBROUTINE PROP 
INCLUDE BASIC 
INCLUDE HEAT 
INCLUDE PROP 
INCLUDE PROCl 
INCLUDE PROC3 


PMU1=0. 

PMUa=0. 

PMU3S0, 

PHU5=0. " ' 

PMU6=0. 

CPTILsO. 

DO 1 I=1«IS 
PHU1=PMU1*PPJ n*BF( H 
VASPPCI »/BF«l » 

PMU2=PMU2*VA*W«m 

PMU3=PMU3*VA 

1 CPTIL=CPTIL*VA*CPF<I) 

DO 2 I=I<IS 
VA=PP(I)/BF(I) 

VB:VA«UH(I> 

VC=PP(I>*BFIIJ 

PMU5=PMU5+VB/ ( 1 ,385-0.029* VC/PMUl) 

2 CONTINUE 

PHU6=CPTIL/I.9e69*,29*Pt1U3 
TPMU=T ( 1 » ** ,659/PHUl 
VMU=2 , 096E-7*TPHU*PHU5 
VISCsVMU 

T C0ND=4 . X65E-7*TPMU*PHU6 
PR=V»U*CPG < 1 » /TCOND/Pm 1 > 

SC=PMU5/PMU2 

3 FORHATI///* molecular TRANSPORT PROPERT1ES»/10X, ‘TEMP ■•«E10.5«*K 
*.2X,*VISC =*,E10.5,* 6M/CM-S*,ZX,» ICONO s»,E10,5,* CAL/CM-S-K* «2X 
••PR =• ,E10,5,2x* ‘SC =*»E10,5) 

RETURN 

END 



u u 


SUBROUTINE PRTLIN(NAME .ARRAY, NVALS> 
DIMENSION ARRAYdl 


WRITE(6. 20000) name, NVALS. ( ARRAY ( I ) , 1=1 ,NVALS ) 
20000 F0RMAT(2X, Ab, 14, 3X, 10E11.4, / IbX, 10E11.4) 
RETURN 
END 
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SUBROUTINE READIN 

READS IN THE CARD(S) INDICATING THE SPECIES AND 

THE ALPF and ALPE VALUt-S THE PROGRAM HILL USE 
INCLUDE BASIC 
INCLUDE PROF 
INCLUDE PROCl 
INCLUDE PROC2 
DIMENSION ISPECC2J 

DATA IEND /3HEND/ 


UMAX = JP 
J = X 


IF(IS ,GT. UMAX) GO TO 200 
URlTE(6f 10) 

10 FORMATIlHltlbX.TOH* SPECIES INITIAL MOLE FRACTIONS.FIRST GUESSEStA 
•NO DIFFUSION FACTORS • t/> 

URITE <6,20090) 

0090 FORMAT)/, 15X, 6HSPEC IE , 16X , 4HALPF , 16X, 9HALPE, lOX 
•.16HDIFFUSI0N FACTOR/) 

DO 160 J = 1, IS 

READ (b,20100)ISPEC,KNY,XF,XE,XbP 


0100 FORMAT(2A9,I2,2E10.S,F10.9) 

URITE(6, 20120) ISPEC , XF , XE , XBP 
0120 fORMAT(16X,2A4,12X,E10,3,10X,E10.3,10X,F10,4) 

DO 140 1=1,2 

IUANT(I+1,J) = ISPEC(I) 

140 CONTINUE 

XALPF(J) = XF 
XALPF(J+IS) = XE 
BF(J)=XBP 
160 CONTINUE 


160 NUANT = IS 
ISP = IS + 1 
ISS = 2 • IS 
GO TO 220 


MOPE than UMAX SPECIES REQUESTED 


200 WRITE (6,20200) UMAX 
20200 FORMAT (43H0»** PROGRAM STOPPED IN READIN -- MORE THAN, Ib, 
$ 22H SPECIES REQUESTED ••• 

S) 

STOP 

220 RETURN 
END 
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SUBROUTINE RERAY(N.CtNNtOiNNNtLS.IS) 

C DIRECT INVERSION PROCEDURE -- C IS REPLACED BY C**-l 

INCLUDE BASIC 

DIMENSION D(JYtl) <SD( JY) tC(JYtl) «L( JY) <S( JYI iLL( JY) «LLL( JY> tLS(l) 

C 

K0UT=6 
N1=N+1 
NP=N+NN 
DO 11 1=1. NP 
LLL(Ii=I 

IF(LS(D) 113,113.112 

112 L(I)=LS(I) 

GO TO 11 

113 L(I)=I 

11 CONTINUE 
IX = “1 

IF(IS+2) 111,109,1X1 ■ 

106 FORMATIllH L ( I ) , 1=1 , 13, 5X , ( 3013 

107 FORHATdSH ( (C ( I , J > , J=1 « 13 , 12H I , (0 <J ) , J=1 , 13 , 6H) , Id , 13 , 15H ) BE 
IFORE RERAY) 

loe F0RMAT(2X, IIEIO . 3/ ( 12X , 10E10«3)| 

109 URITE(KOUT,107) NP.NNN.N 

- WRITE(KOUT,106)NP,ltTir,l=l,NP) ' 

IX = 0 

DO 110 1 = 1, N 

110 WR1TE(KOUT,108) (C(I, J) , J=1,NP) , (0 C I , J ) , J=1 , NNN) 

111 IS=-1 

C TRIANGULATE MATRIX 

DO 15 1=1, N - - - 

DO 160 H=1,NP 

160 S(M)=ABS(C<I,M) ) 

IF(IS) 10,16,16 

16 IS=0 

SD(1)=1. 

GO TO 12 

C REDUCE ROW I BY PRECEEOING ROWS 

16 DO 17 J=2,l 
K=L(J-1) 

DIV=-C(I.K) 

IF(DIV) 161,17,161 

161 C(I,K)=0. 

D0162 H=1,MP 
DIVC=OIV*C( J-1 ,M) 

S(M) = AMAXl(S(r) .ABS(DIVC) ) 

162 C(1,H)=C(1,H)+DIVC 
IF(NNN) 17,17,163 

163 PO 164 M=1,NNN 

164 U ( I ,M )=n ( I ,H ) +DIV*D ( J-1 ,M ) 

17 CONTINUE 

r SEEK MAXIMUM PIVOT 

12 DIV=0. 

K=L ( 1 ) 

J=I 

I'O 13 JJ=I,N 
h=L( JJ) 
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IF(ABS IC(I •«» )-DIV>tS«13tX2l 
121 HIV=ABS (C(ltM)) 

K=M 

J=JJ 

13 CONTINUE 

SDII)=C1V/S(K) _ 

L(J)=L(I) - - - 

L(I )=K 

IF (SOC I )“l»E"fl> 131tl31»14 
C SINGULAR MATRIX RETURN 

131 IS=-I 

WKITE(KOUTtl32) 1«K 
■ “ C(ItK)=l.E+30 

S(K)=1.0 
60 TO 12 

132 F0RMATI115Xt2l3) 

IH OIV=C(I,K) 

C(ItK)=1.0 

KsLLLCJ) 

LLL(J)=LLL(1) 

LLL(I )=K 
LL(K)=I 

C NORMALIZE ROW 

IF(NNN) l*»3,143tl41 

141 DO 142 JrltNNN 

142 D(I.J)=0(ItJ)/DIV 

143 UO lb J=1,NF 

15 C(ItJ)=C(ItJ)/DIV 
IFIIX) lS2tlb0tlb2 

151 F0RHAT(24H PIVOT ROW/COL/RES. RATIO 5 | I4«1H/13,1H/E9.2 « 1H« ) I 
150 URITE(K0UT«151) ( I tL < I ) tSD I 1 ) i 1=1 tNP ) 

c diagonalize matrix 

152 NM=N-1 

DO 20 I=1«NM 
K=LtI*l) 

00 20 jsltl _ _ 

DIW=-CIJtK> 

IF(01V)19»20«19 

19 C(J«K)=0. 

IF(NNN) 191.191tl92 

192 DO 193 MsltNNN 

193 Dl J.H)=D( JtM)+DIV*D( I+ltM) 

191 DO 201 h=ltNP 

201 r.( JtM>rC( J»M)+DIV*C«I + 1»M) 

20 CONTINUE 

c interchange columns 

DO 30 11=1. NP 

1 = II 

21 J=L(I) 

Ltn = l 

IF ( J-nP2.30.22 

22 1F(TS)25,23.25 

23 DO 24 n=l.N 
S«KJ=C(M.I ) 

24 C(M.I )=CtM.J) 

IS=I 

T “ I 
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GO TO 21 

2R IF(IS-J)26.2e.2fe 

26 DO 27 M=1*N 

27 C(M.I )=C(HtJ) 

I=J 

GO TO 21 
26 DO 29 H=1,N 

29 C(M.n=S(H» 

IS = 0 

30 continue 

C INTERCHANGE ROUS 

UO 40 IIslfN 

I=II ■ ■ " 

31 J=LL(I) 

LL(I)=I 

1F(J-I 132,40,32 

32 IF(IS)35,33,35 

33 DO 34 H=1,NP 
StM)=C<l*M) 

34 C(1.M)=C( J.M) 

IF(NNN) 343,343*341 

341 DO 342 M=1,NNN 
SD(H)sD<I.HI 

342 D(I,M)=D(J,M) 

343 IS=I ■ ■ ■ 

I = J 

GO TO 31 

35 IFaS-J)36.38,36 

36 DO 37 M=1.NP 

37 Cl I*M)=C( J.M) 

IF(NNN) 373,373*371 

371 DO 372 H=1*NNN 

372 D(:*M)=0(J*H) 

373 1=J 

GO TO 31 

38 DO 39 M=1*NP 

39 C(1*M)=S(H) ■ 

IF(NNN) 393*393*391 

391 DO 392 H=1*NNN 

392 D(I*M)=SDtH) 

393 IS=0 

40 continue 

IF(IX) 411*409*411 

407 F0RMATI15H I (C ( I * J ) * J=1 * 13 * 12H 1 , ( D ( J ) *d=l * 13 * 6H ) * I>1 * 13 * ISH ) AF 
ITER RERAY 1 

409 WRITE(KOUT*407) NP*NNN*N 
DO 410 1=1, N 

410 URITE (KOUT ,106)(C(I*J1*J=1*NP),(D(1*J) * J=1*NNN) 

411 RETURN ' 

END 
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SUBROUTINE STNPRT ( TG tOO ) 

integer pages 

INCLUDE BASIC 
INCLUDE HEAT 
INCLUDE PROCl 
INCLUDE PR0C2 
INCLUDE PHOC3 

C0HH0N/BLK2/TTI JSD t AKTiDRfDT tERI VT • AMHi CPPHi TTI * HN*HO 

C0MH0N/BLK6/PS( JS. JP2) tPPM< JSl) tW.TZ 
C0MM0N/BKS/PX( JS*JP2I 
C0HM0N/BK4/HH0(JSn »TTW< JSn 
DIMENSION TG(1)« 00(1> 

C-- — - 

ITEMS = 10 

PAGES = (NL ITEMS - 2) / ITEMS 
DO 300 I s li PAGES 

lEND = I • ITEMS ♦ 1 
IBGN = IEND - ITEMS 1 
TEND = AMINl (lEND.NL) 

WRITE(6«20100) 

WHITE (6 • 20120 > (S(J) i J=IB6N< lEND * 

WRITE(6«20140 I (T6(J) • J=Ib6Ni lEND) 

WRItE( 6. 20145) ( TL t J) » JsIBGNt lENO ) 

UKItE(S«20150) 

DO 200 K = 1. IS 

WRITE ( S«20160 )NAMA (K ) « (PXIJ.K) , JsIBGNt IEND) 

200 continue 

WRITE(6«20165) 

DO 260 K = ISP, ISS 

WRITE) 6*20160 )NAMA(K- IS) * (PX( J«K) «J= IBGN* lEND) 

250 continue 

IF(IGE0M.E0.3 .OR. (IGEOH.EO.O .AND. ILOSS.EO.O)) GO TO 300 
URITE(6, 20170) (TT<J) *J= IBGN .TEND) 

WRITE I 6 *20172) (TTW« Jl . J=IBGN , ILND ) 

WRITE (6 *20175) )QO(J) *J= IBGN* IEND I 
300 CONTINUE 

moo FORMAT ( IHl *// * T24,62HAXIAL DISTANCE ALONG MONOLITH * CYLINDER 
t OR PLATE S (CM) ) 

m?n FORMAT (// *T20 . 10(F7.3*3X) ) 

0140 F0RMAT(//,T4, 13H0ULK TEMP (K), T20* 10)1X*F6.1*3X) ) 

0145 FORMAT (// *T4 , 13HW all TEMP (K), T20* 1 0 { IX * F6 . 1 * 3X > ) 

0150 FORMAT)//, T4.13HM0LE FRACTION 

S //.T6,12HHULK SPECIES , / ) 

OlhO FORMAT)/, T11*A4,6X. 1 0 ( 2X , F7 . 5 * I X ) ) 

0165 FORMAT)//, Tb,lPHWALL SPECIES , / ) 

0170 FORMAT ( //,T4 , 7HC00LANT 

/,T4,10HTEMP (K), T20, 10 ( IX , F7.2 , 2X ) ) 

017? FORMAT)//, T4,12HC00LANT TUBE 

», /,T4, IOHTFMP (K)<T20, 10(1X,F7.2,2X) ) 

0175 FORMAT)//. T4,1?HHFAT LOSS TO 

/ , 1 M . 1 ..iirnoi f.f 1 (OAl/S), TRM . 1 0 ( 61 0 .4 ) ) 

RETURN 

END 
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SUHHOUTiNt TLRMS J H3 ♦ B4 t B5 t H6 t B7 , TW t TW2 t I t ARG ) 
DIMENSION Rq ( 1 ) • B5(l)« B6(I)t B7(l) 

DIMENSION TUU)« TW2(1)« AR6(4) 

C 

C 

ARG(l) = BSd-fl) 

ARG(2) = - B4(I<-1) 

ARG<3> s B5(I+1) 

ARG(4) = - <TW(I'»1) * BGCI-d) B7(I + 1)» 

C 

IFII .NE. 1) GO TO 100 ■ 

ARGIll = 0. 

ARG(2> = ARGI2) * B3(I-»1) 

GO TO 200 

100 1F(I .NE. NEOU) GO TO 200 

ARG<2) = AKG(2I * B3<I't’l> 

ARC (3) so, " ■ ■ 

200 continue 

C 

RETURN 

END 


165 



non »-»nnn non 


SUBROUTINE TRID ( SUB .01 AG t SUP iH < N ) 

C *****«***•«** *»*r*»******«*«*****«*«****«**»*«*«»«*« *************** 
C • 

C * 

C * SUBROUTINE TRID SOLVES SYSTEMS OF TRIDIAGONAL LINEAR EQUATIONS 
C • OF THE FORM 

C * SUB<K)*X(K-l)+DlAGIK»*XtK)+SUP(Kl*X(K'fl)=B(K» 

C * BY THE THOMAS ALGORITHM 

C * 

C • DESCRIPTION OF THE ARGUMENTS 

C * ’ 

C • SUB - the N - 1 SUBDIAGONAL COEFFICIENTS 

C * DIAG - THE N DIAGONAL COEFFICIENTS 

C * SUB - THE N - 1 SUPRADIAGONAL COEFFICIENTS 

C » B - INPUT - THE N RIGHT-HAND SIDE CONSTANTS 

C « OUTPUT- THE CALCULATED VALUES OF THE N UNKNOWNS 

c » N - The number of equations 

C » : 

C • I 

C ***»j»««*«*****«****«*«*»*«*««««*«»*««***»«»*#*»«««*«**»**«***«*«**»i 

c 

DIMENSION SUB (N) .DIAG (N) «SUF(N) tB(N) 

IF(N.GT.1»G0 To 10 

solution for N = 1 (A SINGLE EQUATION) 

B(1)=B(1)/DIAG(1) 

RETURN 

ELIMINATION OF TRIDIAGONAL SYSTEM 

0 DO 11 K=2.N 

RATIO = -SUB(K)/DIAG(K-1) 

DIAG(K)=DIAG<K)+RATI0*SUP(I<-1) 

11 B(K)= B(K) + RATI0*B(K-1) 

BACK SUBSTITUTION 

B(NI= B(N)/DIAG(N) 

K=N 

DO 12 1=2, N 

K=K-1 

12 B(K) = (B(K)-SUP(K)*B(K4-1) )/DIAG(K) 

RETURN 

END 
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SUBROUTINE TRIPH 
INCLUUC BASIC 
INCLUDE HEAT 
INCLUDE PROCl 
INCLUDE PR0C2 

DIMENSION SUB(jSl)t DIAG(JS1)« SUP(JS1)« B(JS1I 
DIMENSION ARG(4> 


NEQU = NL - 1 

c ■■ 

DY2 = DY • Dt / 4.0 
C 

DO 100 I s li NEQU 

CALL TERMS (B3tB4«B5tB6fB7«TUtTU2«l«ARG) 

SUB(I) = ARG(l) 

DlAGtl)= ARG(2) ' “ 

SUP(I> = ARGO) 

B( I ) = ARG(4 ) 

C 

100 CONTINUE 
C 

CALL TRlD(SUB*nIAG»SUP.B.NEOU) 

C 

DO 200 1 = It NEQU 
TW2<I + 1) = BCD 

TWKI+l) = TU(I + 1) • (2.0+Bl(Dl)+B2(I + l) ) - TH2(X + 1» - 
$ BKl + l) • TW(l) - TUCl+2) • B2C1 + 1) 

200 CONTINUE 
RETURN 
END 
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SUHROUTINE VIEW (OIAMiNLiStF) 

INCLUOE BASIC 

01 HENS I ON S( JS) «F ( JSl t JSl > tO( JSl t JSl ) 

NLP=NL+1 

Nll=NL-1 

K=DIAH/2. 

J=1 

1 DO 2 I=JiNL 
2s<(SII>-S(J> )/B»**2*2. 

IF (2-1000. > 9.9.10 

9 D(I.J) = (2-SQRT(Z*Z-A.1)/2. "■ 

GO TO 7 

10 0(I.J)=1./Z 

7 D( J.I >=D(I. J) 

2 CONTINUE 
JsJ+1 

IF (J-NU 1.1.3 

3 J=2 

F(2,2)=l,0 

4 DS=ABS(S( J)-S(J-l) )/DIAn*4. 

DO 5 IsJ.NLL 

DSS=ABS(S( 1+1)-S(1) )/0lAH*4. 

FtJ.l+l) = (0(J.T)-D(JVI+n-0(3-l.l)+0(J-l.l + ln/0S 
F(I+l.J)sF( J.l+l»*OS/DSS 

5 CONTINUE 

J=J+1 ■ 

F( J. J)=1.0 

IF (J-NL) 4.6.6 __ 

6 DO 6 1 = 2, NL 

DS=ABS(S(I )-S(I-in/0IAM*4. 
F(I.1)=(0(1.1-1)-0(1.I) )/OS 
F(I.NLP»s(OtNL.I)-D(NL.I-ll)/OS 
DO B Jsl.NLP 

IF (I.EO.J) GO TO 8 _ _ 

F(I,n=F(I.I)-F(I.jr 

8 CONTINUE 

IF(KR7.LE.3)G0 TO 14 
DO 12 1=2. NL 

12 WRITE(6.13) (F( I . J) . Jsl-.NLPJ 

13 FORMAT (8X.10E10.5) 

14 CONTINUE 
RETURN 
END 
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subroutine VIEUl (DIAMC f Dr ) 
iNCLUOr BASIC 
INCLUDE ENLKGY 
INCLUDE HEAT 

COHMON/BLKl/ FL ( JSl , JSl ) •OKI ( JSl ) fORZl JSX) 

DIhENSION FM( JSl) tFQ( JSl) tFPi JSl) «AAN( JSl) 

RR=0IAMC/DT 

PI=3,1416 

DO 40 J=2iNL 

AL=(S( J) -S( J-1) )/DT*2. 

SL=AL**2 
SR=RR^*2 
AH = AL**2 + SR-1 . 

BB=SL-SR+1. 

AN=PI*DIAMC*<S( J)-S( J-1 ) ) 

ANP=PI*DT*(S( J)-S( J-1) ) 

AA=Pl*(0lAMC*»2-DT*»2)/4« 

AAN( J)=AN 

CI = 1 .-1 ./RR+2./(PI*RR)*ATAN<2./AL»S0RT< SR-1, ) ) 

C2=SQRT<4.*SR+SL)/AL 

C3=(4.*(SR-1. ) + <SL/SR)*ISR-2* ) I / ( SL'*^4 . « ( SR-1 • ) ) 

C4=(SR-2, )/SR 

C5=SQRT( ( AB+2. >»»2-4.*SR) 

F( Jf J)=CI-AL/(2.*PI*RR)*(C2*ASIN(C3)-ASIN(C4)4(C2-1. )*PI/2. ) 

FL( Ji J)=l./KR-i./(Pl*RR)»(AC0S(dB/AB)-l./(2*»AL)*(C5*AC0S(BB/(RR*A 
IB) )+B0*ASINa./RR)-PX*AB/2.) ) 

FKM J)=0.5*(l,-F< J% J)-FL( Ji J) ) 

FO( J)=0.5»ANP/AA* (1,-AN/ANP*FL( J* J) ) 

FP( J)=l.-Fn( J)*AAN( J)/AA-FQ( J) 

40 CONTINUE 

DO 41 J=2iNL 
NLP=NL+1 
JJ=J+1 
JM=J-1 

DO 42 K=2fNL 
KP=K-H 

IF(K.E0*J)G0 to 42 
F(JiK)=FM(J)*F M(K)*AAN(K)/AA 
FL< J«K)=FMU)*FQ(K) 

IF(K*GT.JJ)G0 to 43 
IF(K.LT.JH)GO TO 44 
GO TO 42 

43 PP=l.n 
KM=K-1 

DO 4b N=JJ%Kh 

45 PP=P(>*FP(N) 

F( JtK)=F( JfK)*PP 
FL( JtK)=FLC JfK)*PP 
GO TO 42 

44 FP=1.0 

UO 46 N=KP.JM 

46 HP=PH»FP(N) 

F( J.K)rF( J*K)^PP 
FL( J.K)=FL( JfK)*PP 
42 CONTINUE 
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PPsl.O 

1F(J.C0.NL)60 TO 60 
UO 47 N=JJ«NL 
i»7 PP=PP«PP(N) 

60 F( J.NLP»=FHU»*PP 

0H2( J)=FO( J>*PP 

PPsl.O 

IFtJ.C0.2IG0 TO 61 
00 40 NS2.JH 
46 PPsPP*FP(N» 

61 FtJ.DsFMt J)*PP 

ORK JIsFQt J»*PP 

41 CONTINUC 

IF(KR7.tC.3»60 TO 1 
URITC(6«Sli 

51 FORMAT t//lOX,* WALL TO WALL VIEW FACTORS 
00 49 JS2.NL 

49 WRITEt6.50) |F < J.K ) .Ks^NLP » 

50 F0RMAT(10X.10F10.3J ■■ 

URITC(6t53> 

53 FORMAT(//10X« • WALL TO TUBC VIEW FACTORS 
DO 52 JS2.NL 

52 WR1TE(6»50» tFL t J.K > «Ksl «NLP» 

1 CONTINUE _ 

ARAT=D1AMC*01AHC/4.0/DT" 

DO 65 JS2.NL 

20Ss(S(J»-StJ-in*DlAMC/DT 
DO 64 K=2.NL 

64 FL<J.K)=FLtJ.K)*20S/<StK»-S(K-l»> 
ORHJ)sORl<JJ*ARAT/(S(J»-Std-ll » 

65 0 R 2 ( J»sOR2( J|*ARAT/(S(JI-S(J-1I » ■ 

RETURN 

END 


•//) 


•//) 
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APPENDIX A 


BRIEF DEVELOPMENT AND DEMONSTRATION 
OF BIFURCATION APPROXIMATIONS TO DIFFUSION COEFFICIENTS 


Because the development of the set of transport relations used in 
this code is not in the readily available combustion literature, it is 
appropriate to include a brief summary of the developinent in this report. 

The formulation derives much of its utility from the simplification 
introduced when binary diffusion coefficients are approximated with the 
bifurcation relations. These are 



(A-1) 


where the F^ are diffusion factors for individual species which are independent 
of the system of species involved. They could be temperature dependent but 
this flexibility is rarely justified. The accuracy of the approximation 
has been tested for a variety of systems. Two examples taken from Reference 3 
are shown in Tables A-1 and A-2. These results are typical when good (and 
consistent) sets of diffusion data are correlated. Note that these particular 
correlations have arbitrarily taken F« as 1.0 to anchor an otherwise floating 
correlation. 

In order to develop an explicit flux relationship using this .approximation, 
it is appropriate to start with the Stefan-Maxwell relation 


9s " Z 


P^ij 




/ (^) . (^ ) 




s- 


(A-2) 
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TABLE A-1. CORRELATION OF BINARY DIFFUSION COEFFICIENTS FOR A 
HYDROGEN-OXYGEN SYSTEM USING PRESENT METHOD 

TEMPERATURE = 12,000®R, PRESSURE = 1 ATM 


Species 

■* j- 

Pij From 
Kinetic Theory 

(ft^/sec) X 100 

>^i 

Pij 

ent 

(ft^ 

From Pres- 
Correlation 

/sec) X 100 

Error Using 
Present 
Correlation 
(Percent) 

Error If All 
Oij Are As- 
sumed Equal 
(Percent) 


( a ) 

Di ffusion 

coefficients calculated using Lennard-Jones 



potential 

with 

force data 

from Svehla (Ref. 

14) 

H H2 

36.0260 

0.24713 


53.1613 

47.6 

- 63.1 

H H20 

25.9891 



23.7639 

- 8.6 

- 43.9 

H 0 

26.6238 



24.7360 

- 7.1 

- 50.1 

H 02 

22.8038 



19.7757 

-13.3 

- 41.7 

H OH 

26.4341 



24.3147 

- 8.0 

- 49.7 

H2 H20 

17.3862 

0.3720 


15.7877 

- 9.2 

- 23.5 

H2 0 

17.7166 



16.4335 

- 7.2 

- 24.9 

H2 02 

15.0085 



13.1381 

-12.5 

- 11.4 

H2 OH 

17.5759 



16.1537 

- 8.1 

- 24.4 

H20 0 

7.0928 

0.8322 


7.3461 

3.6 

87.4 

H20 02 

5.2795 



5.8730 

11.2 

151 .6 

H20 OH 

6.9078 



7.2210 

4.5 

92.4 

0 02 

5.6458 

0.7995 


6.1132 

8.3 

135.4 

0 OH 

7.2643 



7.5163 

3.5 

82.9 

02 OH 

5.4946 

1.0000 


6.0091 

9.4 

141 .9 

OH 


0.8133 





Average 

Absolute Error 




10.8 

68.6 


( b ) 

Diffusion 

coefficients calculated using values for 



collision 

cross 

.-sections suggested by Svehla 

(Ref. 15) 

H H2 

67.6000 

0.2208 


74.4024 

10.1 

- 77.1 

H H20 

28.3200 



27.0030 

- 4.7 

- 4t.4 

H 0 

27.7200 



30.8482 

11.3 

■ 1 

H 02 

24.5500 



22.5734 

- 8.1 

- 37.1 1 

H OH 

29.5900 



27.5549 

- 6.9 

- 47.8 1 

H2 H20 

19.5800 

0.3034 


19.6568 

0.4 

- 21 .1 

H2 0 

23.6000 



22.4560 

- 4.8 

- 34.5 

H2 02 

17.1900 



16.4323 

- 4.4 

- 10.1 

H2 OH 

20.1600 



20.0586 

- 0.5 

- 23.3 

H20 0 

8.2950 

0.8360 


8.1500 

- 1.7 

86.3 

H20 02 

5.7150 



5.9638 

4.4 

170.4 

H20 OH 

7.1450 



7.2799 

1.9 

116.3 

0 02 

6.8500 

0.7317 


6.8131 

- 0.5 

125.6 

0 OH 

8.6060 



8.3166 

- 3.4 

79.6 

02 OH 

5.5520 

1.0000 


6.0857 

9.6 

178.3 

OH 


0.8192 





Average 

Absolute Error 




4.8 

73.1 
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TABLE A-2. CORRELATION OF BINARY DIFFUSION COEFFICIENTS FOR AN 
OXYGEN-NITROGEN-CARBON-HYDROGEN SYSTEM BASED ON DATA 
OF SVEHLA (Refs. 14 and 15) 

TEMPERATURE = 12,000®R, PRESSURE = 1 ATM 




D-jj From 


Pij From Pres- 

Error Using 

Error If All 

Species 

Kinetic Theory 


ent Correlation 

Present 

Pij Are As- 

i 

j 

(ft^/sec) X 100 


(ft^/sec) X 100 

Correlation 
i Percent) 

sumed Equal 
(Percent) 

0 

02 

6.8500 

0.7393 

6.0528 

-11.6 

- 4.5 

0 

N 

7.3372 


7.6554 

4.3 

- 10.8 

0 

N2 

5.3995 


5.6277 

4.2 

21.1 

0 

CO 

5.4662 


5.6846 

4.0 

19.7 

0 

C02 

4.4638 


4.6277 

3.7 

46.5 

0 

C 

8.0754 


8.3865 

3.9 

- 19.0 

0 

C3 

5.1820 


5.3597 

3.4 

26.2 

0 

CN 

5.3620 


5.5958 

4.4 

22.0 

0 

H 

27.7200 


29.8130 

7.6 

- 76.3 

0 

H2 

23.6000 


20.4311 

-13.4 

- 72.2 

0 

H20 

8.2950 


7.5057 

- 9.5 

- 21.1 

0 

OH 

8.6060 


7.7923 

- 9.5 

- 24.0 

0 

CH4 

5.8190 


6.0848 

4.6 

12.4 

0 

C2H 

4.8947 


5.0977 

4.1 

33.6 

0 

HCN 

4.8625 


5.0401 

3.7 

34.5 

02 

N 

5.6566 

1.0000 

5.6595 

0.1 

15.6 

02 

N2 

3.9611 


4.1604 

5.0 

65.1 

02 

CO 

4.0028 


4.2025 

5,0 

63.4 

02 

C02 

3.1637 


3.4212 

8.1 

106.7 

02 

C 

6.3129 


6.2000 

- 1.8 

3.6 

02 

C3 

3.7100 


3.9624 

6.8 

76.3 

02 

CN 

3.9623 


4.1369 

4.4 

65.1 

02 

H 

24.5500 


22.0402 

-10.2 

- 73.5 

02 

H2 

17.1900 


15.1043 

-12.1 

- 62.0 

02 

H20 

5.7150 


5.5489 

- 2.9 

14.5 

02 

OH 

5.5520 


5.7607 

3.8 

17.8 

02 

CH4 

4.4735 


4.4984 

0.6 

46.2 

02 

C2H 

3.6310 


3.7686 

3.8 

80.1 

02 

HCN 

3.5678 


3.7260 

4.4 

83.3 

N 

N2 

5.4277 

0.7907 

5.2620 

- 3.1 

20.5 

N 

CO 

5.4763 


5.3153 

- 2.9 

19.4 

N 

C02 

4.5136 


4.3270 

- 4.1 

44.9 

N 

C 

7.9727 


7.8416 

- 1.6 

- 17.9 

N 

C3 

5.2069 


5.0115 

- 3.8 

25.6 

C 

CN 

5.3784 


5.2323 

- 2.7 

21.6 

N 

H 

25.5139 


27.8700 

9.3 

- 74.4 

N 

H2 

17.1218 


19.1036 

11.6 

- 61.8 

N 

H20 

6.9743 


7.0181 

0.6 

- 6.2 

N 

OH 

7.1732 


7.2861 

1.6 

- 8.8 

N 

CH4 

5.7836 


5.6895 

- 1.6 

13.1 

, N 

C2H 

4.9083 


4.7665 

- 2.9 

33.3 

N 

HCN 

4.8645 


4.7126 

- 3.1 

34.5 
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TABLE A-2. Continued 


Species 
i j 

Pij From 
Kinetic Theory 

(fAsec) X 100 


Pij From Pres- 
ent Correlation 

(ft^/sec) X 100 

Error Using 
Present 
Correlation 
(Percent) 

Error If All 
Pij Are As- 
sumed Equal 
(Percent) 

N2 

CO 

3.8943 

1.0756 

3.9074 

0.3 

68.0 

N2 

C02 

3.1114 


3.1809 

2.2 

110.2 

N2 

C 

6.0528 


5.7645 

- 4.8 

8.1 

N2 

C3 

3.6214 


3.6840 

1.7 

80.6 

N2 

CN 

3.8603 


3.8463 

- 0.4 

69.4 

N2 

H 

21.3750 


20.4922 

- 4.1 

- 69.4 

N2 

H2 

14.1671 


14.0435 

- 0.9 

- 53.8 

N2 

H20 

5.0300 


5.1591 

2.6 

30.0 

N2 

OH 

5.2629 


5.3561 

1.8 

24.3 

N2 

CH4 

4.3182 - 


4.1824 

- 3.1 

51.5 

N2 

C2H 

3.5367 


3.5039 

- 0.9 

84.9 

N2 

HCN 

3.4655 


3.4643 

- 0.0 

84.8 

CO 

C02 

3.1390 

1.0647 

3.2131 

2.4 

99.2 

CO 

C 

6.1194 


5.8229 

- 4.8 

9.9 

CO 

C3 

3.6584 


3.7213 

1.7 

72.0 

CO 

CN 

3.8938 


3.8853 

- 0.2 

64.8 

CO 

H 

21.6122 


20.6996 

- 4.2 

- 66.9 

CO 

H2 

14.2296 


14.1856 

- 0.3 

- 50.4 

CO 

H20 

5.1001 


5.2113 

2.2 

27.1 

CO 

OH 

5.3305 


5.4103 

1.5 

22.3 

CO 

CH4 

4.3595 


4.2248 

- 3.1 

51.5 

CO 

C2H 

3.5680 


3.5394 

- 0.8 

80.9 

CO 

HCN 

3.5023 


3.4994 

- 0.1 

82.9 

C02 

C 

4.9902 

1.3079 

4.7402 

- 5.0 

35.0 

C02 

C3 

2.8753 


3.0294 

5.4 

111.3 

C02 

CN 

3.1245 


3.1629 

1.2 

102.4 

C02 

H 

18.4881 


16.8510 

- 8.9 

- 59.4 

C02 

H2 

12.2917 


11.5481 

- 6.0 

- 39.1 

C02 

H20 

4.1217 


4.2424 

2.9 

58.7 

C02 

OH 

4.3441 


4.4044 

1.4 

50.6 

C02 

CH4 

3.5835 


3.4393 

- 4.0 

82.5 

C02 

C2H 

2.8685 


2.8813 

0.4 

128.0 

C02 

HCN 

2.7965 


2.8488 

1.9 

133.9 

C 

C3 

5.7767 

0.7218 

5.4901 

- 5.0 

13.2 

C 

CN 

6.0033 


5.7319 

- 4.5 

8.9 

C 

H 

26.1719 


30.5380 

16.7 

- 75.0 

C 

H2 

18.0635 


20.9280 

15.9 

- 63.8 

C 

H20 

7.5630 


7.6883 

1.7 

- 13.5 

C 

OH 

7.9334 


7.9818 

0.6 

- 17.5 

C 

CH4 

6.3330 


6.2328 

- 1.6 

3.3 

C 

C2H 

5.3831 


5.2217 

- 3.0 

21.5 

C 

HCN 

5.3406 


5.1626 

- 3.3 

22.5 
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TABLE A-2. Concluded 


Species 
i j 

Pij From 
Kinetic Theory 

(ft^/sec) X 100 

»^i 

From Pres- 
ent Correlation 

(ft^/sec) X 100 

Error Using 
Present 
Correlation 
(Percent) 

Error If All 
Pij Are As- 
sumed Equal 
(Percent) 

C3 

CN 

3.6276 

1.1293 

3.6632 

1.0 

80.3 

C3 

H 

21.0069 


19.5166 

- 7.1 

- 68.9 

C3 

H2 

13.9792 


13.3749 

- 4.3 

- 53.2 

C3 

H20 

4.8271 


4.9135 

1.8 

35.5 

C3 

OH 

5.0416 


5.1011 

1.2 

29.7 

C3 

CH4 

4.1210 


3.9833 

- 3.3 

58.7 

C3 

C2H 

3.3265 


3.3371 

0.3 

96.6 

C3 

HCN 

3.2583 


3.2994 

1.3 

100.7 

CN 

H 

20.9403 

1.0817 

20.3763 

- 2.7 

- 68.8 

CN 

H2 

13.8853 


13.9641 

0.6 

- 52.9 

CN 

H20 

4.9948 


5.1300 

2.7 

30.9 

CN 

OH 

5.2299 


5.3259 

1.8 

25.1 

CN 

CH4 

4.2953 


4.1588 

- 3.2 

52.3 

CN 

C2H 

3.5330 


3.4841 

- 1.4 

85.1 

CN 

HCN 

3.4626 


3.4447 

- 0.5 

88.9 

H 

H2 

67.6000 

0.2030 

74.3967 

10.1 

- 90.4 

H 

H20 

28.3200 


27.3309 

- 3.5 

- 76.9 

H 

OH 

29.5900 


28.3745 

- 4.1 

- 78.0 

H 

CH4 

20.3467 


22.1568 

8.9 

- 68.0 

H 

C2H 

18.6611 


18.5624 

- 0.5 

- 64.9 

H 

HCN 

18.8560 


18.3525 

- 2.7 

- 65.3 

H2 

H20 

19.5800 

0.2963 

18.7301 

- 4.3 

- 66.7 

H2 

OH 

20.1600 


19.4453 

- 3.5 

- 67.6 

H2 

CH4 

13.7590 


15.1843 

10.4 

- 52.4 

H2 

C2H 

12.5045 


12.7210 

1.7 

- 47.7 

H2 

HCN 

12.5953 


12.5771 

- 0.1 

- 48.1 

H20 

OH 

7.1450 

0.8064 

7.1436 

- 0.0 

- 8.5 

H20 

CH4 

5.4665 


5.5782 

2.0 

19.6 

H20 

C2H 

4.5559 


4.6733 

2.6 

43.6 

H20 

HCN 

4.5242 


4.6205 

2.1 

44.6 

OH 

CH4 

5.6987 

0.7767 

5.7912 

1.6 

14.8 

OH 

C2H 

4.7817 


4.8517 

1.5 

36.8 

OH 

HCN 

4.7388 


4.7969 

1.2 

38.0 

CH4 

C2H 

3.9244 

0.9948 

3.7886 

- 3.5 

66.7 

CH4 

HCN 

3.8677 


3.7457 

- 3.2 

69.1 

C2H 

HCN 

3.1729 

1.1874 

3.1381 

- 1.1 

106.1 

HCN 



1.2009 




Average 

Absolute Error 



3.7 

50.9 
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where is the mole fraction of species i, T is the tem)jerature, is the 
binary diffusion coefficient for species i and j, and is the multi- 
component thermal diffusion coefficient for species i. Substituting the 
approximation for binary diffusion coefficients (Eq. A-1) into the Stefan- 
Maxwell relation (Eq. A-2) and rewriting in terms of mass fractions, yields 



where, for convenience, a total diffusion mass flux has been defined as 
the sum of the molecular and thermal diffusional fluxes. 


'^i = 




T p £n T \ 
i ^3S } 


(A-4) 


Multiplying each side of Eq. (A-3) by summing over all i, and 

noting that the sum of the diffusive fluxes is zero and the sum of the mass 
fractions is unity yields: 


I fi!i . 00 1 . pD I «i 

J 1 T7 3s - j Tj 3S 

Substituting Eq. (A-5) into Eq. (A-3) results in 


(A-5) 


^ y "jiiL Vj 

3S - M, f TJ 3S 00 M, f T- 


(A-6) 


At this point is is convenient to define several new quantities. 






(A-7) 
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U2 = 4. 
0 


(A-8) 

j 

(Y^/Fj^) (dFj/dT) 

(A- 9 ) 

Taking the derivitor of Eq. A-8 yields 



ds ” zL F . ds 

w 

^ M.X.- dF. 
M 

(A- 10 ) 


Introducing Equations A-7 through A-lO into A-6 yields after some rearranging 




+ D 


9T ^ 

■5? ■ 


= .£L_ 


y^Fi 


'^^i . ^i fdM p ^^2 p ,, dTll 

Hi“ TT [hF ■ BT" ■ Isl 


The variation of with temperature as determined with a nine component 
system over a range from 4000° to 16000°R in Reference 3 were found to rarely 
exceed 0.1%. Consequently has been taken as zero and a universal set of 
F^. determined for all species. Correlations of these values with molecular 
weight have been reasonably good as indicated in Figure A-1. The correlating 
equation 





Mi 

26 


0.461 


(A-12) 


is recommended when specific values of F^. are not available from other correlations. 

It is apparent that 15 must represent the temperature and pressure dependence of 

the P. j. Although ^ need have no specific relation to a real diffusion coefficient, 
* \) 

the arbitrary choice of Fq as 1.0 prompts the interpretation of U as the self 


177 



Molecular Weight, 


Figure A-1. Correlation of diffusion factors with molecular weight 










diffusion coefficient of O 2 , Pq q . This would force a precise fit of 

the correlation to P« « but would not in general provide the best overall 

^Z 

correlation. Ignoring this caveat we can evaluate « from the equation 

^ 2.^2 

of Reference 4 as given by Equation 7 of the text: 


U = 2.628 X 10 


-3 


T(T/Mref)^^^ 

Po2 

rerij 


(cm'^/sec) 


(A-13) 


with T in °K, P in atmospheres, and collision cross section, o in A. For 

O 

O 2 as the reference species, o is equal to 3.467 A. Using the data from 
Reference 4, Figure A-2 demonstrates that the integral expression for transport 
properties is approximated by: 


ij 


1.07 [T/(c/k)]’°-^^^ 


where the maximum energy of attraction function, e/k, for O 2 is 106.7 and 
thus : 


D = 0.172 X 10"^ t^-659/p (cm^/gec) 


(A-14) 


The extension of this correlation to thermal conductivity is fully 
developed in the text. However, the system viscosity is also required in 
PROF for evaluating Reynold's number correlations. Use is made of the 
correlation suggested by Buddenberg and Wilke (Reference 16) and endorsed 
by Hirschfelder et al. (Reference 4), namely 

^i^i 

^Mix ^ 




RTu,. 
385 pjif 


1 j ^ij 


(A-15) 
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By introducing the bifurcation relations, taking from Reference 4 the relation 
for pure component viscosity 


”i " Pi^^ii* (A- 16 ) 

assuming - 1.12 (actually varies from 1.10 to 1.14 in the temperature 
range of interest), and adjusting 1.385 to 1.344 for simplification, there 
is obtained 



(A-17) 


The results presented in this Appendix show that the bifurcation approximation 
can result in major simplification in transport property evaluations with only 
minimal loss in accuracy. 
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